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ABSTRACT 


In order to assist administrators of long-term 
facilities in responding to the changing environment 
of and demand for institutional care of the elderly 
and chronically infirm, a unified approach to patient 
classification and nurse staffing has been developed. 
The approach is essentially structure oriented and is 
directed toward the presentation of alternative nurse 
staffing strategies; its basis is the quantification of 
patient nursing needs through intermediate steps of 
assessment and classification. The classification system, 
in turn, has been related to the demand for nursing 
services, aS categorized by care area, for each of 
three homogeneous patient groupings. Effective allo- 
cation and assignment of patient nursing care activities 
have been modeled by use of techniques of mathematical 
programming. 

The ability to predict, with reasonable confi- 
dence, the level of care to which a patient would 
most appropriately be assigned based on an assessment 
of various functioning status items, psycho-social 
Status indicators, and medically defined conditions 
was considered initially. Application of an existing 
nonlinear multiple regression technique to patient 


assessment data drawn from a comprehenSive assessment 


i 





alt 


instrument yielded a means by which patients could be 
classified according to 37 health status indicators. 
The regression method was chosen for its specific 
applicability to problems of fitting a polychotomous, 
ordered response, analogous to the level of care. In 
order to achieve a classification system of reasonable 
size for use in typical long-term care facilities, 
various subsets of the original 37 variables were 
tested for their prediction and recognition powers; 
a subset of 12 variables was shown to possess adequate 
capabilities in this regard. This subset was synthesized 
into an implementable patient classification system 
procedure. 

The determination of an optimal staff mix,. 
allocation of nursing time, and assignment of nurses 
to patient care demands was captured in the Basic 
Staffing Model (BSM), a mixed-integer linear program. 
The BSM maximizes an objective function that combines 
the concepts of the appropriate assignment of nursing 
personnel to specific patient care activities and 
the satisfaction of high priority demands at greater 
than minimum levels, if possible. Model constraints 
include the availability of nursing resources, budgetary 
limitations, legal staffing requirements, and adherence 
to bounded representations of patient nursing care 


activities. Two extensions to the BSM, Model I and 
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Model II, were developed to facilitate the exploration 
of alternative staffing policies. Model I allows for 
the parametric representation of changes in the total 
patient-centered services provided and personnel budget, 
while Model II enables parametric alterations in upper 
and lower bounds on staff and patient demands to be 
examined. 

Two algorithms were derived for the efficient 
solution of both BSM extensions. Recognizing the forms 
of the models, branch and bound methods incorporating 
post-optimality analysis in specialized linear program- 
ming routines were developed. These methods were based 
on the Upper Bounded Dual Simplex algorithm of Wagner; 
their computational characteristics for both general 
parametric mixed-integer linear programs and specific 
examples of applications to nurse staffing problems 


were investigated. 
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Chapter I. Introduction: The Need for Research 


I.l General Discussion 

As the twenty-first century approaches, one of 
the most important components of our multi-billion 
dollar a year national health system, long-term care, 
seems certain to become even more vital. With age being 
the primary, although not exclusive, determinant of entry 
into the long-term care system, demographic projections 
indicate that potential demands will be high. Brody 
[12] notes that while at the turn of the century some 
three million persons (4% of the U.S. population) were 
age 65 or over, we now have approximately 20 million 
persons (10%) in this category. Figures cited by Flagle 
[24] show a 55% increase by 1990 in persons 65 or over, 
based on 1970 census data for the state of Maryland. 
The United Nations estimates that for the period 1970- 
1980 the world's aged population will grow by 38% while 
the total population increases by 28%. 

The response of Federal and state governments to 
this trend has for the most part been in the form of 
specific provisions in the current Medicare/Medicaid 


legislation.* These programs provide support or 


* Titles XVIII and XIX of the Social Security Act. 
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supplementary assistance to those persons in need of 
long-term care and who meet certain eligibility criteria. 
Still being debated is proposed legislation establishing 
a National Health Insurance (NHI) system which presum- 
ably will contain payment structures for some forms of 
care of the chronically ill. Flagle [24] notes that the 
five-year period following the installation of Medicare 
saw a 68% increase in nursing home patients. If this is 
any indication of the underlying elasticity of demand for 
long-term care service, the potential impact of NHI could 
be highly significant, especially if payments are weighted 
towards institutional care rather than possible alterna- 
tive care programs. The relative magnitude of this impact 
can be estimated by noting figures prepared by the Depart- 
ment of Health, Education, and Welfare (HEW) [75] showing 
that there are currently some 16,000 long-term care 
Facilities in the U.S. caring for about 1.1 million per- 
sons, representing approximately 5.2% of the population 
over 65. 

Because of the large amount of federal funds 
being devoted to health care in general, and with the 
probable expenditure of much larger sums in the future, 
the need for accountability and cost-benefit analysis has 
become evident. The legislative reaction has been the 


creation of Professional Standards Review Organizations 





(PSRO) on a state-wide basis to set standards and cri- 
teria for the various component specialties of the 
health delivery system, and to monitor compliance with 
such guidelines. The enabling legislation* specifically 
calls for such standards review organizations in the 
long-term care area, although to date little progress 
has been made towards complying with this mandate. It 
is certain, however, that the decisions ultimately made 
will bear heavily on both the forms of long-term care 
delivery and the quality of care provided. 

The potential impact of population increases, cur- 
rent and proposed publicly-supported health insurance 
systems, and PSRO taken both individually and in concert 
have given rise to problems that must be dealt with by 
national and regional policy-makers. It is incumbent 
upon them to predict the resource requirements generated 
by these factors in a number of alternative care settings 
and to choose the mix of services which will most ade- 
quately serve the client population. Assessment of the 
social and technological requirements for prolonging 
patient independence in home care settings must be carried 
out and weighed against the alternative of institution- 


alization. Secondary effects such as those to be felt 


* 1972 Amendments to the Social Security Act. 





in the labor and education markets for health profes- 
sionals must be recognized. Where increased utilization 
of institutions is projected, plans must be made to in- 
crease existing capacity either by direct government 
involvement or through incentives to the private sector. 
In short, there is much work to be done on the "macro" 


planning level before the predictions become reality. 


1.2 The Need for Research at the Facility Level. 

It should be obvious at this point that regard- 
less of the increasing emphasis on alternative settings 
for long-term care, the future portends greater require- 
ments for quantity and quality of care to be provided by 
facilities. Being licensees of their respective state 
governments, with additional federal certification re- 
quired in order to care for Medicare patients, individu- 
al facilities must be capable of reacting in a positive 
way to the requirements set forth by requlatory agencies. 
In doing so, they will in a sense be translating to a 
"micro" level those broad objectives mentioned earlier. 
This dissertation represents an effort towards assisting 
the individual facility in this task. 

As an historical note the reliance upon an 
institutional setting for care of the aged and chronical- 


ly infirm is a relatively recent phenomenon. Brody [12] 





identifies six distinct periods in the history of long- 
term care in the U.S. The colonial period and the period 
of the 19th and early 20th centuries to 1920 were charac- 
terized by care provided in the home, with relatively 
little dependence on locally-financed public almshouses. 
Between 1920 and 1935 the rise of the eleomosynary insti- 
tution and the creation of publicly-financed care systems 
took place, after which the years 1935-1945 constituted 

a time of reaction to such developments. The postwar 
years to 1965 saw the increased entry of privately-owned 
proprietary nursing homes into the LTC system, while con- 
currently, the major mechanisms for public and private 
support were being molded. The advent of the current 
period is marked by the enactment of Titles XVIII (Medi- 
care) and XIX (Medicaid) of the Social Security Act of 
1965, from which the concepts of a stratified LTC system 
responsive to patient needs are derived. 

Since we shall later have cause to quantify the 
relationships between patient health status and the 
levels of care of the stratified system, we pause to 
clarify some definitions. Because most of this study 
was conducted in the State of Maryland, definitions of 
LTC facilities are provided below as taken from that 
State's Report of the Governor's Commission on Nursing 


Homes [56, p. 117]; roughly analogous interpretations 





are given in other sources (cf. Coggeshall [16]). Note 

that to a large extent, the levels of care correspond 

to the type of facility involved, and we shall therefore 

use the terminology interchangeably.* 

Chronic Disease Hospital: A facility that provides 
medical care to persons incapacitated by illness, 
with the intent and, indeed, mandate to rehabili- 
tate and return them to the community. 

Skilled Nursing Home: A facility that maintains facili- 
ties and staff necessary to render skilled nurs- 
ing care (which means by an RN or by an LPN under 
supervision). 

Intermediate A Facility: A facility that provides long- 
term care to residents whose illness is not acute 
and whose proper care requires no more than 8 
hours per day under supervision of an RN or LPN. 

Intermediate B Facility: A facility that provides 
services which a person normally provides him- 
self, but for which he is now dependent upon 
others because of advanced age, physical, or 
mental limitations. 


Current practice calls for the physical placement of a 


* More detailed discussion of specific regulations per- 
taining to levels of care will be provided later. 





long-term care patient into a facility certified to pro- 
vide that level of care which is deemed appropriate for 
his or her needs. A change in patient requirements thus 
necessitates movement from one facility to another, 
although recent trends indicate an increase in facilities 
with multi-level certifications, thereby eliminating this 
requirement. 

With this background information in mind, we must 
necessarily consider the task of the facility administra- 
tor, upon whom the responsibility for compliance with 
regulatory agencies rests. In broad terms, the adminis- 
trator must be concerned with the adequate and timely pro- 
vision of the three components of institutional life as 
identified by Brody [12]: 

l. Basic maintenance services (shelter, food, 

sanitation); 

2. Medical and paramedical services to foster 
maximum physical health and functioning 
capacities; 

3. Psychosocial components. 

The second and third of these components are highly 
labor-intensive, and their provision accounts for the 
major portion of LTC facility operating costs. Estimates 
provided by the Governor's Commission [56], Schwartz 


[63], and HEW [75] put total labor costs at about 





55%-60% of total operating budgets. The above studies 


indicate that payroll and fringe benefit costs attribut- 


able to nurse staffing make up anywhere from 35% to 50% 


of operating costs. In short, the administrator is con- 


fronted with a budgeting problem which is certain to in- 


volve relatively large sums to be spent on nursing staff. 


Cost, however, is neither the only nor the most 


important consideration in developing a nurse staffing 


program. 


Aydelotte [5] points out five additional fac- 


tors which must be weighed and acted upon in order that 


an acceptable program may be obtained: 


Item 


(2) 


I 


Quality of patient care to be delivered and 
its measurement; 

Characteristics of the patients and their 
care requirements; 

Prediction of the supply of nursing manpower 
required for (1) and (2); 

Logistics of the staffing program and its 
contro" 

Evaluation of the quality of care desired, 
thereby measuring the success of the staffing 


itself. 


is perhaps the key factor among all five, since 


unless a reasonably accurate appraisal of the care needs 


of facility residents can be made, the staffing mix and 





staffing patterns developed in concert with budget con- 
Siderations will be at least partially ineffective, with 
a consequent degradation in quality of care. 

What means, then, are available to administrators 
and their directors of nursing for relating patient 
health status and the implied care needs to a nurse staff- 
ing system? Relying upon the placement status of the 
patient as an assessment tool allows the use of minimum 
staffing levels as specified by state regulations, for 
example, as a guide. The number of hours of bedside care 
per patient of each level, staff-to-patient ratios, and 
Supervisory personnel required could be determined for the 
patient mix that exists, and the minimum or near minimum 
resources maintained. If further funds are available, 
some choice could be made as to additional personnel re- 
quired. Patient placement status, however, may not neces- 
sarily be an accurate reflection of demand for nursing 
services, since many other intervening factors are in- 
volved in the placement decision. Information generated 
in the course of utilization review might give a more 
accurate picture of patient care needs, but relating such 
information to nursing requirements and hence to budgeting 
is likely to be difficult and to occur in an untimely 
fashion. 


The need for maximization of care quality for a 
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given budgetary level, and the relationship between 
patient health status and nurse staffing are thus central 
issues confronting the administrator. The evaluation of 
quality of care suggested in item (5) above is a most 
difficult task. What can be done at present is to either 
assume quality is acceptable with the staffing pattern 
employed or to rely on some appropriate measure of quali- 
ty surrogates. In the former case, a staffing pattern 
which contributes to job satisfaction by assigning nurses 
to tasks commensurate with their training can lead to 
better staff morale. Revans [57,58] has demonstrated 
that staff morale is a good indicator of the quality of 
care provided. In connection with the use of quality of 
care surrogates, Burroughs [14] showed that a number of 
eaSily measured indicators (e.g., number of requests for 
pain medication, number of incontinent patients, length 
of stay, etc.) could be quantitatively related to subjec- 
tive estimates of quality of care by use of discriminant 
analysis. These studies were conducted in acute care 
facilities, however, and additional analysis is needed 


in the LTC area. 


1.3 The Proposed Research 
In order to provide administrators and nursing 


directors, in their roles as decision makers, with useful 





der 


tools for assessing the relative merit of alternative 
Caursessof acwiıon, Bkunlfied appsoacheiss pmoposed. This 
methodology is based on a multi-disciplinary operations 
research point of view and will address the following 
four issues: 

l. Budgetary restrictions; 

2. Assessment of patient care needs; 

3. Nurse staff allocation; and 

4. Quality of care. 

It will be shown that a key factor in such an approach is 
relating patient health status, through classification, 
with nursing care requirements. 

After an initial review of the existing litera- 
ture pertaining to the general problems of patient classi- 
fication and nurse staffing, a classification system for 
long-term care patients into levels of care will be de- 
veloped. The system will initially take as its input in- 
formation available from a comprehensive patient assess- 
ment instrument. It will be shown that a reasonably 
accurate level of care classification can be derived from 
components of the assessment instrument by application of 
a regression-based statistical technique. A proposed 
classification instrument designed for easy use by facili- 
ty staff will be demonstrated. 


We shall then forge the link between patient 
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assessment and classification, and demand for nursing 
services, using information thereby obtained as inputs 
into a mathematical programming model designed to syn- 
thesize the four elements mentioned above. Certain theo- 
retical properties of the model form will be advanced to 
better allow for examination of alternatives. Such ad- 
vances are made in the area of parametric mixed-integer 
linear programming. 

The proposed staffing model takes as its objec- 
tive to be maximized the allocation of nursing care activi- 
ties to the most appropriate nursing skill level avail- 
able and a priority measure on care area/classification com- 
binations. By invoking such a discipline, good quality care 
should be achievable. In order to assure that the model 
is computationally tractable, a grouping of the nursing 
tasks into various "care areas" will be utilized. Atten- 
tion is given to the varying needs of patients for the 
tasks of each care area by basing the model on a subdivi- 
sion of the care areas by the levels of care indicated 
through patient classification. Constraints on personnel 
budget and availability, legal staffing minimums, and 
unit aggregate demands by care area and patient classifi- 
cation are also included. In addition, the latter of 
these three constraint sets includes a priority related 


system for meeting demands. 
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Through the use of mathematical programming 
methods, staffing levels and skill level mix as well as 
allocatıon of direct care demands by skill level will be 
examined under a variety of circumstances. An example 
of how such information may be generated and presented 
will be provided. Finally, we will conclude by noting 
difficulties encountered in the course of the research, 
conclusions reached, and recommendations for further 


study. 





Chapter II. A Review of Existing Methodologies for 
Patient Classification and Nursing 


Resource Allocation 


E rodauac Lion 

In order to place the current study into proper 
perspective, we briefly review some of the major studies 
most closely related to our twin objectives of patient 
classification and nursing resource allocation. The sum- 
maries presented are in no way meant to be an exhaustive 
survey of preceding work; rather, they represent an 
effort to capture the essence of the extensive body of 
significant literature in the field. 

Studies involving classification, and the con- 
comitant nurse staffing problems of allocation and assign- 
ment in an acute care setting will be discussed first, 
since historically they provide the conceptual basis for 
subsequent work in long-term care. An examination of 
work specifically devoted to this latter area will follow. 
Further references and amplification of earlier research 
may be found in the surveys conducted by Aydelotte [5], 


Stimson and Stimson [71], and Shuman and Wolfe [67]. 
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II.2 Patient Classification in an Acute Care Setting 
Connor [18] 

It can be argued that most patient classification 
procedures in use today are based on the fundamental con- 
cepts proposed by Connor in his pioneering studies of 
nursing unit staffing requirements. Essentially, the aim 
of his study was to determine the underlying functional 
relationship between patient care needs and the response 
to these needs in terms of nursing manpower. If patients 
can be assessed and their care needs anticipated, a more 
flexible staffing pattern is possible than the commonly 
used ratios based primarily on census counts. As an ini- 
tial step the amount of direct care provided to a sample 
patient population was measured. It was found that the 
time devoted to direct care for individual patients by 
the nursing staff varied with the degree of self- 
sufficiency of the patient, i.e., the more self-sufficient 
the patient, the less direct care he received. Moreover, 
it was found that the patient population could be classi- 
fied into essentially three groups: 1) a self-care group 
at one extreme consisting of patients who were able to 
take care of most of their needs and made little or no 
demand on nursing staff, 2) a total-care group at the 
other extreme consisting of patients who were often seri- 


ously ill or bedfast and made large demands on nursing 
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staff for such care activities as feeding, bathing, etc., 
and 3) an intermediate care group who required only par- 
tial assistance from the nursing staff. The three group- 
ings were verified by the emergence of a tri-modal dis- 
tribution of measured direct care times, as shown in 
Figure II.1 

A classification scheme was developed, using a 
Boolean logic, that classified patients into self, par- 
tial, and total care categories on the basis of those 
factors that reflected their degree of self-sufficiency. 
Indeed, 1t was these factors that correlated with the 
large demands made on nursing staff in the provision of 
direct care time. The basic classification scheme, and 
the factors involved, is shown in Figure II.2. 

Connor then determined the average number of 
hours of direct care devoted to patients in each of the 
three patient categories by noting the group averages 
for the three populations as shown in Figure II.l. After 
determining the average number of hours devoted to in- 
direct care by additional work sampling on typical 
medical-surgical units, Connor proposed that the total 
nursing hours for a given patient mix on an acute care 
unit could be predicated by first computing a Direct Care 
Index 
(11-1) 
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Figure II.2 


Connor's Combinations o£ Factors for 
Categorization of Patients 


Caregery f. Dolf-sora 


Any of the following combinations checkad 


(a) 


(b 


{cj 


Ambuletory, of up’ ta chalr—Self {without 
assistance] 

Feading self, or requiras food cut 

Bsthing In bathroom, or at bedslde—Peartlal 
Self {con batha self excspt for back end per- 
hops extromitles) 

Ambulatory—wlth assisionce 

Up In cholr—Self 

Bathing In bathroom, or at bedslde—Parilal 
Self 

As In ta and b} wlith 

Vislon Inadequols 

Oxygen theropy 

LV. feeding 

buf no two of thasa factors slmultonaously 


Catogocry If, Parttet er Inrermedtato core 


Any of the following combinations checked 


la) 


Ambulatory—with asslstonce 

Bathing In bathroom, or at bediide— Partlol 
Self 

Feeding—Complele Assistance (except IV. 
feeding) 

optional (does not af- 
fect clossification under 
these conditlons) 


Vislon Inodequote 
Oxygen therapy 


(b) Up ln chalr-—So!? 
Balhing al bodilda—Cormclala Assisionca 
Feeding Self, or requises food cu? of LY. 
feading 
Cnrygen therapy 
Vision Inadeauuulo 
lc) As In (b) wih the following changes 
Un In chair—With Assistance 
Bath at bedside 
Up In chalr— With Azalsionco 
Both at bedside-—-Partlal Self 
Feading—-Comnlata Assistance 
Vision inadequatz 
opllonal 
Oxygen therapy 
Both at bedalde 
Feoding—Self, or requires food cul, or LV. 
feeding 
Vision inadequate 
Oxygen tharopy 
(f) Being Speciaied of Nacessity {patient has 
continuous nursing assistonca to the extant 


optional 


ld 


lo 


Jorilanal 


thot maol rellef must be provided for speclal 
duty nurse) 


NOTE, Aay potlont whe othorwise folls into Cotegories | or 1, bur whe 
Is vador uction therapy of is ia isolation, Incontinen? finclhudiag wound 
droinayo necetsitoting chongo ol bed Haynl, of morhed’y emollenally 
disturbed Ineeds olmort constant obsorvotion, In single troom, ereales 
disturbonce! will be dropped to the neat category. 


Category tlt. Intensive, or "total", care 


All combinations not previously mentioned. 
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where 

N. = number of patients classified as Category I. 

N, = number of patients classified as 

Category II, 
N, = number of patients classified as 
Category III, 
and the coefficients of the equation represent the care 
hours required for each category. The total work load 
can then be found from 
WL = I + 20 hours (MINI. 2 ) 

where the additional factor of 20 hours represents the 
average time required for all other indirect care activi- 
ties. In a subsequent extension, a method for "controlled 
variable staffing" based on predictions from the above 
indexes was demonstrated by noting that the daily aver- 
ages and shortfalls in nursing time demanded, relative to 
the unit's basic staff, could be controlled through use of 
additional float personnel. Young [83] and Wolfe and 
Young [81] provide further insights on the practical 
aspects of the Connor staffing approach in their discus- 
Sion of day-to-day implementation procedures. Although 
essentially sound and more effective than traditional 
techniques used for staffing, the method has drawbacks in 
that it neither indicates the nursing skill levels re- 


quired nor suggests how available time might be allocated 





20 


between patients of various categories. 
Young [84] 

It might be anticipated that a patient classifi- 
cation system, such as that proposed by Connor, might 
not be universally applicable and should therefore be 
applied with some caution. Young, in a study of the 
direct care provided to pediatric patients found that, 
although pediatric patients can be classified into age 
groupings, no classification as to degree of self- 
sufficiency within such groupings was possible. Indeed, 
as opposed to adult patients, the more self-sufficient a 
pediatric patient became, the greater the demand on nurs- 
ing time. Examination of data on direct care nursing time 
revealed a continuous distribution of times with a mean 
of about 90 minutes; i.e., the frequency distribution of 
direct care times was unimodal indicating no distinct 
subpopulations. In order to determine total minutes of 
direct care required it was felt that a multiple linear 
regression model was most appropriate, using as the re- 
gressors the 0-1 variables corresponding to the subcate- 
gorization of 12 factors determined by professional 
nurses. These are shown in Table II.l. 

The numbers in the right-hand column are not to 
be interpreted as time required for the factors indicated; 


they are, in fact, regression coefficients. However, for 
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Table 11.1 
Young's Classification of Pediatric Patients 


Regression Coefficient 
for Determining 


FACTOR DEGREE Direct Care Times 
MOBILITY: Complete Self il 
With Assistance 10 
Bedfast 25 
ACTIVITY: Restrained 20 
Hyperactive 10 
Emotionally Disturbed 10 
Normal Activity 10 
Lethargic 1 
Unconscious 5 
FEEDING: Bottles x 6 5 
x 8 10 
x 12 15 
Feed Complete Assistance 25 
Feed Partial Assistance 9 
Feed Gavage 5 
Self 1 
TE AE Normal (Simple) 1 
Special (Complicated) 2 
MED. FREQ.: Average 1 
Frequently 20 
TREATMENT: Oxygen 20 
Dressing-Compresses 50 
Suction 20 
Tracheotomy 20 
Drainage 15 
Croupette 15 
Pediatric Procedures 20 
INCONTINENT: Diapers Frequently (10,20,30 per day) 1 
Bed Change 10 
TPR: Normal Frequency al 
Special 10 
ISOLATION: Protective 3 
Preventive 3 
FAMILY Assist -10 
MEMBERS: Non-Assist il 
Special Attention Given 2 
SPECIAL OF Special Nurse 2 
NEC: Relief Only IL 
BATHING: Partial Assistance 10 
Complete Assistance 15 


Add constant 


of 40 to total 
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those factors checked for a particular patient, the sum 
of the numbers, when added to a constant of 40 does pro- 
Teele=sthestotal direct care time, in minutes, required for 
that patient. Young found a very high correlation be- 
tween these predicted times and observed times. More 
Significantly, it was again found that patient assessment 
and, in this instance, "continuous" classification was 
much more predictive of patient care required than a 
census count and rule-of-thumb. 
Other Studies 

For the sake of completeness we briefly note 
some additional studies involving hospital inpatient 
classification. The National League of Nursing Education 
[48], Cloussen [15], Stanford [70], and the Hospital 
Association of New York State (HANYS) [30] all have pro- 
posed classification systems based on assessment of a 
patient's physical and emotional needs as well as the 
incidence and/or frequency of various medication or 
treatment regimens. Work by Poland, English, Thronton, 
and Owens [54] also involves the use of patient classi- 
fication; however, since the primary thrust of the study 
is staffing, we defer its consideration to the next 


section. 
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1 > Nurse Staffing in an Acute Care Setting 

Beginning with the concepts proposed by Connor, 
the study of the nurse staffing problem has intensified 
with an ever-increasing degree of sophistication in the 
models proposed. Typically, solutions were sought for 
the following three aspects of staffing as they relate to 


demand for service: 


1. Number and skill levels of nursing personnel 
needed; 
2. The distribution of the time of available 


personnel among direct and indirect care 
activities (the "allocation" problem) ; 

3. The delegation, by skill level, of the con- 

stituent tasks of direct and indirect care 
(the "assignment" problem). 
Each of the methods presented below represents an attempt 
to answer one or more of these questions. 
Wolfe [80], Wolfe and Young [81,82] 

In order to overcome some of the shortcomings of 
Connor's staffing models, Wolfe proposed a methodology 
based on the concept of a "multiple assignment model." 
This approach permits the assignment of a skill level to 
more than one activity as opposed to the classical 
assignment model which imposes one-to-one assignments. 


Wolfe viewed all nursing activities in both the direct 





24 


and indirect care categories as falling into one of 12 
mutually exclusive task complexes as shown in Table 

II.2. A work sampling study was undertaken to determine 
the nursing time required for each complex during an 
eight-hour day. These times were related to patient con- 
dition by regressing them on the number of patients in 
each of Connor's three classifications on the unit. The 
resulting regression coefficients, oe were interpreted 
as time required for the task complex per patient of 
category Nu, i=1,2, or 3; certain coefficients were found 
to be insignificantly different from zero and thus con- 
stant for all or a subset of the classifications. 


A cost, c indicating both salary cost and im- 


ij’ 
puted monetary disutility incurred by assigning skill 
level i to task complex j, was then determined. The lat- 
ter cost component was obtained by application of psycho- 
metric measurement techniques to judgments solicited from 
both nurses and patients. Wolfe's integer programming 


multiple assignment model is, therefore, 


MIN2% (c..-4b.)x.. + Za ay. 
1 AA iMi i 


2 
E (II.3) 
TD: = 0%: l 
ij 


y > 0 and integer 
al 





25 


Table II.2 


Task Complexes Identified by Wolfe 


Technical Tasks I ag 
Technical Tasks II 8. 


Evaluation of Pt. 9. 


Needs and Assign- 
ment 


Supervising and vO). 


Teaching 


Preparatory Care I ll. 


Preparatory Care Zi 


II 


Clerical Tasks I 
Clerical Tasks II 


Clerical Tasks III 


Housekeeping Duties 


Escorting and Emergency 
Errands 


Maintenance, Checking, and 
Ordering of Supplies and 
Equipment 
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where 
A; 1s the hourly wage of skill level i, 
a; 15 the number of hours worked by cach member 
of skill level i (Wolfe assumed aes oo"; 

2 is the assignment variable, a value of l 
indicating the assignment of skill level i to 
task complex j, and 

Y; is a variable indicating the number of person- 

nel of skill level i required. 
Closer examination of (II.3) reveals that the objective 
function calls for the minimization of an assignment cost 
plus a non-productivity cost, assuming the following: 

l. All task complexes must be assigned; 

2. Ample personnel of all skill levels are avail- 
able to meet the requirement Yai 

3. No partial shifts are allowed; 

4. No splitting of task complexes between skill 
levels is allowed. 

Probably because of these assumptions, Wolfe claims not to 
Mave touched upon the allocation and assignment problems 
noted earlier, but rather to have provided a means for 
determining the number and skill levels of nurses required 


on a unit. This is the primary drawback of his methodology. 
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Singer [68] 

Recognizing the stochastic nature of require- 
ments for nursing care as a function of a patient's con- 
dition, Singer attempted to improve the staffing aspects 
of Connor's work by modeling the fluctuations in the num- 
bers of patients in each of the three classifications in 
a time-dependent framework. He was able to do this by 
using a multiple-population, immigration-death process 
model, and provided good predictions of Connor's direct 
care index over relatively short horizons. The model 
did not perform as well for time periods greater than 
about three days. Unfortunately, few indications as to 
how these findings may be put into practice are given. 
Jelinek [32,33] 

In an effort to predict the implications of any 
number of several well known factors on some convenient 
measures of nursing activity, Jelinek proposed the use of 
a multiple linear regression model. A major result of 
such a methodology is isolation of effects through con- 
trolled experimental design. The author reports an appli- 
cation of this technique in regressing variables indicat- 
ing the number of hours per patient-day devoted to direct 
care, indirect care, and non-productive activities on a 
set of independent variables related to patient length of 


stay, health status (as per Connor's classification), type 
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Srnec, Untt census Mand staff size and skill level mix. 
Among the findings reported were decreasing returns to 
scale for direct care activities as a function of avail- 
able staff, with a corresponding increase in non- 
productive activity. The findings also emphasized the 
poor predictive power of patient census relative to nurs- 
ing activity. 

Poland, English, Thornton, and Owens [54] 

The authors offer a normative model for predic- 
tion of demand for direct care by means of an attractive- 
ly designed system called PETO. Through work sampling, 
the relationships between various criterion levels of 
seven patient assessment factors and the number of hours 
of required "physical nursing care" were determined. 
These times were then translated into a point scale for 
each criterion level. Thus, a patient would be assessed 


in accordance with Table II.3, then the point total con- 


verted to "physical care units" (PCU) by means of Table 
11.4. Finally, the unit direct care work load can be 
found by adding all PCU's. Based on nursing care audits, 


the amount of required care delivered on sample units was 
found to be significantly increased after installation of 
the system. Unfortunately, however, PETO failed to indi- 
cate how staffing by skill level as well as allocation of 


resources may be achieved. 





Element 


Diet 


toileting 


Vital signs 
and 
measurements 


Table 11.3 
PETO Elements of Care 
Description 


Feeds self without super- 
vision, or family or 
parent feeds patient. 

Feeds self with supervision 
of staff. 

Feeds self with constant 
presence of staff, or gas- 
trostomy feeding every 4 
hours. 

Total feeding by personnel 
or instructing the parent, 
Or Continuous I. V.; OY 
blood transfusion. 

Tube feedings more frequent- 
ly than every 4 hours. 


Toilets without supervision. 

Toilets with supervision, or 
specimen collection or 
uses bedpan. 

Up to toilet with standby 
supervision, or output 
measurement every hours, 
or daily colostomy irriga- 
tion. 

Incontinent, average output. 

Incontinent with diarrhea, 
or immediate postoperative 
colostomy, or urethrostomy 
or drainage with frequent 
dressing change. 


Routine--daily temperature, 
pulse, and respiration. 
Vital signs every 4 hours, 
or night check every 4 

hours. 

Vital signs monitored, or 
hypothermia, or vital 
Signs every 2 hours. 

Vital signs and observation 
every hour, or vital signs 
monitored, plus 
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Point Value 


12 


da 


12 








Element 


Respiratory 
aids 


Suction 


Cleanliness 


Tabac Liss, cont ca. 
Description 


hypothermia and 
neurological evaluation. 
Blood pressure, pulse, respi- 
ration, and neurological 
evaluation every 30 minutes. 


Bedside humidifier, or blow 
bottle. 

Mist or humidified bassinet 
when sleeping, or cough and 
deep breathe every 2 hours, 
or IPPB without supervision 
every 4 hours. 

Continuous oxygen, or cough 
and deep breathe every hour, 
or continuous assisted venti- 
lation. 

Mechanical respiratory aid, or 
IPPB with supervision every 
4 hours. 

IPPB continuously with inter- 
mittent hand ventilation. 


Routine postoperative standby. 

Nasopharyngeal or oral suction 
as needed 

Tracheostomy suction every 
hour, or nasogastric tube 
irrigation every 2 hours. 

Tracheostomy suction every 30 
minutes, patient responsive. 

Tracheostomy suction every 30 
minutes, patient not re- 
sponsive. 


Bathes self, bed straightened. 

Bathes self with help or 
supervision, daily change of 
bed. 

Bathed and dressed by person- 
nel, or partial bath given, 
daily change of linen. 

Bathed and dressed by person- 
nel, special skin care, 
occupied bed. 


30 


Point Value 


12 
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TapDle II.3 +cont'da. 


Element Description Point Value 
Turning and/ Up in chair with assistance 
or assisted once in 8 hours. 1 
activity Up in chair with assistance 

twice in 8 hours, or walk- 

ing with assistance 2 


Bedfast with assistance in 
turning every 2 hours, or 
up walking with assistance 
of two persons twice in 8 


hours. 4 
Bedfast with assistance in 
turning every hour. 8 
Turning on orthopedic frame 
every hour. 12 
Table II.4 


PETO Conversion of Points into Hours of Care 


Physical 
Care Units Hours of 

Points (PCUs) Nursing Care 

4-11 1 1 
12-19 2 2 
20-27 3 3 
28-35 4 4 
36-43 5 5 
44-51 6 6 

52 59 7 7 
60-67 8 8 
68-75 9 9 
70783 10 10 





32 


Shuman [66] 

Although the work of Shuman is primarily con- 
cerned with regional health system planning models de- 
signed to demonstrate methods of increasing productivity 
and minimizing the cost of providing services, it is 
worthy of further examination for two reasons. First, 
it provides insight into the model-building aspects of 
operations research as applied to health system problems. 
Secondly, an extension of Wolfe's model is presented 
which provides an answer to the allocation question not 
given in the original study. 

The extension is made in the following way: with 
reference to (II.3) the variables Xi are taken to be 
the number of times personnel type i performs task com- 
plex Jj. The Ze are integers. The Yas on the other hand, 
are now assumed to be continuous variables. Finally, an 


additional constraint set 


Zx.. =d, , yj 


is added to the model, where aj represents the number of 
times task complex j is demanded. Obviously, the task 
complex demands may be assigned to several personnel 
types, and the resulting model solution yields an indica- 
tion of personnel allocation. Shuman claims that by 
allowing Y; to be continuous, a part-time staffing policy 


can be inferred from the model solution. This, in fact, 
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may be difficult to do because the optimal values of Y; 
are as unlikely to terminate in half-time units as they 
are in full-time units. Thus, while the extended model 
is certainly an improvement over Wolfe's original model, 
it nevertheless is not completely representative of 

reality. Shuman provides a further modification to 

(II.3) indicating how this model may be used for several 


units simultaneously. 


II.4 Patient Classification in a Long-Term Care Setting 

Classification systems for long-term care have 
evolved essentially from those developed earlier for 
acute care. They have required some modification, mainly 
because of the need to specifically address the problems 
of aging and/or chronic disease. As will be seen, the 
studies have, in general, emphasized the following three 
areas: 

l. Degree of self-sufficiency in daily activities; 

2. Mental and emotional state; 

3. Rehabilitative and instructional needs. 
Katz, et al. [36] 

mho study of —Katz, etal, with respect to the 
recovery process of aged patients, is notable for its 
prescient recognition of the importance of degrees of 


patient independence in activities of daily living (ADL) 
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as they relate to nursing service required. The ADL's 
identified as being most important are 

l. Bathing, 

2. Dressing, 

3. Toileting, 

å. Transferring, 

5. Continence, 

6. Feeding, 
with each further subdivided into a tri-level evaluation 
of independence. An eight category "Index of Indepen- 
dence in Activities of Daily Living" is offered, as shown 
in abbreviated form in Table II.5. This instrument was 
tested on patients representing numerous primary diagnos- 
es. Katz found that almost universal applicability was 
possible regardless of diagnosis. Further, the phenome- 
non of increasing demand for nursing services as an 
ordered function of the index category was observed. The 
existence of factors of primary significance among the 
six in predicting degree of patient dependency was con- 
jectured, but little follow-up study is reported. 
McKnight [41] 

As part of a study conducted to determine the 
amount and type of care received by nursing home patients 
in the Denver area, McKnight developed a set of "Criteria 


for Patient Nursing Need" which was used to specify a 
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Table II.5 


katz, comal Classifications ot Independenceen ADL 


Category Description 
A Independent in feeding, continence, trans- 
ferring, going to toilet, dressing, and 
bathing. 
B Independent in all but one of these functions. 
C Independent in all but bathing, and one 


additional function. 


D Independent in all but bathing, dressing, and 
one additional function. 


E Independent in all but bathing, dressing, go- 
ing to toilet, and one additional function. 


F Independent in all but bathing, dressing, go- 
ing to toilet, transferring, and one addi- 
tional function. 


G Dependent in all six functions. 


Other Dependent in at least two functions, but not 
classifiable as C, D, E, or F. 


three level (minimum care, moderate care, maximum care) 
patient classification system. These levels of classifi- 
cation are essentially analogous to the three levels of 
care Intermediate B, Intermediate A, and Skilled nursing 
care defined in Chapter I. The criteria identified were 
as follows: 

l. Diagnostic description (primary and secondary); 


2. Nursing activities (medications, treatments, 





36 


Special diets, etc.); 

3. Physiological factors (vision, hearing, con- 

tinence, mobility, etc.); 

4. Psycho-social factors (emotional and social 

behavior); 

5. Rehabilitative factors (functional ability, 

both mental and physical). 
Classification was dependent upon the similarity of 
patient condition with several descriptive items speci- 
fied for each criterion and category. 

A major portion of this study was devoted to de- 
termining the amount of care received by patients of the 
various classification levels. Time study analysis was 
undertaken to determine the performance times for 25 
nursing care areas, both by patient classification and 
nursing skill level. In addition, activity sampling was 
performed in several of the nursing homes under study to 
give an indication of how nursing personnel time was allo- 
cated. It should be noted that the Denver work included 
the replication and extension of a study conducted by the 
City of Milwaukee Health Department [47]. 

Both the patient classification and nurse staffing 
aspects of McKnight's efforts will play a large role in 
our work. More will be said concerning this study in 


Chapter III. 
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Salmon, et al. [61,62] 
INP IDS system developed by Salmon, et al. has 
as its major objectives 
l. The forecasting of patient needs on input to 
the long-term care system, and 
2. Placement assessment, both initial and follow- 
Up: 
A rating scale of one to five, one indicating no effort 
required and five indicating maximum effort required, was 
attached to each of the six areas denoted by the acronym 
- Restorative procedures 
- Activities of daily living 
Problem behavior 
- Illness 


- Dependency, general 
- Social service 


NOH DPD 
i 


In applying this system, the authors discuss its successes 
in meeting the second of their objectives, but similar 
outcomes were not obtained for the first goal. No attempt 
was reported concerning the relationship between patient 
classification and staffing levels. 
Burack [13], Burack and Denson [50] 

Noting the need for an interdisciplinary set of 
indicators of geriatric patient status, Burack proposed 
a five-item list of guidelines as follows: 

l. Medical classification (related to serious- 


ness of disease): 
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2. Functioning performance (independence in 
ADL) ; 

3. Projected goals (restoration of activity 
levels); 

4. Therapeutic category (need for professional 


service, medications, or treatments); 
5. Services required (medical, nursing, psy- 

chiatric, rehabilitative, social work). 
Notice that the concept of goal-oriented care is empha- 
sized in this system, a feature that has typically been 
relegated to minor status in other classification systems. 
Burack encourages further subdivision of the basic guide- 
lines by care providers to better reflect their own needs. 
It is unfortunate that no indication of how this classi- 
fication may be translated into nursing action is given. 
Parker [52] 

In the first of two studies designed to elicit a 
mathematically-based geriatric patient classification sys- 
tem, Parker examined the relationship between the 19 
health status indicators and 4 health-care need status 
categories shown in Table II.6. Given that the binary 
peto S represents a “patient profile in terms of the 
indicators and dis i=1,2,3,4, the health care status, it 
was demonstrated that Bayes' Theorem could be used to 


determine the probability of a patient with profile S 





Indicators (=1 if present, =0 if 


l 


10. 


Table II.6 


Used by Parker 


Single (=1), Married 
(=0) 


Drugs prescribed (=1 
if at least one 
drug presceibed) 

Cancer 

Senility 

Vascular lesions of 
the central ner- 


vous system 


Diseases of heart 


Diseases of skin 
Argehr reis 


Decubiti, ulcers, 


fistulas 


Total or occasional 
incontinence 


je: 


12% 


t3. 


14. 


ED 


16. 


19% 


18. 


19: 
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Definitions of the Indicators and Categories 


absent) 

Any social disability 

Immobility--less than 
independent ambula- 
tion on level surface 


Needing to be fed 


Special treatments 
needed 


Major assistance dressing 
and bathing 

Conscious but unable to 
communicate needs 


Aberrant behavior 


Mild confusion 


Severe confusion 


Health-Care Need Status Categories 


1. 


2 


3. 


4. 


Private home 


Personal care, home 


Skilled nursing home 


Chronic disease hospital 
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being in category a; rom 


4 
P(a.|s) = P(S fay) Play) E p(s |a,)P(d,) (11.4) 
where 
P(sja,) = p(s, |a,) i p(s,|S,Ma;)-...*P(S¡9 ETT. s) 


and S; indicates the presence of the pen indicator. Using 
such a rule to classify 1,245 Medicare/Medicaid patients, 
Parker achieved a prediction rate of 70%. Methods were 
suggested for use in extracting a subset of the original 
variables which appeared most important in computing the 
categorization. These "key indicators" were found to be 
(in order of importance) 

Pee oor rty ; 

2. Continence, 

3. Major assistance dressing and bathing, 

4. Any special disability, 

5. Severe confusion, 

6. Special treatments needed. 
Reducing the set of indicators by using only the key 
indicators above, a prediction rate of 70% was also 
achieved. Further comments on this approach will be made 
in Chapter IV. 
Parker and Boyd [53] 

In the second of his two studies, Parker explored 


both the use of and relationships between a discriminant 





41 


analysis approach and a multi-hierarchical cluster analy- 
sis approach to patient classification. Taking as data 

a sample of patient profiles derived from Collaborative 
Patient Assessment Instrument (CPAI)* scores on some 60 
items, stepwise discriminant analysis was first applied 
in order to relate these items to one of six appropriate 
levels of care. An overall prediction rate of 77% was 
achieved. The clustering technique was then applied in 
order that patients with similar profiles might be 
grouped together. Five clusters were ultimately dis- 
covered, roughly analogous to the original six-level 
categorization. Taking cluster centroid vectors as being 
indicative of membership in that cluster, Parker proposed 
that a candidate's scores be added and compared to the 
totals for the five clusters as a rough classification 
rule. 

Key variables identified by each method are shown 
in Table II.7. Those from stepwise discriminant analysis 
were determined by order of entry into the discriminant 
function, while cluster analysis variables were inferred 
from examination of how well each variable explained 


cluster membership (the "B/T" criterion). 


* A more complete discussion of this instrument and its 
development is given in Chapter III. 
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Further elaboration on both the CPAI and the 
statistical techniques employed will be provided in sub- 
sequent chapters, aS well as analysis of the results. 
Note, however, the importance socio-demographic informa- 
tion plays in both findings, as well as the fact that no 
easily applied classification scheme was tested for its 
predictive capabilities. Neither of the results were re- 
lated to staffing procedures. 

Other 

Two other classification studies are noted brief- 
ly. Ashton [3] developed a set of nursing criteria for 
evaluating extended care patients, the criteria being re- 
lated ultimately to a tri-level (self-care, intermediate 
care, intensive care) classification system. Major fac- 
tors identified were continence and mobility, mental and 
social status, independence in the ADL, and clinical nurs- 
ing procedures required. 

Shaughnessy, et al. [65] proposed a classifica- 
tion based on seven major factors (mental status, communi- 
cation, sensory deprivation, physiological disturbances, 
physical dependence, environmental and psycho-social 
items) in order to assist in developing a triage catego- 
rization of patients. The categorization was then relat- 


ed to general personnel allocation guidelines. 





Table II.7 
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Key Variables Determined by Discriminant 


and Cluster Analysis by Parker and Boyd 
Discriminant Analysis Cluster Analysis 
l. Patient's present location 1. Communication of needs 
2. Dressing 2. Present location 


3. Length of time at location 3. 


4. Behavior pattern 4. 
5. Age ae 
6 Stair climbing 6% 
7. Total family income TR 
8. Dentition 8. 
we) Birthplace on 
10. Toileting 10. 
11. Mobility level Ir 


(Abbreviated list) 


II.5 Nurse Staffing in a Long-Term 


With one notable exception, 


Research site 

Appropriate level of care 
Dressing 

Wheeling 

Toileting 

Mobility 

Bathing 

Feeding 


Bowel function 


Care Setting 


staffing studies re- 


lated to long-term care facilities have relied largely on 


the techniques of computer simulation to explore the 


questions of personnel mix, allocation, and assignment. 


We shall discuss three such simulation studies, and a 


mathematical programming approach. 
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For purposes of the discussions which follow, we 
note that computer simulation usually involves the model- 
ing of a complex system by means of interrelated mathe- 
matical functions. Such an abstraction attempts to cap- 
ture the essence of the manner in which various factors 
affect the system both directly and in combination with 
other variables. The model may be cast into either a 
descriptive or normative framework. Wagner [77] points 
out two basic caveats for such an approach: 

l. If uncertainties are inherent in the model 

(e.g., the simulation of demand for service 
based on some probability distribution) re- 
sults are subject to statistical error; 

2. In order to make the system model computa- 

tionally tractable, much abstraction and many 
Simplifying assumptions may be required, with 
the undesirable result of divorcing the model 
from reality. 

It seems clear that the authors cited below were aware of 

such difficulties. 

Turner, et al. [74] 

Noting the need in a nursing home environment for 
the determination of approximate resource levels (nursing, 
equipment, non-medical support, and outside medical 


assitance), and a reasonable scheduling system for such 
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FESDUIACESWESWcouUltWDe allocatéd, Tufner, “et al.*proposed 
two simulation models with these goals in mind. These 
models incorporated a list of 28 direct care activities 
and 22 supporting (or indirect care) activities, and 
showed that about 80% of such activities could be sched- 
uled, the remainder occurring more or less randomly. As 
stated by the authors, the models do not allow for para- 
metric examination of the task performance time estimates 
used. Proposed measures of effectiveness by which vari- 
ous configurations were compared were the following: 

l. Patient waiting time for service; 

2. Specific task/skill level assignment patterns; 

3. Amount of staff idle time; 

4. The attainment of minimum service requirements. 
McKnight and Steorts [43] 

The model proposed by McKnight and Steorts repre- 
sents a significant movement of simulation studies toward 
reality. Quite explicit descriptions of the nursing home 
environment in terms of number and size of units, patient 
mix (as per McKnight's classification system cited earl- 
ier), services offered, scheduling and occurence intervals 
of demands for service, staffing constraints, and task/ 
skill level preferences are considered. 

Beginning with the individual patient care plans 


translated into specific tasks and then aggregated within 
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classification levels, the model generates demands, both 
scheduled and random, with required service times taken 
from a normal curve fitted to task time-study data. To 
each demand is attached a priority, being a function of 
both patient classification and the approach of a cutoff 
time beyond which the demand would be unmet. Specific skill 
level assignments are made on the basis of a first and 
second choice preference structure. 
In addition to generating financial reporting in- 
formation, the model also reports the following: 
1) Percentage of patient demands met, both with 
and without waiting; 
2) Percentage of first and second choice person- 
nel assignments; 
3) Hours of care assigned by task/skill level; 
and 
4) Staff utilization information. 
Thus, from the manager's point of view, the relative 
effectiveness of various staffing configurations might be 
tested by making several model runs and comparing them 
on the basis of the above criteria. 
Hundert [31] 
In a very recent study, Hundert describes a two- 
phase nursing home simulation model whose dual objectives 


are the provision of the maximum amount of care at 
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minimum cost. The system takes as its input the follow- 
ing parameters: 

1) Schedule of services--exact times or time 
intervals for task performance; 

2) Server preferences--a first and second choice 
by task and patient classification; 

3) Task priorities--set initially by the user, 
then updated in much the same way as in the 
McKnight and Steorts model; 

4) Patient demand profile--a translation of the 
patient care plans in terms of the probabili- 
ty that a patient would require one of the 
specified services, with frequency and sched- 
uling information also provided. 

Through the use of time-study data, Hundert modeled ser- 
vice time distributions as gamma functions, which the 
simulation model can use to infer required task perform- 
ance times. 

Based on the output of this portion of the simu- 
lation, the cost phase of the model predicts average cost 
per patient day by means of linear regression estimates 
derived from historical data. Not surprisingly, it was 
found that nursing salaries, facility capacity and age, 
and patient census were the most important factors in the 


cost model. 
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The entire model was validated and used to deter- 
mine the optimal unit size, staffing mix and levels, 
and service schedules for a proposed 180-bed nursing home. 
Acceptable results are reported for the first portion of 
the model, but cost estimates were deemed to be in ques- 
Cion. 

Liebman [38] and Liebman, et al. [39] 

The work of Liebman represents a significant de- 
parture from earlier nurse staffing studies in the long- 
term care field in its reliance on the use of integer 
programming allocation and assignment models. In an 
attempt to develop "team effectiveness profiles" for com- 
parison of various nursing team configurations in a 
hospital-based extended care unit (ECU), a short-term 
(or daily) assignment model and a long-term allocation 
model were proposed. Taking i as the index of individual 
patients, j the index of skill levels, and k the index of 
nursing tasks, the short-term model is formulated as 


1 
MAX = A Ciktikřijk 


cx 
Er 2 Zik =m Vi,k 
J ` (II.6) 
BE tikřijk T ’j’ J 
= 0, 1 
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and the long-term model is 


1 
MAX ork 22 SKY ik 
jk 
S E Sy,, = 1, Vk 
or Fe) 
y 
£ T Yik 3 j 
k 
> 
er 0 
The variables on in (II.6) indicate the assignment 


pattern of nursing skill levels to specific tasks and 
patients, where en implies the assignment is made. In 
(II.7), the Y 4% denote the fraction of demand for task k 
allocated to skill level j. The ean are skill level-task 
preference measures solicited from professional nurses by 
means of a Q-sort psychometric measurement technique. 
Although interval scaled measures were thereby obtained, 
equally acceptable results uSing corresponding ordinal re- 
lationships were reported. The tik and Ti are times re- 
quired to perform specific tasks on an individual patient 
level or aggregated over all patients, respectively. 

These were determined by asking nurses how long each pro- 
cedure usually took. Nursing time available to devote to 
direct care activities, Sj, was determined by aggregating 


task performance times for each skill level, again accord- 


ing to current assignment procedures. The constants c* 
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and c** are normalizing factors. Both models require 
that all tasks demanded must be accomplished. A com- 
binatorial method relying on complete exhaustion of 
available time supplies was used to solve (II.6), while 
(II.7) was easily formulated as a network flow problem. 

The validity of the model was tested by comparing 
the most effective assignment plans and team configura- 
tions computed with existing practices on the unit. Pro- 
fessional nurses were asked to state their preference 
between pairs of such results, with the model solutions 
usually being favored. 

While this study represents a significant effort, 
it does leave room for improvement. First, its concepts 
are definitely hospital-oriented, albeit for an ECU, and 
the findings may therefore pose some difficulties for im- 
plementation in the typical long-term care facility. Task 
performance times and the method of computing the E tend 
to cause the solution to reflect heavily the status quo. 
Additionally, it is implicitly assumed that existing 
patient care plans may be readily formulated for use in 
the model in terms of the nursing task list presented. 
Finally, no provision is made for determining how con- 
strained nursing resources can affect the amount of direct 
care given, and how choices related to task priorities 


for various patient types might be made and evaluated. 
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Mehta [44] 

We briefly note a recent study by Mehta, which 
essentially quantified existing staffing and assignment 
procedures in sample facilities and categorized activi- 
ties as direct care, indirect care, travel, communication, 
and personal time. By work sampling, the time devoted to 
each task was obtained and an indication of the skill 
levels performing each task was given. The staff mix was 
determined by simply dividing the aggregate time spent by 
each skill level by seven hours to obtain full-time 
equivalent personnel. The use of a float staff was 
recommended to handle time for which no accounting could 
be made by the above procedure. Not surprisingly, the 
staffing recommendations that were provided closely re- 
semble the staff mix that existed on the units at the 


time of the study. 


II.6 Work of the "Four University" Group 

We conclude our survey by noting the contribu- 
tions made in several areas of the long-term care field by 
the so-called "four university" group. Established in the 
late 1960s, the study team consisted of those researchers 
listed in Table II.8 with primary research interests as 
indicated. Recognizing the many schemes for long-term 


care patient classification already in existence, the 
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Table IT-8 


The "Four University" Research Group 


Group Location Principal Investi- Purpose of Primary 
w gator Interest 

Case-Western Dr. S. Katz* Studies of course of 
Reserve Dept. Dr. M. Stroud** illness in chronic 
Preven. Med. disease 

Harvard Center Dr. Paul Donsen Evaluation of programs 
for Community Mrs. Ellen Jones for care of the 
Health and chronically ill 
Medical Care 

Hospital Assn. Mr. L. Danehy Community resource 

of New York allocation 

State 

Johns Hopkins Dr. C. Flagle Facilitating decision 
Sch. of making in management 
Hygiene & Pub. of patients 


Hlth. Opera- 
tions Research 
Division 


e Now at Michigan State University 
xx Now at University of Pennsylvania 


group saw its main task as the development of an assess- 
ment system that would not only capture the important 
aspects of then current systems, but also establish a 
"common language" of patient classification that wald be 
interdisciplinary in scope. The result was the Collabor- 
ative Patient Assessment Instrument (CPAI) alluded to 
earlier, along with its accompanying user's manual [35]. 


A detailed discussion of the CPAI will be undertaken in 
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the next chapter. For the present, it is sufficient to 
call attention to the statement of Jones [34], which 
points out the significance of the CPAI as a patient- 
oriented rather than process-or service-oriented system. 
Although the original group has officially dis- 
banded, the individual research teams have continued to 
pursue their interests in the LTC field in a collabora- 
tive manner. Among those studies currently underway are 
the improvement of the assessment tool as a means for im- 
proving the quality of patient assessment, and additional 
work on the relationships between patient classification, 


demands for care, and facility nurse staffing. 


II.7 Summary 

This chapter has attempted to describe briefly 
some relevant studies concerning patient classification 
and nurse staffing in both the acute and long-term care 
environments. Commonalities in conceptual approaches, 
significant improvements, and areas for further explora- 
tion have also been highlighted. It is against this 
background that the models and procedures proposed in this 
dissertation are presented. It is hoped that the his- 
torical perspective that has been provided will prove 


helpful in evaluating the new methodologies. 





Chapter III. Development of the Methodology 


ee Introduction 

Having discussed the history and current state of 
the art of nurse staffing methodologies we turn now to 
the development of the proposed models. Recall that we 
essentially are seeking a process that will assist a LTC 
facility administrator in solving personnel budgeting 
problems by presenting alternative courses of action for 
consideration. Unlike some previously designed systems 
of this sort, however, the approach suggested here relies 
heavily on the concept of patient needs as related to de- 
mand for nurSing service. As a result, supplementary in- 
formation is concurrently produced that will benefit the 
director of nursing in allocating the predicted nursing 
resources. 

After an initial discussion of the general prob- 
lem to be solved, we will demonstrate how patient needs 
may be related to nursing activity through classification. 
The insights obtained will lead directly to the construc- 
tion of a mathematical programming model and two modified 
formulations derived from it. These models will serve as 
the basis for solving the budgeting and nurse staffing 


problems. 


54 
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111.2 The Administrator's Problem 
As was pointed out in Chapter I, the LTC facility 
administrator, as the primary managerial decision-maker, 
is responsible for the total allocation of available re- 
sources so that the essentials of institutional life may 
be provided. Given that the bulk of operating expenses 
are attributable to personnel costs, especially those of 
nursing staff, we chose to concentrate primarily on the 
manner in which such expenses could be related to the 
provision of service. Recall, too, that the five con- 
siderations identified by Aydelotte [5] (quality of care, 
patient characteristics, staffing predictions, staffing 
logistics, and quality measurement) were postulated as 
being of consequence in any budgeting process.* Addi- 
tionally, certain exogenous variables must be considered 
by the administrator in reaching an acceptable decision 
as to personnel budget and the staffing mix it would imply. 
These variables include the following: 
l. State and federal requirements for minimum 
staffing and skill level coverage (e.g., R.N. 
On duty at least one shift per day, adherence 
tota staff to patient ratio of 1:25, etc.) ; 
2. Reimbursement rates for Medicare/Medicaid 
patients; 


3. Demand intensity implied by the facility's 


*These items were discussed in Chapter I. 
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patient mix. 

The synthesis of all of the above-mentioned fac- 
tors into a normative, decision-making framework, along 
with such intangibles as the desire to maintain a good 
professional reputation, is part of what Baloff, Aber- 
nathy, and Hershey [8] identify as the "aggregate budget- 
ing process." The decision variables involved in this 
process are the budget to be allocated to nurse staffing 
and the staffing levels and skill mix to be maintained. 
Essentially, then, the administrator's problem may be said 
to be 

How to fix budget and staffing levels in a manner 

responsive to patient needs and other factors, 


both for the short run and as a result of future 
alterations in the exogenous and/or endogenous 


considerations mentioned earlier (e.g., legal re- 
quirements, reimbursement rates, patient mix, 
etc.). 


EGE, ecan [8, pea] succinctly capture the essence of 
why this problem is non-trivial: 


...The difficulties of nurse staffing planning are 
partially caused by two insurmountable problems 
for which no complete remedies exist: (1) the 
impossibility of developing perfect demand fore- 
casts for the future, and (2) the inability of 
specifying exactly and unequivocally the manner 

in which demand should be translated into required 
staff. 
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To these points we might also add the inherent pitfalls 
in the intermediate step of relating patient health status 
to demand. 

Recognizing the major areas of difficulty in nurse 
staffing planning, as noted above, it is perhaps wise to 
view the administrator's problem in terms of the resolu- 
tion of these problems. Certain of the exogenous changes 
to which a response is desired may be captured in the 
model which translates demand into required staff. The 
relationship between patient health status and demand, 
along with the development of future demand forecasts, 
may be viewed as a separate but interrelated entity. 

While we do not propose to predict future patient status 
in a manner analogous to Singer [68], we may rely on the 
relative stability of the patient mix within LTC facili- 
ties to lessen the need for such a procedure. Then, too, 
a method of assessing patient health status which may be 
readily applied at frequent intervals decreases the neces- 
sity for long-range demand forecasting, since it is 
usually the case that such predictions are required in the 
absence of better information. With this overall solution 
framework in mind, we turn now to the development of a 


demand-related patient classification system. 
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III.3 A Patient Classification System 

The lack of a truly useful patient classifica- 
tion methodology for our purposes in the long-term care set- 
ting has been previously noted. Such systems as do exist 
fall short of the goal of relating patient health status 
to demand for nursing service for essentially two reasons: 
(1) health status is, in fact, not adequately reflected, 
and (2) demand predictions are overly simplified and not 
readily applicable to the aggregate budgeting process. By 
way of example of the first shortcoming, we cite the use 
of patient "level of care" placement. Since the place- 
ment decision must, of necessity, include the considera- 
tion of several other factors in addition to patient 
health status, these factors may at times be so important 
as to alter any decision based solely on health status. 
Then, too, the time lags involved in the reassignment of 
a patient whose health status has changed can cause the 
current placement to be a poor predictor of nursing needs. 
The second shortcoming is exemplified by the PETO system, 
where the predicted demands for nursing services are 
stated in terms which do not allow for determination of 
the skill level mix as it relates to budgeting. It should 
be noted that other methods, such as those proposed by 
Burack [13,50] and Katz [36] were never shown to be re- 


lated to demand for nursing service. 
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Perhaps the ideal vehicle for reaching our stated 
goal is the intermediate step of a patient care plan. 
Having evolved from the traditional nursing care plan by 
including patient data of interdisciplinary interest, 
the patient care plan in its most recent form is usually 
presented in a "problem-oriented" format, and relies on 
a thorough consideration of patient status. Such a 
structure requires recording of both subjective and ob- 
jective problem identification and assessment, a plan of 
approach, and the specific actions required to ameliorate 
the difficulties. Often included are expected goals to 
be attained as a result of treatment. The implementation 
of this approach at the nursing level is shown sche- 


matically by Bower [10, p.10] in Figure III.1. 


Figure III.1 


Bower's Process of Planning Nursing Care 


Phase I 
Assessment 









Phase II Phase III Phase IV 
Problem Formulating Plans Planning for 
Identification for Nursing Action Evaluation 





mplementation 
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It is ot note that this form of patient care plan in- 
corporates a one-to-one correspondence between patient 
status and nursing demands. Aggregated over all patients 
in a facility, these demands would provide a reasonable 
indication of required staffing. There are, however, 
three drawbacks to this approach: (1) skill level prefer- 
ences for task assignments are often omitted, (2) exami- 
nation of alternative budget levels and skill level mixes 
would be difficult, and (3) the complex relationship 
between the myriad of possible nursing actions and per- 
formance times would be difficult to identify and stan- 
dardize. On amore mundane level, the problem-oriented 
patient care plan has yet to be fully implemented ina 
Significant number of long-term care facilities. Indirect 
use of care plans through translation of indicated activi- 
ties into tasks listed on a predetermined schedule has 
been suggested by Liebman [38,39] and in the several simu- 
lated models cited in Chapter II. Such procedures, how- 
ever, can be extremely time-consuming, and can require a 
great deal of professional nursing judgment to implement. 
Given the state of the art in patient classifi- 
cation systems and the inherent difficulties in the use 
of patient care plans for budgeting purposes, there is a 
need for a classification system emphasizing ease of 


application and comprehensive assessment and demand 
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prediction. That is, the categorization of patients into 
distinct, homogeneous groups should take place with the 
following criteria in mind: 

1. The classification decision should have its 
roots in a comprehensive assessment of a 
patient with respect to functioning status, 
impairments, medical status, psycho-social 
Status, and medically defined conditions; 

2. The classification groupings must bear some 
well-defined relationship to demand for nurs- 
ing services. 

In keeping with the first of these guidelines, the 
proposed classification system is to be based on the 
assessment information contained in the Collaborative 
Patient Assessment Instrument (CPAI) developed by the 
"four university" group.* A listing of the classifica- 
tion descriptors is given in Table III.l, and an illustra- 
tion of the instrument itself is shown in Appendix A. As 
can be seen from the descriptors included, the overall 
theme of the assessment is problem-oriented instead of 
diagnostically dependent. Jones [35] further identifies 


the following four characteristics of the descriptors: 


* The organization of this group was discussed in 
Chapter II. 





Table III.1 
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CPAI Classification Descriptors 


Identifying and sociodemographic items: 


Date of classification 

Interviewer or Classifier 

Program or project identification 

Patient's record or study number 

Social security number or health 
insurance benefit number 

Birth date 

Birthplace 

Sex 

Race 

Marital status 

Religious preference 

Residence address 

Patient's location at time of 
assessment 

Length of time at location 

Usual living arrangements 

Number of living children 

Education 

Usual occupation 

Employment status 

Total family income 

Health care coverage 


Functioning status items: 
Mobility level 


Transferring 

Walking 

Wheeling 

Stair climbing 

Bathing 

Dressing 
Eating/feeding 
Toileting 

Bowel function 

Bladder function 
Communication of needs 
Orientation as to time, place, person 
Behavior pattern 


Impairment items: 


Sight impairment 
Hearing impairment 
Speech impairment 
Fractures and disloca- 
tions 
Joint motion disorders 
Joint pain and swelling 
Missing limbs 
Paralysis/paresis 
Dentition 


Medical status: risk factor 


measurements: 

Height 

Weight 

Blood pressure 

Blood cholesterol 

Kidney function measure- 

ment: BUN 
Albuminuria 
Cigarette smoking 


Certain medically defined 

conditions:* 

Alcoholism 

Anemia 

Angina and/or myocardial 
infarction 

Arthritis 

Cardiac arrhythmias 

Congestive heart failure 

Diabetes mellitus 

Drug abuse 

Hypertension 

Malignancy 

Mental illness 

Neurologic disorders 

Respiratory disease, 
chronic 


* In addition, all medical diagnoses may be listed. 
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l. Patient orientation: the individual is 
describpef las he or shelactually is, not in 
terms of services being rendered or physical 
location; 

2. Multidimensionality: The scope of the assess- 
ment is a broad one; 

3. Objectivity: little subjective or interpre- 
tive analysis is required, hence results are 
reproducible; 

4. Relevancy: the items included have been shown 
to be strongly related to the major issues of 
morbidity and mortality. 

In addition, the subcategorizations provided with the 
functioning status, impairment, and medically defined 
items serve to indicate the patient's degree of dependency 
with respect to each descriptor. The incorporation of 

the suggested coding provided in the CPAI User's Manual 
[35] thus yields an ordinally scaled set of variables for 
purposes of analysis. 

As part of their developmental work on the CPAI, 
three of the research groups undertook the field testing 
of the instrument in various locations. As a result, a 
total of 623 assessments on Medicare/Medicaid patients 
was obtained, the sample consisting of those patients 


whose placement levels were as shown in Table III.2. As 
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Table ITI.2 


CPAI Sample Categorized by Level of Care and Research Group 


Research Group 






Level of Care HANYS Johns Hopkins Harvard | Total 
Private Residence PE AA 47 
Rented Room(s) — ise! 2 
Domiciliary/Personal i 

Care Facilit 18 
ECF/Supportive Nursing u 100 
core | el 
Chronic Disease/ we 

Rehabilitation Hospital 295 306 

Total ae a 295 623 


an additional and highly important part of this study, 
criteria for judging a patient's most appropriate level 
of care were established and utilized by various health 
professionals to obtain an indication of the divergence 
between legal placement and placement based on judgment 
of patient needs. The resulting categorization, which 
illustrates the observed divergence, is shown in Table 
ieee. 3. 

The appropriate placement decision was made, for 


the most part, on the basis of the patient's health and 
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Table III.3 


CPAI Sample Categorized by 
Appropriate Level of Care and Research Group 


Appropriate Level of Care HANYS Johns Hopkins Harvard | Total 


Private Residence/Room/ 1 59 4 64 
Home with or without 

Home Care 

ICF/Supportive Nursing 

Care Facilit 34 77 54 165 
Skilled Nursing Home 

(ECF/ECU) 91 45 205 
Chronic Disease/ 

Rehabilitation Hospital 157 163 


Other Hospital 


Total 202 | 29s | 62 


mental status. Where there were some overwhelmingly im- 





portant factors among the patient's sociodemographic 
characteristics (e.g. relatives living at home, insurance 
arrangements, etc.), the placement decision may have been 
altered. Similarly, borderline placement questions were 
usually resolved by recourse to this information. By and 
large, however, the adjudged appropriate level of care 
provides the high degree of association with patient 


status (obtained through assessment) that we seek in 
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keeping with the first criterion for classification. 

Quantification of the relationship between this 
judgmental variable and the multivariate indication of 
health status provided by the CPAI can therefore become 
the basis for an operable patient classification system. 
Due to the nature of the data involved, and in keeping 
with the need for a readily-applied methodology, a re- 
gression-based multivariate statistical approach has been 
utilized. Further explanation of specific methods, 
rationale for the approach, detailed description of the 
application of the method of choice, and the resulting 
classification system are provided in Chapter IV. 

There remains the consideration of the link 
between the derived classification system and nursing 
service demands. When others have attempted to quantify 
these demands it has typically been the case that the 
time necessary to perform the several tasks identified 
with the functions of nursing are presented. These per- 
formance times are usually, although not necessarily, 
obtained from time-study experiments, some typical ex- 
amples being given in Liebman [38,39], McKnight and 
Steorts [43], and Hundert [31]. Unfortunately, the task 
listings usually contain anywhere from 65 to 75 items 
which, when combined with the various patient classifi- 


cations and nursing skill levels ina staffing model, can 
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lead to models too large for most computers to solve. 

In a nursing home research study conducted in 
Denver, Colorado, however, McKnight [41] proposed the 
grouping of tasks typically performed in long-term care 
facilities into 25 "care areas" as shown in Table III.4. 
Notice that in addition to tasks which have traditionally 
been identified as "direct care" items (i.e., performed 
in the presence of the patient), there is also involved 
an area identified as "associated nursing care functions." 
The constituent tasks of this grouping, such as charting 
and rounds, are certainly more patient-oriented than 
other "indirect care" tasks performed by nurses. For this 
reason McKnight [42] identifies the 25 care areas as 
"patient centered" instead of as direct care groupings. 

A major portion of the Colorado work was devoted 
to time studies of nursing activities in 14 homes general- 
ly judged as giving "good" quality care. Specifically, 
the amount of time devoted by nursing personnel to the 
constituent tasks of the 25 care areas over a 3 1/2-day 
period for 195 sample patients was obtained. The findings 
forte haiole 711.5 indicate the average time, per 
occurrence, for each care area by patient classification 
(minimum care, moderate care, maximum care). Although 
this is a common way of presenting this type of result, 


Such data are difficult to apply without prior knowledge 
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Table III.4 
Nursing Care Areas Identified by McKnight 


ELIMINATION 
we Tolmer 
2. COMMODE 
3. BEDPAN 
4. URINAL 


CATHETERIZATION (AND RELATED PROCEDURES) 


1. CATHETERIZATION 

2. INSERTION OF FOLEY CATHETER 
3. BLADOER IRRIGATION 

4. EMPTYING DRAINAGE BOTTLE 


DAILY CARE 
l. MORNING CARE 
2. EVENING CARE 
3. REST PERIOD 
4. NIGHT CHECK 


BATHING 
1. COMPLETE BED BATH 
2. PARTIAL BED BATH 
3. PARTIAL BATH 
4. SHOWER 
5. COMPLETE TUB BATH 


BED MAKING 
1. OCCUPIED 
2. UNOCCUPIED 
3. STRAIGHTENING BEDS 


CARE OF SKIN 
t, PREVENTIVE SKIN CARE 
2. TREATMENT OF DECUBITUS 
3. GENTLE MASSAGE 


CARE OF HAIR AND FACE 
1. COMBING AND BRAIDING 
2. SHAMPOOING 
3. CUTTING 
4. SETTING HAIR 
5. SCALP TREATMENT 
6. SHAVING 


ORAL HYGIENE 
1. BRUSHING TEETH 
ce. DENTURE CARE 
3. MOUTH CARE 
4. CARE OF LIPS 


CARE OF NAILS 
1. FINGERNAILS 
2. TOENAILS 
3. MANICURE OR PEDICURE 


REHABILITATION ano RECREATIONAL 
1. PASSIVE EXERCISE 
2. SUPERVISED EXCRCISE 
3. DIVERSIONAL 


FOOD ano NOURISHMENT 
. BREAKFAST 
2. LUNCH 
3. DINNER 
4. EXTRA FEEDINGS 
S. WATER 





12, 


132 


14. 


15. 


16. 


17. 


18. 


19. 


20. 


21. 


22. 


23 


24. 


23, 


68 


MEDICATIONS 
1. ORAL 
2. INTRAMUSCULAR 
3. SUBCUTANEOUS 
4. !&STILLATSION..OF. DROPS 
5. INTRAVENOUS 
O. RECTAL 
7. JNUNCTION 


CARDINAL SYMPTOMS 
1. MOUTH TEMPERATURE 
2. RECTAL TEMPERATURE 
3. RADIAL PULSE 
4. BLOOD PRESSURE 


THERAPEUTIC MEASURES 
l. HOT WATER BOTTLE 
@. Ice BAG 6. Ultra-Sound 


3. COMPRESS 7. Diathermy 


4. HOT SOAK 
5, dE camp 8. Whirlpool 
9. Bandages 


URINE TESTS 
1. CLINITEST 
2. ACETONE TEST 
3. TESTAPE OR URISTIX 


ENEMA 
1. CLEANSING 
2. RETENTION 
3. REMOVAL OF FECAL IMPACTION 
4. INSERTION OF COLON TUBE 


CARE OF COLOSTOMY or ILELOSTOMY 
1. CHANGE OF DRESSING 
2. IRRIGATE IN BED 
3. IRRIGATE IN TOILET ROOM 


ASSOCIATED NURSING CARE FUNCTION 
1. TELEPHONE CONTACT OR REQUEST 
2. PERSONAL CONTACT 
3. ASSISTANCE WITH PROCEDURES 
4. CHARTING 


WEIGHING PATIENT 
1. AT BEDSIDE 
2. |N WEIGHING ROOM 


GAVAGING PATIENT 


FEMALE PROCEDURES 
1. VAGINAL DOUCHE 
2. PERINEAL CARE 


OXYGEN THERAPY 
1. CANNULA 
2. CATHETER 
3. MASK 
4. Te NT 


CARE OF CRITICALLY ILL 
1. ACUTE CONOITION 
2. COMATOSE CONDITION 
3. DYING CONOITION 


PREPARING BODY AFTER DEATH 


MISCELLANEOUS 
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of the existing individual patient care plan; that is, 

an indication of the potential number of occurrences of 

a care area for each patient. Since our aim here is to 
relate demand to homogeneous patient categories, per- 
formance times for a "typical" or "average" patient of 
each group would be the more appropriate measures. In 
order to derive such information, the source data for the 
Colorado study were obtained, and aggregate times by care 
area for all patients found. These times include an 
allowance for the procedural steps of each activity with- 
in a care area that were not performed but were deter- 
mined by professional judgment to be necessary. The use 
Of time information in this form is appropriate to the 
normative models to be proposed. At this point, the 
aggregates were divided by the number of patients in each 
category, and then divided again by 3 1/2 days. The 
resulting information is thus interpreted as the average 
required performance times to be devoted to each care 
area for an "average" patient of each category per 24-hour 
day. Note that we are implicitly assuming that the pat- 
tern of demands as demonstrated in the sample data is a 
reasonable representation of the pattern that would exist 
in populations of similarly classified patients. Com- 
pletely random selection of the patient sample for the 


Original study lends credence to this assumption. The 
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average required care area performance times by patient 
classification are given in Table III.6. As a matter of 
interest, the average total nursing time per patient day 
obtained by summing the care area performance times is 
.61 hours for minimum care patients, 1.55 hours for moderate 
care patients, and 2.77 hours for maximum care patients. 
Recall that in Chapter II we noted that McKnight's 
tri-level classification system may be interpreted as 
being analogous to the three levels of care denoted as 
Intermediate B, Intermediate A, and Skilled Nursing care. 
For the appropriate level of care used in our classifica- 
tion system the analogous levels as shown in Figure 
III.3 are Personal Care Facility, ICF/Supportive Nursing 
Facility and Skilled Nursing Home (ECF/ECU). The result 
we seek in keeping with the second of our classification 
system criteria is thus immediate. The implicit demand 
for patient-centered nursing services based on patient 
health status can be inferred by first classifying the 
patient into one of the three categories (Int. A, Int. B, 
Skilled), and then utilizing the average required per- 
formance time information available from Table III.6 for 
the corresponding classification. Since this time/demand 
data is for an "average" patient, the result is most 
powerful when used to determine demands by care area and 


patient classification for the residents of a large unit 
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Table III.6 
Average Required Performance Times Per Day (in Minutes) 


by Care Area and Patient Classification 
(Based on McKnight's Colorado Study Data) 


Patient Classification 


Care Area Minimum Moderate Maximum 
Y (N=64) (N=95) (N=36) 
Pr Elimination PoS 4.94 8.24 
2. Catheterization 202 1.00 5.84 
3. Daily Care Z230 LI 43.79 
4. Bathing 259 8.09 I3 84 
5. Bed Making 3.46 l 6.41 
@ Care of Skin .05 1.79 Pisao 
vou tare of Hair and Face LO 2.89 3202 
8. Oral Hygiene 92 #28 1.20 
9. Care of Nails .13 oD 35 
10. Rehabilitation and 
Recreation Be lol 6.68 
11. Food and Nourishment 8.56 17235 30.32 
12. Medications 5.99 9.27 10.45 
13. Cardinal Symptoms 61 290 290 
14. Therapeutic Measures .04 .44 1292 
15. Urine Tests 17 49 .18 
16. Enema EZ 1.00 1.00 
17. Care of Colostomy or 
Ileostomy x ZW x 
18. Associated Nursing 
Care Function 4.13 a za] 
19. Weighing Patient .04 .04 x 
20. Gavaging Patient x zur, * 
21. Female Procedures * .02 * 
22. Oxygen Therapy 12 os 10 
fre Care of Critically 111 20 db % 
24. Preparing Body After 
Death * * * 
25. Miscellaneous 5.05 5.50 OT 


* Not observed. 
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or an entire facility. The rationale is that the aver- 
aging process used in the derivation of required per- 
formance times would be replicated if a large enough 
patient population were considered at the facility level. 
The need for a demand-related classification sys- 
tem has thus been met. Note that Since we are consider- 
ing a facility oriented problem the "home care" classifi- 
cation is not considered. On the other hand, neither do 
we consider the "chronic care" category, mainly because 
of the lack of supporting time data. Nevertheless, the 
proposed system is well suited for the majority of long- 
term care facilities and is in consonence with the trend 
toward multi-level certification. In the next section 
we Shall demonstrate how the demands derived by this 
method may be allocated among available nursing person- 


nel. 


III.4 A Basic Staffing Model 

Given that we now have a means whereby patient 
needs may be translated into demand for patient-centered 
nursing activities through the intermediate steps of 
classification and aggregation, we now examine how this 
demand is related to staffing. This issue is the heart 
of the facility administrator's problem. Recall that the 


basic decision variables under the administrator's 
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control are the budget to be allocated to nursing staff 
and the number and skill level mix of nurses to be main- 
tained. An additional decision to be reached is how the 
available nursing resources may be assigned to various 
nursing functions. This latter concern is usually with- 
in the purview of the Director of Nursing. 

Decisions must be made, however, in an environ- 
ment constrained by several factors. Alluded to earlier, 
these include the following: 

1. An upper limit on funds which may be devoted 

to staffing; 


2. Legal staffing minimums and guidelines; 


3. Satisfaction of patient needs through nursing 
intervention; 
4. Limitations on the amount of nursing time 


available to devote to patient care due to 
other staff responsibilities; 
5. Assurance that patients with high priority 

nursing needs are not neglected. 
With all such restrictions in mind, decision variables 
may be fixed at those levels which act to optimize what- 
ever criterion is considered most important. If dollar 
cost of nurse staffing is chosen, for example, it is sub- 
mitted that a minimum cost solution would not be diffi- 


cult to find by using legal constraints as a basic 





75 


quideline and heuristic procedures thereafter to meet as 
many other restrictions as desired. We propose, how- 
ever, a normative approach using as a criterion the 
appropriateness of assignment of nursing personnel as a 
surrogate measure of quality of care. 

The techniques of mathematical programming are 
well-suited to model building of this type; that is, the 
optimization of a criterion function subject to certain 
constraints, both of which may be expressed mathematical- 
ly. Taking personnel budget as a parameter and staff 
numbers, skill level mix, and assignments as variables, 
the postulated model is cast into the mixed-integer lin- 


ear programming form 


MAX CX 

S.t. Ax =b 
x. > O and integer, VjenN, (TT) 
x, 2 0, VjeN, 





TỌ 


Specifically, the Basic Staffing Model (BSM) is given as 


MAX SSS CT tS pP. ; 
ir ijk -ijk Sk jkP5x (ON 
s.t. $S,N, £ Budget GETTA) 
al 
ZE Xizk  A,n,= 0, Vi (III.A) 
Jk 
Ž Xijk ae OF Nik (¿LTD 
L.<n.<U. and integer, Vi (TIEG) 
i- i-i 
ole 0, XN i,j,k (TEER) 
De Pye ak: V3,k (TIA 8) 
where 
i=1,2,...,1 indexes the various nursing personnel 
skill levels, j=1,2,...,J indexes the care areas, and 
k=1,2,...,K indexes the levels of patient classification. 


Further, the following quantities are defined: 


Baek = a constant measuring the appropriateness per unit 
time of having skill level i assigned to care area 
Ifor pk ient class k: C5 4420 We, j,k. 

BER = a dichotomous (0,1) constant, with a value of 1 


denoting the effectiveness per unit time of a care 
area/patient classification combination deemed to 
have a high relative priority, and a value of 0 


denoting a low priority combination. 
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are = a variable which indicates the amount of 
time devoted by skill level i to care area 
j and patient class k in a model solution. 
S, = a constant representing the dollar cost per 
unit time for skill level i. 

Die = a variable delineating the number of mem- 
bers of skill level i to be included in the 
staff mix. 

a. = a constant indicating the number of time 
units available for assignment to patient 
centered activities for skill level i (i.e., 
total available time less time needed for 
other nursing and personal activities). 

= a variable specifying the amount of time 
allocated by the model solution to care area 
j for patient class k. 

L., U. = upper and lower restriction constants on 
the number of personnel of skill level i 
which may be included in the model. 

m Bein = upper and lower restriction constants on the 

amount of time to be allocated to care area 
J for patient class k. 
With this information available, we examine the constitu- 


ent parts of the Basic Staffing Model to illustrate how 


each relates to the administrator's problem. 
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The objective function (III.2) indicates the 
maximization of the sum of both a time-related appropri- 
ateness measure over all possible skill level/care area/ 
classification combinations and a time-related priority 
measure. The meld of skill levels with jobs produced 
serves as a surrogate measure of quality of care in keep- 
ing with the findings of Revans [57,58]. More will be 
said on the priority aspect below. In addition, Liebman 
[38,39] demonstrated that an objective function given in 
terms of performance times produced superior assignment 
patterns to objective functions without such a feature 
when used in models of this type. 

Constraint CEL 3) simply assures that the upper 
limit on personnel budget will not be exceeded. In order 
to insure against aSSigning more time than is available 
for any skill level, the constraints (III.4) restrict 
all possible assignments of skill level i to be less than 
or equal to the time available for patient-centered 
activities. 

The constraint sets (III.5) and (III.8) when con- 
Sidered together serve a dual purpose. The first of these 
is to assure that when the assigned times of all skill 
Mewels to care areas and patient classifications are 
added, patient demands for all combinations are satisfied. 


Note, however, that the demands themselves are model 
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variables, albeit with certain minimum times which must 
be met and maximum times that cannot be exceeded. The 
maxima and minima for every care area/classification/ 
demand combination serve as a starting point for the con- 
sideration of priority of patient needs. 

That some activities are given priority over 
others in the resource constrained environment of a long- 
term care facility is a well-known fact. Little and 
Carnevali [40, p. 53] suggest that attention should be 
eeno n.n those problems which are most relevant to 
the patient's well being, keeping in mind the realities 
of staffing levels and the abilities of those who con- 
trıbute to patient care." They go on to suggest the 
existence of a hierarchy of patient needs within which 
physiological problems would be considered first, fol- 
lowed by safety, love, esteem, and self-actualization. 
Bower [10], on the other hand, postulates a hierarchy of 
needs wherein problems threatening a patient's life, 
dignity, or integrity are considered first; with problems 
which may destructively change the individual considered 
next; and those affecting the developmental growth of the 
individual following. It is emphasized, however, that 
priority setting should not negate the importance of 
those demands considered to be of lower priority, but 


rather to provide a means by which demands may be 
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satisfied efficiently. 

With these ideas in mind, we can incorporate a 
rudimentary priority system into the Basic Staffing 
Model in the following manner. Assume that we can 
dichotomize the set of care area/patient classification 
combinations into a set containing those combinations of 
relatively high priority and a set with combinations of 
lower priority. Guided by the average performance times 
as shown in Table III.5, the upper and lower bounds on 
Pak are set in a manner indicative of the priority class 
for (j,k). That is, relatively high priority combina- 
tions might have their lower bounds set at or above the 
average time, and the upper bound set at some greater 
value. Lower priority items, on the other hand, might 
have their lower bounds set below the average time, and 
the upper bound at or near the average. Such a framework 
conceptually insures that, in a situation of constrained 
resources, high priority combinations will receive an 
adequate allocation of nursing time, while lower priority 
combinations will be allocated at least their lower bound. 

Given the dichotomy of care area/patient classi- 
fication combinations and the loosely defined priority 
system it implies, a modification to the objective func- 
tion as shown in (III.2) is possible. That is, we would 


like to insure that should additional nursing resources 
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be available after satısfaction of the lower bounds, 
they would be assigned to the higher priority combina- 
tions. Thus, if we set Pa = 1] for the high priority 
combinations, the objective function would include a 
measure reflecting the importance of assigning addition- 
al time as is available to high priority combinations. 
In our rudimentary system, note that no further pri- 
ority ranking is given. Such a determination might be 
the subject of future research. The entire process 
described above is obviously highly subjective, as is 
the entire concept of priority setting. It is, however, 
believed to be a reasonable way to handle a difficult 
problem uSing the power of the mathematical model. 
Finally, constraint set (III.6) allows the 
specification of staffing level restrictions, again by 
using upper and lower bounds. Minimum coverage by skill 
level may be obtained by using the appropriate lower 
bounds. Staff-to-patient ratios may also be incorporated 
here. The upper bounds can be set to reflect the exi- 
gencies of the hiring prospects for a skill level, or 
Simply to restrict the allowable number of nurses of each 
skill level the model may assign. Note that the varia- 
bles n, have been restricted to be integers. By comput- 
ing the availabilities a, on the basis of an eight-hour 


shift, the integrality restrictions imply the assumption 
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of no part-time personnel. Although this is the 
manner in which we propose to use the model, an exten- 
sion to a part-time personnel policy is feasible. For 
example, if the a, were computed for a four-hour shift, 
the number of half-time personnel would be obtained, it 
then being an easy step to determine the full- and part- 
time staffing indicated. 

An additional and very important aspect of the 
Basic Staffing Model is its adaptability for use ina 
variety of circumstances. Parameters may be set in such 
a way that staffing and assignment patterns for any given 
unit and/or shift may be calculated, assuming the appro- 
priate information is available. Its greatest potential, 
however, 1S realized in aiding the solution of the aggre- 
gate budgeting problem. In fact, due to the nature of 
the available performance time data, it is best applied 
using a 24-hour day as the basic time frame. Taking the 
parameters as to reflect per shift availability, total 
daily staffing requirements and gross assignments of 
skill levels to care areas and patient classes can be 
derived. Information presented in this fashion allows 
for the flexible scheduling of the daily routine accord- 
ing to shift and individual personnel by the director of 
nursing in any way deemed appropriate for a particular 


facility. 
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III.5 Extensions to the Basic Staffing Model 

Assuming that the Basic-Staffing Model is an 
accurate normative form of the administrator's problem, 
an optimal solution will yield sufficient information to 
enable specification of the best feasible staffing mix; 
that 1s, subject to the assumed constraints on budget, 
patient mix, and so on. Although this optimal solution 
is important of itself, however, the model from which it 
was obtained allows of no information concerning alterna- 
tive courses of action. That such information is of con- 
sequence to the administrator should be quite clear. 

The potential impact on staffing of the internally and 
externally controlled factors alluded to in our discus- 
sion of the administrator's problem must be weighed in 
reaching an acceptable solution. Then, too, the effects 
of alterations in data obtained through subjective evalu- 
ation (e.g., priority estimates) on the solution can be 
of interest. 

Specifically, we propose to demonstrate a metho- 
dology that will provide insights into the Basic Staff- 
ing Model for the following four areas of concern: 

l. The effect of altering the total number of 

hours of patient-centered activities pro- 
vided (the "service level") on the objective 


function, staff mix, and required budget; 
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2. The impact of changes in the budget level 
on the objective function, staff mix, and 
service level; 

3. The effect on the objective function, staff 

mix, and required budget of alterations in 
the priority-related demand estimates; 

4. The impact on the objective function, staff 

mix, and required budget of changes in regu- 
lations concerning minimum staffing by skill 
level and/or staff-to-patient ratios. 

Note that we may examine three subsidiary prob- 
lems by means of the methods developed in connection with 
the first area. First, problems which are suboptimal 
with respect to the Basic Staffing Model optimal solu- 
tion (for a given set of parameters) may be generated by 
restricting service levels to be less than that attained 
by the BSM optimum. The resulting staff mix and budget 
required to support the mix for each such problem can 
assist the administrator in recognizing a suboptimal situ- 
ation in the existing facility staffing pattern. In 
addition, similar information may be helpful in planning 
a budget cut, since it would be easy to see the effect on 
total service and the objective function resulting from 
a decrease in the required budget. In this instance we 


simply interpret the result of the first subsidiary 
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problem in a different fashion. Finally, by relaxing 
the budget constraint and solving with alternative ser- 
vice levels, the effect of more funds on staff mix and 
the objective function may be examined. This allows for 
the consideration of the cost of providing increased 
care. A study of the results of the second area of con- 
cern allows the administrator to directly assess the 
impact of budgetary alterations. Both the service level 
and budgetary approach can provide valuable information 
for the specification of alternative courses of action, 
depending upon the orientation desired. 

In order to obtain information with respect to 
the service level, the Basic Staffing Model is modified 
so that we may explicitly restrict the service level to 
be less than or equal to a constant, S. The resulting 


model, called Model I, is given as follows: 
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MAS C.., PE Pa Bs. 
Er PIK LJK Ik qk* 7k 


S.t. Z S.N} < Budget 
i 


Zx o m ain, <0, Vi 
JK rK au 
= = V- 
u Jk 
sE Pi £S (TIAS!) 
rE 
ie <n, U. and integer, Vi 
i= i2 i 
Sen > 0, SViesek 
Dik € Bjx € Ujk’ SR 


All quantities in Model I are as in the BSM. The con- 
straint (III.9) is the necessary change. Obviously, we 
might solve Model I for each value of S, or budget level, 
of interest, such a tactic necessitating the solution of 
several mixed-integer linear programs. Obtaining solu- 
tions in this manner, however, can be a costly and time- 
consuming process. For this reason, we shall develop a 
methodology that produces solutions much more efficient- 
ly by exploiting the structure of the model. We defer 
consideration of this development to Chapter V. 

We next turn to the consideration of obtaining 
information related to the final two areas of concern. 
Recall that the demand for nurSing services for each care 


area and patient classification is essentially 
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determined by the model with êf jk 2 variable. Attached 
to each ao however, are upper and lower bounds deter- 
mined through priority related considerations. For a 
given mix of Int. B, Int. A and SNF patients determined 
by classification, the bounds on each of the three sets 
of care areas (indexed by k) are multiplied by the re- 
spective number of patients in each category. Thus, the 
resulting bounds yield, in the aggregate, the maximum 
and minimum demands for each j,k combination. Obviously, 
then, the third area of concern may be modeled by con- 
sidering alterations in all or a subset of these bound- 
ing restrictions, given in the BSM as (III.8) 

On the other hand, the final area of concern is 
usually related to the alterations in the minimum re- 
quired number of personnel of various skill levels. Such 
changes may be incorporated by appropriate adjustments 
in the constraint set (III.6) of the BSM. Combining both 
of these concepts, which require the examination of 
alternative configurations of upper and lower bounds, we 


propose Model II as a reasonable characterization: 
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MAX S23 c..,%.., + SSP.) Ps 
ijk ijk 1Jk ik jk E Jk 


Se ee A 
| ees 


zz Koay 7 any € 0, Vi 

Jk 

maa a S 

r; 

ZZ Pok < S ETT. 10) 

Jk 
L.+8(L.*-L.)Zn.ZU.+A(U,*-U.) and (EIN) 
i i i — i- i i i 

integer, Vi 
ak > 0, Wein, tek 


K = or ete mr . 
Likte (biy Da Bj LU sx (Us, ix? Vj ee) 
g Ae[0,1] 
In constraint sets (III.11) and (III.12) the upper and 


wer pound quantities L.*, U.*, L.,*, and U.,* are 
1 a Jk Jk 


chosen so that in varying 8 and/or on the interval 
[0,1] the desired range of consideration for the respec- 
tive bounds is achieved. Constraint (III.9) of Model I 
has again been included as (III.10) in this most general 
Harm of Model II. In this instance, we may fix S a priori 
and examine alterations in bounds at any given service 
level. 

In an analogous manner to that for Model I, the 


parameters @ and A could be fixed at values of interest, 
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and the resulting mixed-integer linear programs solved. 
Again, however, such a method would be tedious, necessi- 
tating a search for a more practical means of solution. 
The derivation of an algorithm for efficient solution 
of problems in the form of Model II will also be under- 
taken in Chapter vr 

In summary, the extensions to the BSM incorpor- 
ated into Model I and Model II capture the essence of 
the administrator's problem. Assuming that both models 
may be efficiciently solved, we have indicated how vari- 
ous interpretations of the results can produce valuable 
information both for the administrator and the director 
of nursing. In fact, the BSM need never be solved ex- 
plicitly, its results being easily obtained, for example, 
as the initial problem in either Model I or Model II. 
An example of how the entire unified approach to staffing 
and the aggregate budget process, beginning with patient 
classification and ending with examination of the results 


of the two models, will be given in Chapter VI. 


111.6 Data Requirements and Availability 

In order to properly implement Models I and II, 
data for the following parameters are required: 

1. Appropriateness measures, Cisk' 


2. Availability constants, asi 
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3. Priority-related upper and lower bounds on 


ind Unas 


demand for nursing service, Lin 4k 


de Priority imdieators, Pay! 
5. Upper and lower bounds on the staff mix by 
skill level, E. and Us; 

6. Salary costs by skill level, S.i 

7. Personnel budget allowance. 
Since items 5, 6, and 7 are likely to be highly depen- 
dent on local conditions and regulations we will not 
discuss them further here. We shall, however, indicate 
sample values for these parameters in connection with an 
example of the methodology to be presented in Chapter VI. 


The appropriateness measures, must be de- 


44k! 
termined for every combination of skill level i perform- 
ing the activities of care area j for patients of 
classification k. Obviously, such information is likely 
to be highly subjective, depending primarily on any given 
individual's perception of the appropriateness criterion. 
Having chosen to examine the issue of quality of care 
within our models, the appropriateness measures should 
reflect the relative efficiencies of assigning nurses 
with various degrees of professional expertise to the 
care of patients whose needs may differ considerably. 


Patient needs are obviously a function of both the over- 


all degree of dependence (as captured in the 
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classification) and the specific care area in question. 
The concept of effective assignment subsumes the basic 
assumption that assignment of underqualified and over- 
qualified personnel with respect to a particular care 
area/classification combination is to be avoided. 
Strictly speaking, the quantification of sub- 
jective perceptions of appropriateness should be execut- 
ed so that interval-scaled Saar are obtained for use in 
our models. One such technique applicable in this in- 
stance is the Q-sort method demonstrated by Whiting [79]. 
This system essentially requires an evaluator to rank 
order the combinations of personnel and procedure from 
most appropriate to least appropriate by placing speci- 
fied numbers of assignments into several discrete cate- 
gories. The choice of intervals and required numbers 
per interval can lead to a discretized normal distribu- 
tion, from which the interval ranking data are obtained. 
Liebman [38,39] applied this technique in deriving 
appropriateness measures for her staffing models, addi- 
tionally comparing several transformations of the inter- 
val data for their relative properties. It was found 
that the set of ordinally scaled appropriateness coef- 
ficients produced results which were more in keeping 
with the evaluator's original perceptions of effective- 


ness than results oktained using other transformations. 
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Consequently, Liebman recommended the use of ordinally 
scaled preference measures in future staffing studies. 

In order to derive such ordinal data for our 
models, the skill level preferences as reported by 
McKnight and Steorts [43] were utilized. Developed by 
a panel of three professional nurses, preferences are 
reported by patient nursing need as reflected in the 
tri-level classification (minimum, moderate, maximum 
care) for 71 specific nursing tasks. First and second 
choices of skill level usually are given. For purposes 
of completeness, the skill level category definitions 
which guided the panel are given below: 

RN: Registered Professional Nurse - one who has 
successfully completed a program in an 
approved school of nursing and who is 
licensed to practice nursing as a registered 
nurse in the state in which she is employed. 

LPN: Licensed Practical Nurse - (a) one who is 
a graduate of an approved school of prac- 
tical nursing and is duly licensed under the 
provisions of the Practical Nurse Act of the 
state in which she is employed; (b) one who 
has received a license by waiver or examina- 
tion as provided for in the Practical Nurse 


Act of the state in which she is employed. 
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NA: Nursing Assistant - one who is employed as 
an "auxiliary personnel" in the nursing 
service of a health facility to assist the 
nurse. These persons are employed and 
trained to perform tasks which involve 
specified services for patients as delegat- 
ed by the professional nurse and performed 
under the supervision of a professional 
nurse or licensed practical nurse. 

In order that compatibility between these results and 
results by care areas might be obtained, the tasks sub- 
sumed by each care area were identified, and the prefer- 
ences for each task examined collectively. Although in 
most cases the overall care area preference structure 

was apparent, professional judgment* was employed to 
indicate the structure for those cases that had no readi- 
ly apparent solution. Scoring a first preference as 2 
and a second preference as l, the resultant appropri- 
ateness measures are shown in Table III.7. 

The second set of parameters of the staffing 
models which must be obtained are the availability con- 
stants, a.. The manner of presentation of these guanti- 
ties is certainly dependent on the time frame and employ- 
ment status (full or part time) chosen for consideration. 


Recall, however, that we shall assume a 24-hour time 





* Consultation with Eleanor McKnight, R.N., M.P.H. 





Table III.7 


Ordinal Skill Level Preferences by Care Area and 
Nursing Need, Based on the Study of McKnight and Steorts [43] 


Care Area Minimum Moderate Maximum 
Elimination NA-2 NA-2 EEN 2 
EPN-1 LPN-1 RN-1 
mee catheterization RN -2 RN -2 RN -2 
Daily Care NA-2 NA-2 RN-2 
LPN-1 LPN-1 LPN-1 
4. Bathing NA-2 NA-2 RN-2 
LPN-1 LPN-1 LPN-1 
5. Bedmaking NA-2 NA-2 RN -2 
LPN-1 LPN-1 LPN-1 
6. Care of Skin NA-2 LPN -2 RN -2 
LPN-1 RN-1 LPN= 1 
7. Care of Hair and Face NA-2 NA-2 NA-2 
LPN-1 LPN-1 LPN-1 
8. Oral Hygiene NA-2 NA-2 RN -2 
LPN-1 LPN-1 MENE 
9. Care of Nails LPN-2 LBNSZ2 LPN -2 
NA-1 NA-1 NA-1 
10. Rehabilitation and NA-2 NA -2 RN-2 
Recreation LPN-1 LPN-1 LPN-1 
ll. Food and Nourishment NA-2 NA -2 NA-2 
12. Medications RN-2 RN-2 RN-2 
LPN-1 LPN-1 LPN-1 
13. Cardinal Symptoms LPN-2 LPN-2 LPN-2 
RN-1 RN-1 RN-1 
14. Therapeutic Measures RN-2 RN-2 RN-2 
LPN-1 LPN-1 LPN-1 
15. Urine Tests LPN-2 LPN-2 BEN-2 
RN-1 RN-1 RN-1 
16. Enema NA-2 NA-2 LPN -2 
MEN=1] LPN-1 RN-1 
17. Care of Colostomy or LPN -2 LPN -2 LPN-2 
I leostomy RN -1 RN- 1 RN -1 
18. Associated Nursing RN -2 RN-2 RN-2 
Care Functions LPN-1 LPN-1 LPN-1 
19. Weighing Patient NA-2 NA-2 RN-2 
LPN-1 BEN ad LPN-1 
20. Gavaging Patient RN -2 RN -2 RN -2 
2l. Female Procedures NA -2 NA-2 LPN-2 
LPN-1 LPN-1 NA-1 
22. Oxygen Therapy RN -2 RN -2 RN -2 
Care Of Critically Ill RN -2 RN-2 RN-2 
LPN-1 LPN-1 LPN-1 
24. Death of Patient NA-2 NA-2 NA-2 
25. Miscellaneous NA-2 NA-2 NA-2 
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span for the model with no part-time personnel allowed. 
The models then determine the number of nurses, by 
skill level, working an 8-hour shift, with shift specifi- 
cation for individuals being left to the discretion of 
the director of nursing. We would like to know, there- 
fore, the portion of a shift, averaged over all three 
shifts, that personnel of each skill level usually 
have available to devote to patient-centered activities. 
An indication of these availabilities was ob- 
tained by examination of the results of a nursing work 
sampling study concurrently undertaken by McKnight [42] 
in connection with her Colorado work. Data were obtained 
in four homes of various sizes related to the partition- 
ing of nursing time into five mutually exclusive activi- 


ty areas: 


l. Patient centered; 

2. Personnel centered (training, staff develop- 
ment) ; 

3. Unit centered (environment, supplies, etc.); 

4. Home centered; 


5. Other centered (personal, standby, relief). 
Patient-centered activities being more or less analogous 


to those tasks subsumed by the 25 care areas, the 
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percentage of time, averaged over all shifts, for these 
activities was obtained. When considering these data 
categorized by nursing skill level, both charge and 
staff RN's and LPN's were included. The resulting per- 
centages and availabilities per eight-hour shift are 


shown in Table III.8. Two factors were noted in the 


Table III.8 


Percentages of Time (Averaged Over Three Shifts) and 
Number of Hours Available Per 8-Hour Shift for 
Patient-Centered Activities 


Skill Level Percentage Availability/8-hours 
RN 69.9 5.6 
LPN 72.0 577 
NA 56.6 4.5 

derivation of these results. First, home size (as indi- 


cated by the number of units) had an obvious effect on 
the distribution of professional nursing time. That is, 
in the smaller facilities professional nurses typically 
gave more time to patient-centered activities and less 
to supervisory work than in the larger facilities. 
Second, a significant portion of the non-professional 
nurses' time was being given to unit and other centered 


Berıtıvies. 
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We next consider the upper and lower bounds on 
the demand variables a Recall that in addition to 
allowing flexibility in meetings patient demands, the 
limits are set with consideration being given to the pri- 
ority of certain care areas over others within and be- 
tween patient classification levels. Using the average 
daily required performance time data shown in Table 
III.6, professional assistance* was again enlisted to 
derive a sample set of limits. Adjustments in these 
bounds may be made by utilizing the properties of Model 
II, but the bounds proposed at least serve as a basis 
for examination of adjusted estimates. 

One factor in attempting to consider priorities 
as related to the bounds was the number of patients, by 
classification, who actually received care in connection 
with the activities of each care area. Assuming the 
distribution of activities among this patient sample is 
fairly representative of the distribution among patients 
of similar nursing needs, this input more or less re- 
flects the priorities existing on the units at the time 
of the origınal study. Taken together with more norma- 
tive concepts of priorities of care, these considerations 


formed part of the basis for setting bounds. 


* Consultation with Eleanor McKnight, R.N., M.P.H. 
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Table III.9 below indicates the percentage of 
patients of each category receiving care and the upper 
and lower bounds for each care area. Note that informa- 
tion for care areas 20 and 24 is omitted due to lack of 
observed activities. For other care area/classification 
combinations not observed in the original study but for 
which partial information did exist, an estimate was 
made. 

By abstracting the information contained in Table 
III.9 the average and extreme daily care times by patient 
classification may be derived. These results are pre- 
sented in Table III.10. As can be seen, the low and 
average times are nearly equal, with a definite tendency 
towards allowing for greater staffing requirements. 

Finally, we present an indication of the priority 
structure for care area/patient classification combina- 
tions in the last column of Table III.9. Recall that 
de: = l] indicates a relatively high priority, and 


Er = 0 a lower relative priority. 
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Table III.9 


Upper and Lower Bound on Bay (in minutes), Percentage of 


Patients in Each Category Receiving Care, and Priority 
Constants p., (N . =64, N 95, N_  =36) 
Ik min mod max 






































Care Area Class, Patients Average Lower Upper Pik 

1. Elimination min 14 ‚53 238 .66 0 
mod 60 4.94 Seno Soo 0 
max 72 8.24 9.00 9.50 1 

MN theterization «min Z . L .02 .».oO. 00d oO. 
mod 6 1.0 iS 1,50 0 
max 33 5.84 5.50 7.00 I 

3. Daily Care min 89 2.36 2.25 3.50 re 
mod 72 ZI 25 23,00 26.00 1 
max 89 43.79 40.00 45.00 0 
4. Bathing min 11 2.59 2:43 3.25 Ih 
mod za 8.09 7.50 9.00 1 
max 31 13,94 13.50 17.00 1 
S. Bedmaking min 55 3.46 3.00 4.00 0 
mod 81 6.11 Se 7.50 0 
max 69 6.41 6.00 8.50 1 
6. Care of Skin min 2 ‚05 05 .06 0 
mod 25 1.79 1,50 275 0 
max 78 11.10 11.00 15.00 1 
7. Care of Hair min 28 175 .60 .85 0 
and Face mod 79 2.89 2525 4.00 0 
max 94 OA 3.00 5,00 1 
8. Oral Hygiene min 9 ‚02 .01 -05 0 
mod 29 .58 .50 1.00 0 
max 56 1.20 1.00 2.50 0 
9. Care of Nails min 9 15 . 10 15 0 
mod 17 35 .30 Bas 0 
max 19 235 ou O 0 
10. Rehabilitation min 53 1.31 ZO 250 1 
and mod 69 4.61 4.50 7.50 1 
Recreation max 69 6.68 6.50 10.00 l 
11. Food and min 98 8.56 7.50 9.50 0 
Nourishment mod 100 17.35 16.50 20.00 0 
max 100 36.32 35.00 45.00 0 
12. Medications min 91 5,00 5.00 7.50 1 
mod 88 9.27 8.50 10.00 0 
max 82 10.45 10,00 14,00 1 
moe cardinal min 8 .61 50 eo 0 
Symptoms mod 17 .90 T5 1072 0 
max 28 2.90 2.50 4,50 0 
14, Therapeutic min 3 .04 03 .06 0 
Measures mod 8 „44 .40 215 0 
max 11 1.92 1.50 3.00 0 
15. Urine Tests min 2 17 115 mou 0 
mod 8 ‚49 ‚40 nS 0 
max 8 .18 SS ‚50 0 
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Table IFI.9, Cont. 


Care Area Class. % Patients Average Lower Upper Pik 
16. Enema min Z 712 .10 -29 0 
mod 8 1.00 2/5 oo 0 
max 17 1.00 who 2.00 0 

E Care of min 0 0 0 0 0 
Colostomy or mod 1 no 1275 4025 0 
Ileostomy max 0 0 125 aa 0 
18. Associated min 100 4.13 3.25 Deco g 
Nursing Care mod 100 Gill 3.195 2229 0 
Functions max 100 621 6.00 8.50 0 

19. Weighing Patient min 9 04 205 i 
mod 5 .04 03 05 0 
max 0 0 . 05 1:0 0 





20. Gavaging Patient - Omitted 








21. Female Procedures min 0 0 .01 504 0 
mod 1 02 -01 04 0 
max 0 0 01 04 0 
22. Oxygen Therapy min 2 2 05 rs 0 
mod 2 09 07 10 0 
max 6 16 ales 30 l 
23. Care of min 2 20 0 25 0 
erıtıcally Ill mod 1 Pr 10 30 1 
max 0 0 lS 50 l 

24. Death of Patient - Omitted 
25. Miscellaneous min 100 5405 4.50 AO 0 
mod 100 5.50 5.00 6.50 1 
max 100 6.50 6.50 750 I 


Table 111.10 


Average and Extreme Daily (24-hour) Patient Staffing 
Requirements (in hours) by Patient Classification 


Patient Classification Lower Bound Average Upper Bound 
Minimum 792 01 ás 
Moderate IO 1,55 1799 


Maximum 2:65 Zur Sr: 





Chapter IV. Patient Classification for Long-Term Care 


fie Introduction 

",..Patterns are the means by which we interpret 

the world." (Meisel [45, p. 1]) 

Having explored in general terms in the forego- 
ing chapter the question of patient classification in the 
long-term care setting, we now move to a more detailed 
treatment of the topic. First and foremost, we seek to 
demonstrate a methodology that will be attractive compu- 
tationally both for analytical purposes and for future 
implementation. The application of such a method to the 
best data currently available and the judicious use of 
the results obtained will lead us to a system which meets 
the criteria for an acceptable patient classification 
system established earlier. Recall that the essence of 
these guidelines was to categorize patients based on an 
assessment of their individual medical, psycho-social, 
and functioning status, and then to relate the classi- 
fication to demand for nursing services. As an addendum 
to the above, the classification system sought should be 
applicable in as efficient and least time-consuming a 
manner as possible. Having previously demonstrated how 


any tri-level system may be related to demand, we 
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concentrate here on derivation of a specific classifi- 
cation methodology. 

The results of previous studies discussed in 
Chapter II are, for the most part, lacking with respect 
to our criteria. Most are obviously unsuitable for re- 
lating the aforementioned factors to a prescribed level 
of care and thence to demand for service. In those 
studies where the relationships could be inferred, either 
new data have subsequently been obtained or the opera- 
tional characteristics of the method with respect to 
facility level use were not explored in depth. None of 
this is said by way of criticism of work that has pre- 
ceded this study; rather, it tends to demonstrate that 
no classification system will be universally acclaimed 
for all possible uses, and so we proceed. 

Initially, some general ideas on the generic 
class of problems known as pattern recognition will be 
discussed, followed by a brief overview of specific tech- 
nigues applicable in solving such problems. After demon- 
strating the derivation of an existing technique due to 
Walker and Duncan [78] and a discussion of its pertinence 
to the current study, application of the method to the 
data will be undertaken, and a readily applied procedure 
will be proposed. Comparisons with past studies and 


some recommendations for future research will conclude 
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cehe chapter. 


IV.2 The Pattern Recognition Problem 

Researchers in the behavioral, social, and 
physical sciences are often confronted with large amounts 
of data which are typically multivariate in nature. 
Fisher, Kronmal, and Diehr [23] point out that the medi- 
cal sciences also tend to generate information of this 
type, characteristically even more multivariate due to 
the nature of the field. All of these groups share the 
common goals of gaining an understanding of the underly- 
ing mechanisms of the animate or inanimate system which 
produced the data and becoming adept at analyzing new 
data in terms of historical precedents. Both goals are 
the province of the so-called pattern recognition prob- 
lem. 

In an excellent treatment of the subject of 
pattern recognition, Meisel [45] points out that the 
major aspects of any such study are the detection of 
possible regularities in the data and the synthesis of 
such information as may be available into a decision rule 
to be used in analyzing future data. This is the "recog- 
nition" aspect of the problem. Note, however, that the 
aplicacion of the rule need” not Only occur ex post 


facto, but may be made an integral part of the 
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derivation of the rule itself. This latter contingen- 
cy is called "learning and recognition concurrently," 

the former being simply, "learning before recognition." 
Both cases are shown schematically in Figures IV.1 and 


IV.2 which have been reproduced from Meisel [45, p. 6]. 


Figure IV.l 


Learning Before Recognition 


Labeled Pattern Processing Decision Rule 
Samples | 


/ 


New Pattern Decision Rule Classification 





Sample (Fixed) 


Figure IV.2 


Learning and Recognition Concurrently 


Decision Rule Classification 
(Adaptive) 


Error Detector 
GTeache rni 












Pattern Samples 
(Applied Sequentially) 


Rule 
Modifier 


The concepts of learning and recognition are identified 







by Sebestyen [64, p. 5] as being the key points of pat- 
tern recognition. The former is said to be, "...the 
estimation of the probability densities describing the 
distribution of the set of samples in the vector space," 


While the latter is, "...based on the evaluation of the 
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already learned conditional probability density of each 
class at the point in the vector space that represents 
the new input to be classified." 

Meisel goes on to point out, however, that 
typically the decision rule obtained via either of the 
above schemes is of academic interest only, primarily 
because the dimensionality of the problem (corresponding 
to the number of variables under consideration) is too 
large for practical usage. Certainly the advent of the 
digital computer has made work in such dimensions tract- 
able, but reduction of the data to a lower dimensionality 
is a pleasing concept from a computational as well as 
implementational point of view. This is especially true 
when the cost of obtaining information for further study 
or installation of the results is high. Some definitions 
will assist us in understanding the concepts presented. 
Definitions [45] 

Measurement Space: The space with dimensionality 
corresponding to the number of variables under 
consideration initially. 

Pattern Space: A finite-dimensional space of 
relatively low dimension compared to the 
measurement space. The corresponding varia- 
bles contain sufficient information to per- 


form the classification process. 
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Feature Selection: (Preprocessing): The deter- 
mination of the components of the pattern 
space. 

Pattern Classification: The derivation of a 
decision rule to classify an unlabeled sample 
in the pattern space. 

Given the modified framework above, the pattern recogni- 
tion problem may now be viewed as a combination of 

feature selection and pattern classification. Figure 1V.3 
taken from Meisel [45, p. 7] shows in flow chart form the 


revised pattern recognition process. 


Figure 1V.3 


Stages in the Derivation of the Decision Rule 


Physical System 


Measurement Space 


Preprocessino 
Preprocessina 


Reduced Pattern Space 


Decision Rule 


The process of feature selection is therefore 


central to the concept of pattern recognition. In terms 
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of the current study it is obvious that feature selec- 
tion will be an appropriate approach to solving the 
problems posed by the second of our criteria. Meisel 
points out several general attributes which an acceptable 
set of features must possess. These include the follow- 
ing: 
l. Low dimensionality; 
2. Sufficient information to yield recognition 
of patterns; 
3. Geometric consistency (preservation of order 
and distance concepts); 
4. Consistency of features throughout the 
samples. 
Even having determined a set or sets of features satis- 
fying Meisel's criteria we have yet to solve the prob- 
lem, for there almost certainly exists a dominance of 
some feature sets over others in terms of how well they 
perform the ultimate pattern recognition task. The de- 
termination of the best feature set falls within the pur- 
view of the "optimal subset (feature) selection" problem 
and will be discussed in more detail later in this chap- 
ter. 
A final concept which Meisel discusses is the 
quantification of the decision rule obtained from our 


process. In general, two methods are available. The 
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first is the specification of a boundary function such 
as 

b(x) = 0 | (MEL) 
Pieces is taken to be EneWector of features for an 
unlabeled point. The decision rule simply stated is to 
assign the point to group one if (IV.1) is non-negative 
or to group two if (IV.1) is non-positive. Points fall- 
ing on the boundary may be assigned arbitrarily. For N 
groups, the derivation of 


N(N-1) 
2 


Gin) 
such boundaries is shown to be sufficient. The second 
approach uses vector functions of the features, one for 
each of the groups, called Pp, (x), Po (X) reses Pp (x) - The 
rule is to classify the new point as a member of group 
if 

p; (x) > ps (x), Nai (AS) 
This method is most often associated with the technique 
of discriminant analysis which will be discussed below. 
It will also be shown to be a more general concept 
applicable to a number of multivariate techniques. 

Before concluding the discussion of pattern 
recognition, a few additional words on the reduction of 
dimensionality problems are in order. Most of the 


multivariate statistical techniques in general use today, 
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such as principal components analysis, factor analysis, 
discriminant analysis, cluster analysis, and multiple 
linear regression embody the concept of feature selec- 
tion as either their primary aim or a secondary method. 
More detail on these methods will be provided in the 
following section, but for now it is sufficient to say 
that there exists a number of mechanical procedures to 
deal with the problem in an efficient computational 
manner. Fisher, Kronmal, and Diehr [23], however, em- 
phasize the point that, "Since the goal is to classify, 
and not to demonstrate methodology, one need not re- 
strict the reduction of data to objective methods" [23, 
p. 381]. They further state that, "Rather than depend 
upon the sheer computational power of a computer, all 
available scientific knowledge should be used to reduce 
the task" [23, p. 382]. In short, the authors argue 
that beyond the techniques mentioned above, the applica- 
tion of intuition, expert opinion, and experience gained 
in previous studies is often the best way to extract 
features. One caveat that must be attached to their 
argument, however, is to be wary of confusing the pre- 
dictive powers of a feature set with a causal relation- 
ship to the classification. Care must be taken that this 
point is understood by those from whom intuitive informa- 


tion is solicited. In any event, the use of informed 
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Opinion is invaluable not only in direct feature ex- 
traction, but also in interpretation of results obtained 


by mechanical processes. 


IV.3 A Review of Multivariate Statistical Techniques 

We have seen that the patient classification 
problem posed earlier is of the pattern recognition form, 
specifically requiring the steps of feature extraction 
and pattern classification. The several statistical 
techniques that lend themselves either directly or in- 
directly to the solution of the problem each have certain 
attributes that make them possible candidates for our 
use. However, the necessary assumptions and restric- 
tions each possesses should make us wary of using them 
indiscriminantely, and for this reason we will examine 
each for its strengths and weaknesses vis-a-vis the cur- 
rent problem. 

The method of principal components analysis is 
an often-used approach to feature extraction. The ess- 
ence of the technique is the extraction from the set of 
original variables a subset of lower dimension that 
nevertheless explains a significant proportion of the 
total variation in the data sample. Beginning with no 
distributional assumption on the vector of variables, x, 


Lawley and Maxwell [37] show that if we take the matrix: 
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A as the variance-covariance matrix of the sample with 
eigenvalues 


Ayr Agnes A, (1V.4) 


where A is (p x p), and if we also have the associated 


eigenvectors of A as 


JE JES PRE AS GIVES) 
such that 

ui jun 0 Wi Aj (IV.6) 

mijl 


then the principal components of a new vector y can be 


written as 


Yrer Mol : (cua?) 


Assuming that 


Ely) =Ac (V:S) 
Var(y) =A 
then we may write 
p p 
a= LM pai y E Zipi (IV.9) 


where the z, are uncorrelated, implying that the y; can 
be represented as a weighted combination of orthogonal 
vectors u, . Noting that we can write 

URLS AS oy . 10) 
where U is the matrix whose columns are the eigenvectors 


AM and A. is a diagonal matrix witha, , ieee owiehe the 
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diagonal, and further assuming that thea, have been 
ordered by decreasing magnitude, the following result 
is obtained. With 
p 
tr(A) =S A. Ciel) 
ot, 
Bing a good characterization of the variance of y, then 


P k 
1= 


AT ZN, k< p (TVRZ) 
i i=l 
small implies that the k corresponding uncorrelated 


variables 


AZ e (IV.13) 


explain most of the variation in the sample, and hence 
are a good set of features to use. If we are willing 

to make the distributional assumption of multivariate 
Bermatrıey orex then tests of the hypothesis that the 
remaining p-k eigenvalues are essentially equal but non- 
zero may be made to enable the conclusion that the 
residual variation is split equally among the remaining 
variables. Although principal components analysis pro- 
vides a set of features it does not allow for the direct 
completion of the classification process, since no 
dependent variables are involved. It would, therefore, 
probably have to be coupled with another method (e.g., 
discriminant analysis) which would classify based on 


principal components only. Additionally, the method is 
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not invariant with respect to difference in scales for 
the variables of the sample, a drawback that presents 
difficulties in the case of the current study. 

A method that does not suffer from this last 
distraction is factor analysis, an extension of the con- 
cepts of principal components. According to Harman 
[29], the classical factor analytic model may be written 
as 


=a. +a. te. us, a . 
RS AE O u j=1 Aig: (IV.14) 


where n is the number of observed variables which are 
linearly described in terms of m (<n) common factors Fe 
and a unigue factor wc The u are called factor load- 
ings. The method thus represents the observed variables 
in terms of a weighted linear combination of common latent 
variables (factors) of lesser dimensionality than the 
original space, plus a specific factor accounting for 
variability which cannot be attributed to the common 
factors. In matrix terms we may write (IV.14) as 
Z=AF + € iy 215) 
where A is the matrix of factor loadings, F is the vector 
of common factors, and € an error vector which here in- 
cludes both the specific variability and any unrelia- 
bility. Now obviously one must specify in advance the 
number of factors, k, believed to give the most reason- 


able form to (IV.15), one possible approach being to 
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increase k until no appreciable alteration in the amount 
of variation explained is forthcoming. With this in 
mnd, vercan determine in (IV.15) frowa form of 

2=ANM£ +Y a iv 216) 
where, according to Morrison [46], 2 is the variance- 
covariance matrix of the sample and Y is the diagonal 
matrix with the errors Em on the diagonal. Lawley and 
Maxwell [37] recommend the use of the method of maximum 
likelihood to determine N and Y. They go on to point 
out, however, that for k?l we may not specify the matrix 
ÑN uniquely, and therein lies both the flexibility and 
difficulty with this technique. 

We would like the latent variables to allow an 
intuitively appealing interpretation in terms of the 
Original variables, with the sale representing the 
amount of variation of an original variable explained by 
the corresponding latent variable. To this end, and 
taking advantage of the non-uniqueness property off, a 
series of transformations on A, may be undertaken to 
yield a reasonable form. These transformations are 


called rotations and may be described mathematically by 
Ai +1 AiR E 


where Rey is the appropriate rotation matrix. Natural- 


1 


ly there is a certain amount of subjectivity in choos- 


ing R; 41” causing Harman [29] to speak of the art 
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involved in the method at this point. Some commonly 
accepted structures do exist, however, that seem reason- 
able from an interpretive standpoint. One of these is 
rotation to "positive manifold," which simply stated is 
the application of the linear transformation to make all 
elements of A non-negative. The rationale is that if we 
choose to use the correlation matrix in place of 2 in 
(IV.16), then we would not expect the loadings generally 
to be negative -in any reasonable A. A second possible 
transformation is rotation to "simple structure," where- 
in there are only a few non-zero elements in any row 

and column of A, and the pattern of loadings is signifi- 
cantly different for each factor. Comrey [17] points 
out that such a structure allows for more intuitive in- 
terpretation of the latent variables in terms of the 
original sample variables. A final type of transforma- 
tion iS a procrustean rotation, wherein an attempt is 
made to force the loading matrix to fit some precon- 
ceived target form developed from a subjective estima- 
tion of the pattern of loadings. Numerous computational 
schemes have been devised to perform these transforma- 
tions, the major differences in them being the ultimate 
structure sought and the orthogonality or obliqueness of 
the factors themselves. 


Certainly factor analysis is a reasonable 
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technique for feature extraction, the features being the 
latent variables. The method is highly subjective, how- 
ever, and like principal components analysis does not 
lead directly to an operational classification without 

an intermediate application of another method. In gen- 
eral, no distributional assumptions are made with regard 
to the sample, but applications of some of the advanced 
techniques suggested by Lawley and Maxwell [37] for test- 
ing hypotheses concerning the number of factors do re- 
quire such an assumption. 

A third possible technique for use in the pattern 
recognition problem is cluster analysis, a method very 
popular in biological taxonomy for a number of years. 
Meisel [45] differentiates between clustering methods 
and pattern classification on the basis of the fact that 
the samples used in the former are unlabeled, the goal 
being to determine natural groupings within the sample 
that might correspond to natural groupings in the popu- 
lation. Tryon and Bailey [73] call cluster analysis a 


"poor man's factor analysis," its methods being based 
more upon simple logic than esoteric mathematics. The 
clustering of variables into groups has an analogous 
interpretation to that of factor analysis, in that the 


clusters roughly correspond to the latent variables and 


as such require interpretation. Anderberg [2] suggests 
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that after deciding on the data to be used and variables 
to be included, a choice of clustering data units or 
variables must be made, followed by the choice of simi- 
larity measures, clustering criteria, and the number of 
clusters desired. As in principal components analysis 
some care must be taken with the scaling of variables to 
avoid nonsensical results. A usual technique is 
standardization to zero mean and unit variance of all 
variables. At this point there is a wide divergence in 
the specific algorithms available, depending largely upon 
whether one wishes to cluster variables or data units. 
Initially, we determine a convenient similarity measure 
to use. In the case of clustering of variables, the cor- 
relation matrix iS computationally attractive. It 

should be noted that the specific type of correlation 
used will be a function of the scaling of the variables 
as ratio, interval, ordinal, or nominal. Anderberg sug- 
gests methods by which diverse mixtures of scales can be 
made compatible. For data units, convenient similarity 
measures (which are the elements of a similarity matrix S) 
are Euclidean distance 


n 


PE a u ey) 


2,1/2 
de 


(139.18) 


where n is the number of data units, or the "city block" 


distance function 
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both of which satisfy the necessary condition that the 
measure be a metric. Taking our lead from Anderberg, 
the major classes of algorithms are the hierarchical and 
nonhierarchical methods, the latter being useful only 
in clustering data units. Among the hierarchical 
methods, agglomerative and devisive techniques are the 
most prominent. Agglomerative algorithms move from 
clusters of one unit each to a final stage where one 
cluster of all n units is formed. Intermediate stages 
are identified by various numbers of clusters. The de- 
visive algorithms perform in the opposite fashion. For 
the sake of illustration, a basic agglomerative algor- 
ithm due to Anderberg [2, p. 133] is presented: 

Step l: Begin with n clusters each consisting 
of exactly one entity. Let the clusters be 
labeled 1 through n. 

Step 2: Search the similarity matrix S for 
the most similar pair of clusters. Let the 
chosen clusters be labeled p and q and let 
their associated similarity be Dee, pra. 

Step 3. Reduce the number of clusters by 1 
through merger of clusters p and q. Label 


the product of the merger q and update the 
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Similarity matrix entries in order to reflect 
the revised similarities between cluster g 
and all other existing clusters. Delete the 


row and column of S pertaining to cluster p. 


Step 4. Perform steps 2 and 3 a total of n-l 


times (at which point all entitites will be 
in one cluster). At each stage record the 
identity of the clusters which are merged 
and the value of the similarity between them 
in order to have a complete record of the 


results. 


A widely used variant of this algorithm chooses 


Euclidean distance as the similarity measure and the 


centroids of the new clusters formed as the points of 


reference. 


DeviSive algorithms cause successive splits 


in the set of data units until n clusters of cardinality 


one remain. 


Little will be said about nonhierarchical 


methods, save for the fact that, in general, they allow 


for the specification in advance of the ultimate number 


of clusters desired. For additional insights into com- 


parative measures of performance and application of the 


various algorithms, the reader is referred to Gower [27] 


and Slagle, Chang, and Lee [69]. 


In retrospect, a cluster analysis on the varia- 


bles of the sample would lead to feature extraction, 
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while analysis on the data units would generate natural 
classes. Attempts to cluster data units corresponding 
to a sample of patient assessment records from the CPAI 
were made in connection with the current study. 
Although appealing, this approach was found to be unac- 
ceptable both from the standpoint of compatability of 
results and the cluster membership derived from the data 
analysis. Additionally, the problem of functionally re- 
lating the features to the classes remains as in the 
previous methodologies explored. Specifically, clear 
discrimination between the clusters obtained in terms 
of level of care of member units could not be obtained. 
This effectively rendered the results unusable for deter- 
mining classification rules. As was noted in Chapter 
II, Parker and Boyd [53] were able to achieve acceptable 
results using a superset of the-data currently under 
study by adding the centroid scores of those variables 
found to be most important from the "B/T criterion" and 
relating the sums to cluster membership. A retrospec- 
tive study on the data using this rule is not reported, 
however. Attempts in connection with the current study 
to apply such a rule were found to yield erratic and 
unacceptable results. 

A fourth branch of statistical methodology for 


Our consideration is discriminant analysis, by which one 
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may ascribe a prior dependent variable to correspond to 
Eneo ledielgassificdation of each data vector x and 
derive awillinear function of x, 

y= wx , (20) 
ime =veetonu Ww lS a weighting vector and must be deter- 
mined from the data. It can be shown in the dichotomous 
case that 

w = (só) "ta (IV.21) 
where d = Xy “X, and S is the sample variance-covariance 
matrix will yield a function (IV.20) such that the vari- 
ance between the two groups in question is maximized. 
Such a criterion is reasonable from the standpoint of 
choosing w such that maximum discrimination may be made 
between groups one and two. The method has extensions 
to the k group case where k-1 functions of the form 
(IV.20) must be specified. The classification rule which 
may be derived from discriminant analysis in the general 
ease 1S tO assign the data unit, x, to group i if the 
Beene of (1V.20) at x 1s closer to (IV.20) evaluated at 
the centroid of group i than the value of (IV.20) at the 
centroid of group j, j#i. 

Discriminant analysis thus provides an accept- 

able classification rule which may be readily applied in 


accordance with our second criterion. The question of 


feature selection is not so straightforward. One 
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approach is to consider the magnitude of the elements of 
Pec oeingeineiCative er the effect cach variable has on 
the ultimate classification. This tactic may only be 
useful up to that point at which there exists an inordi- 
nate number of variables with similar Ws. How to choose 
among them is still an open question, with the optimal 
subset selection problem also coming into play. An 
additional note on the method is that any ordering infor- 
mation concerning the dependent variables is ignored, 
thus eliminating its potential use. The results of 
Parker and Boyd [53] reported in Chapter II are never- 
theless an encouraging application of the method, albeit 
some 60 variables were used to effect the classification. 
The widely used techniques of multiple linear 
regression (MLR) have met with great success, under the 
appropriate conditions, in handling the problems of 
classification and feature selection. The basic model is 
y= xp+eé€ (IV.22) 
Mete y 1s the vector of dependent variables (n x 1), X 
the data matrix Í[nx(p+1)], f the vector of parameters to 
Brzestimated [(p+l) xl], and € the error vector (n x 1). 
The method of maximum likelihood can be applied to obtain 


estimates of f, called b such that 


b = xx) tx y (IV.23) 


the model (1V.22) thus being fitted and providing a good 
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classification tool. The theory of MLR has been suf- 


ficiently advanced to allow us to make some reasonable 


attempts at feature selection once (IV.23) has been ob- 


tained. Draper and Smith [22] point out four ways in 


which this may be done. In each case an attempt is made 


to find a smallest subset of variables without signifi- 


cantly degrading the predictive power of the model. The 


four methods are the following: 


1% 


Examine all possible combination of varia- 


a ONS „P such combinations in 
individual regression models. 

Backward elimination - Solve the model in p 
variables. Test the hypothesis that each b, 
obtained is significantly different from 

zero using an F- or t-test. Eliminate those 
found to be insignificant, resolve and repeat 
until all coefficients are found to be sig- 
nificant. 

Forward selection - Test each variable before 
entering it into the model using a partial 


F-test to see if it accounts for a signifi- 


o s e 2 o. 
cant amount of variation using R as the cri- 


terion. Stop when no such variables can be 
found. 
Stepwise regression - Similar to forward 


selection, except at each step of the process 
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all variables in the model are reexamined 
to see if they are still significant using 
the partial F-test, as well as performing 
the function in (3). Variables which do 
not meet the level of Significance are 
either dropped or not added as the case may 
be. 

With its predictive capabilities and selection 
possibilities, MLR is the most appealing method examined 
to this point. The assumptions of multivariate normality 
of the data can often be justified using Central Limit 
Theorem arguments for large sample sizes. The disturb- 
ing fact, however, is that the dependent variable is 
assumed to be a continuous, linear function of x, a 
supposition that is actually unwarranted in a situation 
where we are attempting to fit a truly polychotomous 
response. Although the ordering information on the de- 
pendent variable will now be important, something more 
than an ordinal relationship will be implicitly assumed. 
That is, typically we might assign the value 1 and 2 to 
the response variable for two groups, with the obvious 
ordinal relationship implied. But MLR takes such an 
assignment as an interval measure, when in actuality 
interval values (if they could be obtained) might be 1 


and 10, for example. Additionally, the assumed 
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eane ney ore Can often lead to "fuzzy" classifica- 
tions such as 1.4 between 1 and 2, in which case one is 
hard pressed to make the assignment rule mathematically 
tractable. We conclude, therefore, that MLR probably 
comes closest to a reasonable technique for our use, but 
its limitations for present purposes are serious enough 
to warrant further exploration. 

As a final point, we reiterate the Bayesian 
approach of Parker [52] for dichotomous independent 
variables. Recall that, from Bayes' Theorem, the basic 
classification rule is determined from 


P(S}d,)P(d,) 
p(d,|s)= era) (IV.24) 


zers] d.)P(d,) 


with dis i=1,...,m, indicating the categories and S 
being the binary vector representing a patient profile 
on the n independent variables. Further, 


p(s|a,)=P(s,]4,) ME EOS E (Lv. 2>) 


Seeerwtceits tO assign the patient wath profile S to 
category i if 


p(d,|s) > P(a,| S) jži . (IV.26) 


Parker went on to perform a feature selection procedure 
using two ad hoc measures that indicated the ability of 
the original variables to separate the groups, and his 


results were quite reasonable. The Bayesian approach is 
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certainly useful in capturing the discrete nature of the 
ordered dependent variable, and methods such as those 
proposed by Parker to deal with feature selection seem 
apropos. That the BayeSian approach could be extended 
for use with polychotomous independent variables is 
obvious, but the ordinal information in the variables 
would not be fully exploited, and the computations would 
be tedious though easy. In summary, the method seems 


Pe applied eo bDinamyex, a fact emphasized by Parker. 


IV.4 The Method of Walker and Duncan 

The foregoing section serves to point out the 
basic attributes and deficiencies of the methods sur- 
veyed. A major point that could be gleaned from the 
summaries presented is that if the strengths of discri- 
minant analysis in the classification realm and the 
techniques of feature extraction developed in conjunc- 
tion with multiple linear regression could be melded 
into a single method embodying each, a satisfactory 
answer to the pattern recognition problem currently 
under consideration might be achieved. It is submitted 
that such a method is one developed by Walker and 
Duncan [78]. 

The essential feature of this method is the 


modeling of the polychotomous response variable in terms 
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of the logistic function 


Damit Hm ey 27) 
where we immediately assume the multivariate case with 
eng the vector of observations of the independent 
variables and $ the vector of parameters to be estimated. 
Note that since we assume a polychotomous response, the 
model (IV.27) is a reasonable representation of reality. 
The graph of p vs. x, given for illustrative purposes as 
Figure IV.4, is a symmetric sigmoid curve with asymtotes 
at zero and one. This is entirely in keeping with the 
notion that the observed response is an "either-or" 
affair; the individual observed either belongs to a par- 


ticular category or does not. 


Figure IV.4 


Logistic Function 
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The logistic function has a long history of use 
in the field of bioassay along with the integrated 


normal curve, f(a), where 


x! 
F(a) = eN u exp (-5AC) du, 0 LX 200 (IV 26) 


v2r )- 

From (IV.28) the probit, P, may be defined as the normal 
deviate. We shall not dwell on the bioassay applications 
of the logistic function, but further amplification is 
provided by Ashton [4]. The impetus for the Walker and 
Duncan work was an analysis of data arising from a ten- 
year retrospective study of the medical status of 5,000 
individuals. The objective was to determine the pre- 
dictive relationship between ten variables and the inci- 
dence of various degrees of coronary heart disease. In 
a later application of the method, Talbert [72] was suc- 
cessful in modeling the stagewise development of ova in 
vitro with a polychotomous response variable. The 
authors indicate a wide range of additional applications 
in a number of diverse fields. 

A highly interesting connection, demonstrated by 
Cornfield [19], may be drawn between the form (IV.27) and 
a particular derivation of the discriminant function by 
a Bayesian approach. Considering an individual to have 
a prior probability p of being in group 1 and a (l-p) 


probability of being in group 2; and if £, (x) and f, (x) 
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are the likelihood functions for the two groups, then 
Bayes' Theorem gives the posterior probability of being 


in group 2 as 





(1-p) £, (x) 
P, = TI-p)E,(altpr (x) ae 
which yields 
= 1 
Po ER) (TV. 30) 
eee ee 
1-p f, (x) 


Under the usual multivariate normal assumptions of dis- 


criminant analysis £, (x) divided by f, (x) gives 


f, (x) Se et See eee A: en, 
x) me 2 ook My MQ EZ MoI tK'E (me) “MH i 





implying that (IV.30) is of the form 


_ 1 
P2 © Ifexpl- (arx p) O) 
which is the same as (IV.27). The major difference be- 


tween the discriminant function and the logistic function 
approach, as pointed out by Cox [20], is that the latter 
is not so highly dependent on the distributional assump- 
tion except under some very special cases. He goes on to 
state that, 
If it is operationally meaningful to consider 
Prob (y=1) and Prob (y=0) as defined for each 


fixed x and having approximately the form [of a 
linearized version of (IV.27) taking (p)/P5) 1, 
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and if y=l, 0 are the only two possibilities, 

the approach from [the linear form] seems 

preferable to that based on the model of two 

normal x- populations [discriminant analysis], 
because it assumes less. This applies whether 
the analysis is made to interpret data, or 
whether the object is to 'classify' a new indi- 

vidual. [20, p. 64] 

Note that the logistic function (IV.27) as given 
may not be easily transformed to a linear form. Sucha 
nonlinear function is called "intrinsically nonlinear" 
by Draper and Smith [22]. This has led many authors to 
use the transform mentioned by Cox, namely 

P2 

ln ==A+ x'ß (BY 33) 

SE 
often called the "log odds" or logit. The approach taken 
by Walker and Duncan is to use instead the function 
(IV.27) directly via a Taylor Series first-order approxi- 
mation, then to demonstrate a method whereby the normal 
equations may be solved by a recursive procedure. This 
represents a departure from the usual methods of fitting 
nonlinear models, which are usually based either on 
iterative schemes to solve the linearized approximation 
to (IV.23) or steepest descent algorithms that use local 
gradients as their basic device. The derivation that 
follows is for the case of a trichotomous response vari- 
able, since this is the model we will ultimately utilize, 


and is based on the derivation given in the original 


work. 
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We begin by assuming that with each of our N 
observations there is associated a response variable 


which may be captured as 


P] = 0 if the observation does not fall into 
class l, 

= 1 if the observation falls into class l, 

Py = 0 if the observation does not fall into 
class 2, 

= 1 if the observation falls into class 2, 


Pa = l - Pp, - Py 


= 9) if the observation does not fall into 
class 3, 
= 1 if the observation falls into class 3. 


Essentially we shall fit two functions simultaneously, 
corresponding to 


-1 
P,=E(p,)=£,=f (a, ,P,x)=(L+exp (-4,-X'B) (IV.34) 


and 


-] 
P]+P,=E(p,+p,)=f,=f (A, ,B,x)=(l+exp(-4,-x'ß)) (IV. 35) 


from which 


Pz = 1 - (Pı+P,) eV 236) 


can be determined. The assumption made is that 


De SE (7) 


1l 2 i 


which allows us to write £, and f, with the same 


"slope," P- and illustrates that the ordinal nature of 
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the groups is important. What Ys, group lis at a 
lower level than group 2, both of which are at a lower 
level than group 3. In the application that follows, 
this assumption is borne out by noting that a patient 
requiring Intermediate B level care requires less care 


than one needing Intermediate A care, both requiring 


less care than one needing Skilled Nursing care. The 
model for the nen observation is therefore 
Pin = Flt) Bx) +€ 1, (IV.38) 


Pin * Pan = ee) en 
We now take as a linear, first-order approximation to 
(IV.38) the first two terms of the Taylor Series expan- 
sions of each, yielding 
dEl pr, SEA, 181%,)2 


m7, ê 
=f Re re + 5 e (TV. 39) 


Pin 


- Srs Brinda, en e 
p, +p, =f (,,6,x_) +————_ +. —_————_- 
Pan Aa pn TR, SE 
where EKET and 6 are assumed values of the respec- 
tive parameters. From (1V.39) we derive so-called 

" : : ” de de 

working observations Yin and Y5p 2s 

GEA Bea. f(a, 


Lina 
s (IV.40) 
Sen ¿A Fe: e 





SS 


The partials may be computed from (IV.27) and found to 


be Jd £ gf 
sd] In Ins Sd Sn 2n 
(IV.41) 
ih: Q es Q j=1,2 
= Seer. 7 = P x J=1l,2,...,8 
5 f; ro ln an 56; Sie om JN 
where 
s = the number of variables plus two, 
Pi, and Par are as given by (1V.34) and (IV.36) 


respectively, and 

Qn and Q3, are one minus Pin and Pont respectively. 
It should be obvious at this point that since the indi- 
vidual equations of the model (IV.38) are interrelated, 
some sort of weighting corresponding to a correlation 
must be given to the data matrix X. It can be shown 
that such a weight matrix will be the inverse of the 
diagonal matrix of variance and covariances of the 


errors, whose elements are the 2 x 2 matrices: 


"_ =V le) = pe 0 Aaa 
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We may now specify the model in terms of the working 


observations y* as 


where 


I< 
+ 
| 


|D 


TI 


11 


N1 


N1 


12 


12 


N2 


N2 


(IVA ) 
(IV.44) 
(IV .45) 
Xis 
Xis 
(IV.46) 
“Ns 
X 
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11 
= 
S TN 
© on 
A, 0 
22 
A= : (IV .48) 
0 "A 
n 
where 
P1n%1n 9 
ne E (IV.49) 
n 0 P3n93n A PART EN 
in accordance with the form of (IV.40). The analogous 


equation to (IV.23) is therefore 


1 


A = 
@ = (X'AWAX) “AW'Wy*, (IV.50) 
and we may also write 
A 
v(0) = (x'awax)”? (IV.51) 


The derivation to this point has brought us to the form 
(IV.50) which, as was mentioned earlier, might be solved 
using an iterative algorithm. Walker and Duncan, how- 

ever, suggest an alternative procedure based on a recur- 
sive scheme that updates estimates of Ba = 8. and vr at 


the consideration of each observation. Without proof, 


the appropriate relationships are 
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an -1 i 
ne] nelönnn "nen ne Be 
where 
D =A x'V x A +w + Give. 53) 
n n-n n-l-nn n 
and A 
2 p SCE 
= +V Men In In (TV.54) 
>n *n-1 n-l-n n n K A 
Pin j Pon a en 
and where 
A A 
Pin and Poy are (IV.27) appropriately evalu- 


ated with the best current estimates of pP being 
used. 
This essentially completes the derivation save for some 
computational considerations. 

First of all, ıt can be seen that the scope of 
one iteration in this method encompasses the ultimate 
use of all available data, with the estimates of f hope- 
fully stabilizing within its course. A convergence 
eueck based on the final b and the initial guess bo is 
made after each such iteration, a useful form being 


oe 5 207 
Zn <§ İzi 2;...;S (IV 55) 


i 
where § is a tolerance value, say w. Additionally, 
the authors demonstrate that at any stage of an itera- 
mon the effect of the initial guesses of b and V may 


be removed by application of the following, where k denotes 
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the observation at which the correction takes place and 


starred quantities indicate the current values. 


-1 -]., -1 
Vi = N _ ya ) 
(IV.56) 
b. = V (yx tpx- yl ) 
-R k* k =k oO —O 


As a final computational feature, an initial centering 
of the data by subtraction of the arithmetic mean of each 
variable assists in obtaining orthogonality between the 


estimates of e and d,. Presumably, if strong colineari- 


1 
ty among the variables is suspected, methods which yield 
orthogonal columns of the data matrix X could be applied; 
these were not attempted in this study. 

Upon fitting the model (IV.38) by obtaining 
estimates d *, A5, and ps we would naturally desire to 
have a classification rule and some measures of the 
effectiveness of the model which are either self-contained 
or rule-dependent. Towards the end of deriving the 
rule, note that the following posterior probabilities of 


an individual X being in each of the three groups may be 


computed after fitting the model: 


pt_ = (1 + exp(-a¥ - x'@*)) 7+ (1V.57) 
E EA A = (L+exp (exar)? (IV.58) 
pz, = 1 -expladx - xp)" (IV.59) 
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We therefore propose the following simple classifica- 
tion rule: 
For an individual represented by the vector of 


observations Xa” classify as being in group i if 


* * a: TS 
pi, 7 Paw for all j#ži, j=1,2,3 (IV.60) 
where Pin arenrgiven Dy (IV.57), (LV. 59), (IV.59). 


In order to check the merit of the classification 
rule, a classification analysis procedure suggested by 
Parker [52,53] may be utilized. This essentially in- 
volves the comparison of the prior and posterior classi- 
fications, individual by individual, agreement consti- 
tuting a "hit" while disagreement implies a "miss." Note 
that the prior classification is assumed to be 


MAX ’ Dz Sy 
i Pi 


in 

and represents professional judgment as to a patient's 
appropriate level of care; and the posterior classifica- 
tion is 


MAX pt ’ =] feo 


a in 
and represents the model's prediction of a patient's 
level of care. As a means of displaying the dispersion 
of misses, a classification matrix can be formed by 
entering in the appropriate cells the number of indi- 


viduals with the given comparison of prior and posterior 


classifications. Additional information available 





139 


includes both a measure of the prediction and recogni- 
tion powers of the model in terms of the prediction and 
recognition rates, respectively. The prediction rate 
gives the percentage of patients with prior classifica- 
tion i among all patients assigned by the model to 
classification i. The recognition rate, on the other 
hand, is the percentage of patients with prior classifi- 
cation i which was asserted by the model to have the same 
posterior classification. The form of the classifica- 


tion matrix is shown in Figure IV.5. 


Figure IV.5 


Classification Matrix as Suggested by Parker 


Posterior Classification 


Recognition Rate: 


I a/a+b+c 
Prior Classification 2 e/dtetf 
3 i/gthìhti 





Prediction Rate: a/ałd+g e/bteth i/ctfti 


As an added measure, the uSual type of analysis 
of variance (ANOVA) may be performed. We based our calcu- 
lations on the linear approximations to the functions 


feli xX, e”) and £(A*,, X, g*) provided by the working 
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observations (IV.40) evaluated atajo y, G*. We denote 
the (2N x 1) vector of such approximations as y* and the 
(2N x 1) vector of errors 
ES (y-Y) (IV.61) 

where 

y is simply the (2N x 1) vector of pairs consist- 

ing of the logistic functions evaluated at 

Af, p*, and x and d 5, Br, and x: 


Recalling that p is the number of variables under con- 


sideration, we write the ANOVA in the following format: 


Figure IV.6 
ANOVA for Logistic Model 


ANOVA (Uncorrected for Mean) 


Source Degrees of Freedom Sum of Squares Mean Square F Ratio 


Regression Dia. p* 'XWy* SSR/df MSR/MSE 
Error 2N - p - 2 E'WE SSE/df 
Total 2N-1 6* 'XWy*+E WE 


Note that the ANOVA has been derived without a correction 

ene mean Of yi because the concept of a mean element 

value is obscured by the pair-wise nature of the elements 
themselves. Nevertheless, the information available 

from Figure IV.6 provides a good indication of explained 


variation. Observe, however, that due to the fact that a 
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linear approximation to the logistic is being fitted by 
eua node the yector yi is dependent upon the estimates 
derived and will not, in general, be invariant with re- 
spect to alterations in the number of variables, p, for 
example. This will have implications in computing the 
regression and total sums of squares, since in the ordi- 
nary linear regression case the quantity used is y, 
which is invariant. 

It is perhaps wise at this point to pause and 
emphasize a few general points concerning the use of re- 
gression-type models in classificatory applications. The 
first such item was mentioned earlier but bears repeat- 
ing: the relationships we seek between the independent 
variables (regressors) and the dependent variable (re- 
sponse) are of a predictive nature versus a causative 
one. In general, the functional specifications of each 
of these relationships will differ, and this is often a 
source of confusion for the unwary user of the results. 
In terms of our attempt to classify long-term care 
patients this caveat implies that the outcomes of the 
study should not be viewed as relating those factors 
deemed to be of consequence as directly causing a patient, 
of necessity, to be placed in a particular level cf care. 
Rather, the preferred interpretation is that for the 


data available the quantitatively salient factors were 





able to explain a significant portion of the variation 
caused by attempting to obtain congruence between a very 
large (although finite) number of possible vectors denot- 
ing patient status and three levels of care. It is 
therefore assumed in all such studies that the state of 
the world with regard to long-term care patient status 

is fairly regular, thus allowing the derived model to be 
used for predictive purposes, with the added caution given 
by Draper and Smith [22, p. 241] to, "...restrict pre- 
diction to the region of x-space from which the original 
data were obtained." 

An additional consideration in this regard con- 
cerns that portion of the variation mentioned above which 
the model cannot explain. A more general specification 
of (IV.38) might be as follows: 


Pry = GO, Bx 78, ,85,---, 8, 7m) (IV.62) 


EN u; Pon = g(t, Bx 76, 195,---,8, 7m) 


n 
where parameters previously used are as before, but g 
denotes a general functional form which is also para- 
metric on Os i=2,...,K which are unknown to us, and m 
is a vector of measurement errors in the regresscrs. A 
model such as (IV.62) can be used to illustrate the 
assumption that the unexplained variation is due to 
factors which are obscured to the researcher and to 


measurement error. The model applied is based upon such 





an assumption. Blalock [9] summarizes and extends these 
concepts in a concise fashion by noting that a form simi- 
lar to (IV.62), but without m, is usually acceptable in 
cases where the relative magnitude of error due to un- 
known factors is assumed to be greater than that attribut- 
able to measurement error. He goes on to point out that 
such a form is also a reasonable specification for non- 
experimentally derived data; that is, data which are ob- 
tained without the benefit of experimental design to con- 
trol for effects. Methods for isolating the effects of 
measurement errors and errors due to unknown factors 
(when both are assumed to be present) are relatively 


undeveloped. 


IV.5 The Variable Subset Selection Problem 

Near the outset of the chapter the question of 
feature extraction was emphasized as being an integral 
part of the pattern recognition process. Mention was 
also made of the problem involving the identification of 
the optimal set of features of size n, the so-called 
optimal subset selection problem. Some additional com- 
ments on each of these are in order. 

After having fitted the model using all of the 
appropriate and available data at our disposal, it be- 


comes necessary to reformulate our procedure somewhat so 
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that the classification function may be placed ina 
useful framework for the ultimate consumer. Usually 

this reformulation is characterized by a reduction in 

the number of original variables used to effect the model 
fit. There are two main reasons for this. First, if 

the marginal cost of obtaining each additional piece of 
data is relatively high it behooves us to present as low- 
cost a solution as can be obtained. The second consider- 
ation is that it is entirely possible that certain of 

the original variables might actually be obscuring the 
predictive power of the model by interacting in a nega- 
tive way with other variables. The elimination of such 
variables could therefore lead to an improved model if 
the prediction rate is the criterion used to choose be- 
tween alternatives. 

Theoretically, at least, we may postulate the 
existence of a subset of the variables of size n which 
would yield the best fit (in terms of some criterion) 
among all other subsets of size n. In addition, if we 
assume that the marginal cost of obtaining each variable 
subsequent to the first few is constant, then n could be 
interpreted to be the maximum number of variables that 
could be obtained within budgetary limits. The problem 
posed by the attempt to find such subsets is called an 


optimal subset selection problem, the criterion in the 
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more general situation usually being one determined from 
information theoretic concepts. 

For the case of ordinary multiple linear re- 
gression, Boyce, Farhi, and Weischedel [11] have proposed 
an enumerative algorithm which determines the variables 
constituting an optimal subset of size n based on the 
criterion of percentage of variance explained. This 
essentially combinatorial problem is solved by estab- 
lishing bounds to implicitly evaluate certain combina- 
tions of variables which are dominated by others. The 
search is bounded and exhaustive. Further, the authors 
claim that the method is superior to the feature selec- 
tion methods posed earlier in the chapter due to its 
ability to find the global optimum versus local optima. 
Although attractive, the algorithm does not appear to be 
readily applicable to the non-linear regression problem, 
primarily because of the dependence of the regression 
and total sums of squares on the estimated parameters via 
the vector y*. When considering problems of size n, for 
example, the total sum of squares would not be invari- 
ant, but would depend upon the fit itself and thus would 
lead to difficulties in establishing the necessary bound- 
ing relationships. 

In the absence of such refined methodology, we 


have attempted to establish a reasonable list of 
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candidate subsets for examination both by following the 
Beconmendat onggof Fisher, et al. [23] and invoking the 
technique of backward elimination from ordinary MLR. By 
Ang so we wrlimot, in all likelihood, find the opti- 
mal subsets of whatever sizes are chosen as being within 
reason, but as will be shown below the criteria do not 


appear to be highly sensitive. 


IV.6 Results 

We turn now to the specific application of the 
methodology to the problem of patient classification. 
Recalling that we initially seek a recognition of pat- 
terns in the data as they relate to the response variable, 
and then to select features in order to devise an opera- 
tional classification system, we proceed as follows. 

Making use of the available data from the Col- 
laborative Patient Assessment Instrument (CPAI) which 
was discussed in Chapter III, the following modifications 
were made. First, the original sample of 623 patient 
records was reduced by 295 through elimination of the 
data obtained by the Harvard research group. The ration- 
ale for this action was twofold. Discussions with re- 
searchers in the study which produced the CPAI indicated 
that the criteria used by the experts in judging the 


appropriate level of care for the eliminated data were 
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different from those used by the other two groups. 
Parker [53] noted a difference between the Harvard data 
and the data from the other two groups. After applica- 
tion of discriminant analysis to the CPAI data with 
judgment of "appropriate level of care" categorized by 
research groups as the dependent variable, a high de- 
gree of group variability was found between the research 
groups. This phenomenon was not, however , as pro- 
nounced between the Hopkins and HANYS data. 

A further reduction of 22 individuals was made 
to eliminate patients deemed to have an appropriate care 
level of "chronic care" or patients with incomplete 
records. Finally, on the recommendation of professional 
persons the record sets of patients with an appropriate 
care level of "Domiciliary" and "Intermediate B" were 
merged since these patients are likely to be quite simi- 
lar with respect to medical, functional, and psycho- 
social status, other factors intervening in their ulti- 
mate level of care assignment. The final sample used in 
the analysis is depicted in Table IV.1l 

Turning our attention to the set of variables 
to be used, recall that Table III.1 of Chapter III listed 
some 60 items included in the original CPAI, falling into 
the broad categories of socio-demographic data, function- 


ing status items, impairments, and medical status items. 
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Table IV.1 


Derived Sample for Classification Study 


Appropriate Level of Care 
2 Int. A SNF Total 


Sample Size 66 110 130 306 


It will also be remembered that we are seeking a post- 
assessment patient classification system whose purpose 
is to assist in predicting the demand for the level of 
nursing services in the long-term care setting. Of equal 
import is the potential value of the classification pro- 
cess for the placement of patients into an appropriate 
level of care, noting that other factors will also, of 
necessity, be considered. Since the main thrust of both 
purposes is related to what may broadly be termed a 
patient's health status, it is obvious that certain of 
the CPAI items do not and should not be considered in 
our system. Such variables fall mainly under the cate- 
gory of socio-demographics and as such were eliminated 
from the study. Other variables eliminated were those 
found to have such a low incidence of variation they 
could be considered to be constant over all data, and 
those items that were missing in a significant propor- 
tion of the records. The results of other studies cited 


in Chapter II tend to indicate, with the exception of 
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the socio-demographic variables, that the variables 
eliminated should not be of major predictive value. 
Table IV.2 lists those variables that were ultimately 
included. In terms of the broad categorization used 
earlier, we can see that items 1-15 are functioning status 
items, 16-25 impairment items, and 26-37 medical status 
items, thus representing a desirable set of variables 
for the stated purposes of our study. 

The method of Walker and Duncan described in 
Section 4 was coded in FORTRAN IV for execution on the 
DEC System 10. Using the appropriate level of care as 
the dependent variable with its three possible values of 
Intermediate B, Intermediate A, and Skilled Nursing 
Facility, a run was made using all 37 variables listed 
in Table IV.2. After completion of the calculations, an 
ANOVA and classification analysis were performed. In 
the classification analysis, "assigned level of care" 
refers to the assignment made by the model, while "appro- 
priate level of care" refers to the judgments of pro- 
fessionals as to the level of care the patient should be 
receiving. Table IV.3a shows the resulting estimates of 
Al, ado, and B by variable number, as well as the ratio of 
f: to S;, where s, is the sample standard deviation of 
the estimate. Table IV.3b gives the ANOVA, and Table 


IV.3c the classification analysis. As a further 





Table .2 


Variables Selected for Classification Study 


1. Mobility Level 19. Fractures and Dislocations 
2. Transferring 20. Joint Motion (Upper) 

3 Walking 21. Joint Motion (Lower) 

4. Wheeling 22. Joint Pain and Swelling 
5. Stair Climbing 23. Missing Limbs (Lower) 

6. Bathing 24. Paralysis/Paresis 

7. Dressing 25. Dentition 

8. Eating/Feeding 26. Cigarette Smoking 

9. Toileting 27. Anemia 

10. Bowel Function 28. Angina and/or M.I. 
11. Bladder Function 29. Arthritis 

12. Orientation 30. Cardiac Arrhythmias 
13. Communication of Needs 31. Congestive Heart Failure 
14. Behavior Pattern (Freq.) 32. Diabetes Mellitus 

15. Behavior Pattern (Type) 33. Hypertension, Essential 


(Incl. Malignant Hypertension) 
16. Sight Impairment 
34. Malignancy 
17. Hearing Impairment 
35. Mental Illness 
18. Speech Impairment 
36. Neurological Disorder(s) 


37. Chronic Respiratory Disease 


amplification of Table IV.3, note that the F-test for 
evidence of overall regression is highly significant 

with p<.001. Also, the overall recognition rate was 

67% while the overall prediction rate was about 70%. 
Notice, too, that all but about 5% of the misses occurred 
one step off the main diagonal, with a noticeable skew 
towards a higher level of assignment. The assumed model 


obviously had the most difficulty in predicting placement 





Results of Application of Walker and Duncan Method 


Varia 


R e 
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Table IV.3a 


* * j 
ble BY WER Variable 

Oa 17 
5.200 -- 18 
-0.310 -3.100 19 
-0.127 -0.693 20 
0.274 0.584 2) 
0.120 1.198 22 
-0.441 -4.413 23 
-0.439 -2.197 24 
70322 -1.609 25 
AO -2.712 26 
0.358 1.791 27 
-0.114 -0.805 28 
-0.008 -0.083 29 
-0.197 -0.657 30 
-0.048 -0.127 31 
AS OO TZ 32 
0.618 2.185 33 
-0.204 -1.018 34 

35 

36 

37 


Tabte rsp 


1% 


for 37 Variables 


* * 
BS Bils, 
0.125 0.885 
-0.153 -0.884 
E 
-0.238 -0.840 
-0.019 -0.079 
0.038 0.170 
0.223 1.287 
26185-00664 
-0.170 -1.200 
0.604 2.134 
0.723 1.417 
-0.345 -0.629 
-0.297 -0.939 
0.472 0.776 
-0,104 -0.259 
-0.424 -1.134 
0.307 0.852 
-0.279 -0.595 
-1.810 -2.419 
TES 
-1.095 -2.283 


ANOVA for Application with 37 Variables 


Source df SS MS 
on 36 3774.02 104.83 

Error zo 055.74 2.36 

Total 305 4407.76 2 


F 
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Table IV.3c 
Classification Matrix for Application with 37 Variables 
Assigned Level of Care 


Int. B Int mA SNF Total Recognition Rate 





Tot: 66 621 
Appropriate Int. 110 -627 
Level of 
Care SNF 130 -767 
Total 52 119 135 306 .6/2 
Prediction Rate Be .580 740 2700 
of "Intermediate A" patients. Some further comments on 


these results will be made later. 

Having seen the results made available by fit- 
ting the logistic function model for all variables under 
consideration, we now turn to the subset selection prob- 
lem. Using both the intuitive and statistical ration- 
ales explained earlier, a total of six subsets of the 
original 37 varıables were derived and are listed ın 
Table IV.4. The first two of these sets were formed from 
results reported by Parker and Boyd [53] in the applica- 
tion of discriminant and cluster analysis, respectively. 
Variables which were included in those studies but not 
in the present one were naturally excluded, and some 
additions were made to the second subset based on the 


first. The next two sets were determined from an 
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Table IV.4 


Subsets of the Original Variables Chosen for Analysis 


Subset I (Parker and Boyd Discriminant Analysis with modifications) 

1. Mobility Level 7. Communication 

2. Wheeling 8. Behavior (Type) 

3. Bathing 9. Fractures and Dislocations 

4. Dressing 10. Paralysis/Paresis 

5. Toileting 11. Dentition 

6. Orientation 12, Neurological Disorder(s) 
Subset II (Parker and Boyd Cluster Analysis with modifications) 

1. Mobility Level 7. Bowel Function 

2. Wheeling 8. Communication 

3. Bathing 9. Anemia 

4, Dressing 10. Mental Illness 

5, Eating/Feeding 11. Neurological Disorder(s) 

6. Toileting 12. Chronic Respiratory Disease 
Subset III (Cluster Analysis One) 

1. Mobility Level 7. Toileting 

2. Transferring 8. Bowel Function 

3. Walking 9. Bladder Function 

4. Wheeling 10. Orientation 

5. Bathing 11. Communication 

6. Eating/Feeding 12. Speech Impairment 
Subset IV (Cluster Analysis Two) 

1. Transferring 7. Toileting 

2. Walking 8. Bowel Function 

3. Wheeling 9, Bladder Function 

4. Bathing 10. Orientation 

5. Dressing 11. Communication 

6. Eating/Feeding 12. Mental Illness 
Subset V (Backward Elimination, a=.10) 

1. Mobility Level 7. Toileting 

2. Walking 8. Behavior (Type) 

3. Stair Climbing 9. Cigarettes 

4. Bathing 10. Anemia 

5. Dressing 11. Mental Illness 

6. Eating/Feeding 12. Chronic Respiratory Disease 
Subset VI (Backward Elimination, a=.05) 

1, Mobility Level 6. Behavior (Type) 

2. Stair Climbing 7. Cigarettes 

3. Bathing 8. Mental Illness 

4. Eating/Feeding 9. Chronic Respiratory 

5. Toileting Disease 
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application of Parker and Boyd's hierarchical centroid 
clustering algorithm to combined samples of Intermedi- 
ate and Skilled care patients, the variables being 
chosen by the B/T ratio criterion. The final two sub- 
sets were formed after application of a t-test of sig- 
nificance of the estimated regression coefficients in 
the 37-variable problem, one being taken at&=.1i0, the 


other ated =.05. The appropriate test statistic is: 


x - p? 
Eee a Moa) 
oe 
1 
for which we test 
O ae 
O 
HP; A 0 


Using the two-tailed t-test, under H the t-statistic 


0’ 
values are given as R $/s,; in Table IV.3a. Notice that 
all but the final subset contain 12 variables, a fact 
which was arrived at both by chance and design. The de- 
cision was made based on the fact that each subset 
originally contained about 12 variables, and the reduc- 
tion in data collection effort by two-thirds obtained by 
choosing 12 variables seemed reasonable. The 9-variable 
Subset represents an attempt to further reduce the num- 


ber of variables based on the encouraging result obtained 


by use of subset V. 





15.5 


For the sake of compactness, the results of 
fitting the model for each of the subsets are given in 
Table IV.5. Results of each ANOVA are summarized in 
Table IV.6. The ANOVA for each of the subsets shows 
evidence of overall regression to be strong, with 
p <.001 in each case. The ratio SSR/SST provides a 
gross indication of an upper bound for Ri since, in 
general, 

2 SSR - (Correction for Mean) 
R = —_—_—_— (IV .64) 


SST - (Correction for Mean) 


and since SSR < SST, then 


INS SR 2 
Ra = Sop 2 R , (14:65) 
If we were to choose R as the criterion for deciding 


which subset was best, subsets I, V, and VI might be 
worthy of further consideration, although all six yield 
approximately equal RS. Such a measure is not, however, 
in keeping with the spirit of the analysis; that is, 
derivation of a method which has reasonably high predic- 
tive validity for classification purposes. We therefore 
present the classification matrices as Tables IV.7 
through IV.12 to assist us in reaching a decision. To 
facilitate the decision process, the following measures 
of performance may be abstracted from the classification 


matrices for presentation in Table IV.13: 





Estimated B: from 


Table IV. 


Logistic Model for 


5 


Pe 


the Six Subsets of Table IV.4 


SUBSETS 
Variable I TE III IV V 
Mobility Level 00372 -0.404 -0. 360 -- -0.294 
Transferring -- -- -0.224 -0.304 -- 
Walking -- -- 0.051 -0.006 0.126 
Wheeling 0.071 0.105 0.109 0.098 -- 
Stair Climbing -- -- -- -- -0.433 
Bathing -0.664 -0.630 -0.672 -0.707 -0.591 
Dressing -0.493 -0.357 - - -0.427 -0.368 
Eating/Feeding -- -0.494 -0.534 -0.488 -0.557 
Toileting 0.557 0.491 0.362 0.432 0.428 
Bowel Function -- -0.116 -0.083 -0.097 -- 
Bladder Fcn. -- -- 00271 0.028 -- 
Orientation -0.275 -- -0.138 -0.071 -- 
Communication -0.285 -0.287 -0.315 -0.286 -- 
Behavior (Type) 0.204 -- -- -- 05,285 
Speech Impair. -- -- -0.008 -- -- 
Fractures -0.211 -- -- -- -- 
Paralysis -0.094 -- -- -- -- 
Dentition -0.202 -- -- -- -- 
Cigarettes -- -- -- -- 0.644 
Anemia -- 0.560 -- -- 0.628 
Mental Illness -- -1.376 -- -1.000 -1.812 
Neurol, Dis. -0.131 -0.243 -- -- -- 
Chron. Res. Dis. -- -0.868 -- -- -1.007 
a7 1.563 1.191 0.932 0.829 1583317 
a5 4.008 3.706 3.354 OZ 4.506 
Table IV.6 
ANOVA Summary for the Six Subsets of Table IV.4 
SUBSETS 
Quantity I II Tri IV 
MSR 246.86 205.98 216.79 2AL7 67 262855 
MSE 2.06 RS EAS, 2.18 ERE 
F 119.75 81.64 94.71 100.02 112.91 
df 11,294 11,294 11,294 11,294 11,294 
SSR 07155 2205.76 -2384.71 2394.42 2885.64 
SST 3521.6 3007.51 305771 3034.28 3568.69 
SSR/SST 88 la a78 79 .81 


VI 


353.57 
DA 
157.97 
8,297 
2828.54 
3493.31 
Neil 
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1. Overall prediction rate, 


2. Overall recognition rate, 
3. Percentage of misses > 1 step off main 
diagonal. 


With the information made available by these 
tables, certain conclusions may be reached concerning a 
choice of subset. On all but the third of the perform- 
ance measures proposed, the fifth subset is shown superior 
to all other subsets examined as well as the original 
37-variable problem. This last point tends to lend 
credence to the hypothesis of negative interaction among 
some of the variables which appears as error due to the 
specification of non-interaction in the original model. 
The non-interactive logistic function model is probably 
a better representation of reality for subsets V and VI. 
Subset VI performed about as well as the original model 
in terms of the first two measures. The original model 
does appear to give a higher percentage of "near misses" 
than the reduced forms, but certainly not so significant 
an amount as to outweigh its additional workload. Final- 
ly, it can be observed that all seven attempts at fitting 
the model tended to give a slight preponderance of misses 
wherein a patient was deemed to be properly placed in a 
higher level of care than prior judgment indicated. This 


can be seen by examining the totals for the assigned 
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Table IV.7 


Classification Matrix for Subset I 


Assigned Level of Care 





Int. B Int. A SNF Total Recognition Rate 
66 .500 
> 1. 
En 130 „730 
Total Sl 120 135 
Prediction Rate 647 S25 708 
Table 1V.8 
Classification Matrix for Subset DI 
Assigned Level of Care 
Int. B Int. A SNF Total Recognition Rate 
Int. B 66 .621 
Appropriate Int. A 110 .609 
Level of 
Care SNF 130 754 
Total 





Prediction Rate .774 .568 2726 





Table IV.9 


Classification Matrix for Subset III 


Assigned Level of Care 





Int. B Int. A SNF Total 
Int. B 66 
Appropriate Int. A 110 
Level of 
Care SNF 130 
Total 
Prediction Rate .767 .536 .710 


Table IV.10 
Classification Matrix for Subset IV 
Assigned Level of Care 


Int. B Int. A SNF Total 


Int. B 66 
Appropriate Int. A 110 
Level of 
Care SNF 130 
Total 





Prediction Rate . 708 .534 .693 


139 


Recognition Rate 
.500 
.609 


. 754 


Recognition Rate 
SS 
TES 


.746 





Table IV.11 


Classification Matrix for Subset V 


Assigned Level of Care 





Int. B Int. A SNF Total 
Int. B 66 
Appropriate Int. A 110 
Level of 
Care SNF 130 
Total 
Prediction Rate .800 .626 . 729 


Table IV.12 


Classification Matrix for Subset VI 


Assigned Level of Care 





Int. B Int. A SNF Total 
Int BB 66 
Appropriate Int. A 110 
Level of 
Care SNF 130 
Total 
Prediction Rate . 796 . 600 711 


Table IV.13 


160 


Recognition Rate 
.667 
.609 


.808 


Recognition Rate 


Summary of Performance Measures for the Six Subsets 


Measure I 11 111 1V 
(1) .624 .689 .671 .645 
(2) .598 .661 021 .611 
(3) 9.6 6.0 6.5 8.9 


.652 

. 600 

ATA 
V VI 
. 718 702 
.695 .676 


7.8 8.3 
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levels of care as compared with the totals for the 
appropriate levels of care. The net effect of this 
characteristic of the system is to make the classifica- 
tion a bit conservative when taken over the patient 
population of a facility. But since the ultimate use of 
this scheme is in predicting nursing demands, the nature 
of the errors made could lead to slight over-staffing, a 
situation perhaps more tolerable from the patient view- 
point than the reverse. This implies, of course, that 
the assumed cost of misclassification is equal whether 
the misclassification is to a higher or lower level of 
care. 

It is worthy of mention that since further dis- 
cussions with CPAI research personnel indicated the con- 
sideration of factors other than health status in adjudi- 
cating appropriate level of care (e.g., insurance 
arrangements, relatives at home, etc.), some attempt was 
made to analyze the records of those patients who missed 
on several of the runs. Only a general impression that 
the health status might be closer to the assigned versus 
appropriate level of care estimates could be obtained, 
however, due to unavailability of professional persons 
to reassess records in light of the current study goals. 

We conclude, therefore, by noting that subsets V 


and VI are the groupings of choice for consideration in 
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an operational patient classification system. Recall 
poar noe lanm is made that these are the optimal subsets 
of size 12 and 9, respectively, in terms of any cri- 
terion; rather, they represent reasonable approximations 
due to an apparent insensitivity of the model among both 
intuitively and statistically derived subsets. What now 
remains is the translation of the fruits of the analysis 
to an instrument engineered for use by professional 
nurses in a typical long-term care facility. 

-Factors to be considered in the design of the 
patient classification form include the assurance that 
the same assessment scales will be used as in the CPAI 
data, that easily followed directions for use will be 
available, that a clear format be developed,and that a 
simple way to reproduce the pattern recognition function 
served by the logistic model be provided. The first of 
these items may ca by including a table of 
assessment scores with the form, these being abstracted 
from the CPAI User's Manual [35]. The next two are im- 
plementation problems for which we propose a possible 
example. The last item is dealt with by use of easily 
followed formulae requiring addition, subtraction, multi- 
plication and three table look-ups in a table provided. 
In this last regard, the sensitivity of the model to a 


rounding of its estimated coefficients to one decimal 
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place was tested for subsets V and VI with a degradation 
of no more than 2% being observed in either the predic- 
tion or recognition rates. Hence, this more easily used 
format 1S recommended. 

Figure IV.7 below illustrates the proposed 
patient classification form based on subset V. If it is 
desired to use subset VI instead, appropriate substitu- 
tions may be made uSing the estimates provided in Table 
IV.5. “Note that the entire range of possible values of 
the logistic function exponent is represented in the 
"table of values." 

As a final note, the classification form of 
Figure IV.7 may be further modified, if space permits, 
to include the display of assessment and classification 
over time by duplicating the fourth and fifth columns 


with an appropriate notation for the date. 


IV.7 Final Remarks 

The analysis of this chapter has led us through 
the general pattern recognition problem to a specific 
model from which the essentials of the recognition prob- 
lem could be drawn, and thence to subset selection and 
the presentation of a proposed patient classification 
instrument based on medical and functioning status con- 


Siderations only. The aggregate numbers of patients 
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Fipure LV 
Patient Classification Form 
Name 1.D.# Room 


Directions } Assess the patient with respect to the items below in 
accordance with the attached rating scale. Multiply the 
score obtained by the factor indicated and sum the results. 
Find Total 1 and Total 2 as shown, then consult the 
enclosed table to obtain values for Class B, Class A, and 
Class S. The. maximum of these three numbers yields the 
Estimated Classification (Int. B, Int. A, SNF) 


= 


tem # Description Factor Assessment Factor x Assessment 


Mobility -0.3 


| sacning 1.0.6 
| tatingy/feeding [06 | 
[tottering [oe] 
| senavior pe) | os) 
| cigarertes | osl 
u 

as) oo 

T 


menia 100 
Mental Illness 
12 EnronzRes.D13. 


Total 1 = Total + 1.8 


pa = 
= O SI 


Total 


Total 2 = Total + 4.5 


Class B = Value of Total 1 


Class A = Value of Total 2 - Value of Total 1 


Class S = 1 - Value of Total 2 


1 


Class 
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Figure IV.7, continued 


Assessment Rating Scales 


(Note: HA means Human Assistance, MA means Mechanical Assistance) 
[from: Patient Classification for Long-Term Care: User's Manual ] 


Item Description 
1 Mobility Level 


2 Walking 


3 Stair Climbing 


4 Bathing 


5 Dressing 


6  Eating/Feeding 


Score 


Ww N eRe © 


EH 


Ww N- O RX y N e O NOA Ns O 


Ww Ne O 


mm & WN & © 


Interpretation 


Gets outside without either MA or HA 
Gets outside with MA but without HA 

Gets outside with HA, with or without MA 
Confined to institution, moves about 
without HA or MA 

Confined to institution and moves about 
with MA but without HA 

Confined to institution, moves about 
with HA, with or without MA 

Confined to bed or chair 


Walks without MA or HA 

Walks with MA but without HA 
Walks with HA but withoug MA 
Walks with HA and MA 

Does not walk (bed or chair) 
Does not walk (bedfast) 


Goes up and down flight of stairs without 
HA or MA 

Goes with MA but not HA 

Goes with HA but not MA 

Goes with MA and HA 

Does not perform 


Bathes self without assistance 

Bathes self with MA only 

Bathes self with HA, with or without MA 
Does not bathe self 


Dresses self without assistance 

Dresses self with MA only 

Dresses self with HA, with or without MA 
Does not dress self 


Feeds self without assistance 

Feeds self with MA only 

Feeds self with HA, with or without MA 
Does not feed self, is fed 

Does not feed self, tube fed 

Does not feed self; parenteral administra- 
tion of fluids 





Item 





10 


11 


12 
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Figure IV.7, continued 


Assessment Rating Scales, Continued 


Description 


Toileting 


Behavior (Type) 


Cigarette Smoking 


Anemia 


Mental Illness 


Chronic Respira- 
tory Disease 


Score 





NN Fe © NN Fe © 


Nro 


O 


Interpretation 


Uses toilet room without assistance 

Uses toilet room with MA only 

Uses toilet room with HA, with or without MA 
Does not use toilet room 


Appropriate 

Inappropriate, wandering or passive 
Inappropriate, abusive or aggresive 
Inappropriate, other 


Not current smoking 

Smoking less than 15 cigarettes/day 
15-24 cigarettes/day 

25 or more/day 


Not present 
Present (Aplastic, pernicious, folic 
acid, sickle cell, or combinations) 


Not present 
Present 


Not present 

Present (Asthma, chronic bronchitis, 
emphysema, removal of one lung (pneumonec- 
tomy) ) 





Total 1 or 
Total 2 


oe 
-4, 
Ae 
en. 
Ty. 
SA 
le 
BA. 
any 
=A) 
TE 
Ez. 
Day 
-3. 
Ee 
>. 
= 5% 
27 
-3. 
a 
-3. 
ER 
1M 
Er 
20% 
i 


-2 


22. 


-2 


2. 
-2. 
u 
le 
3 
ie 
ml. 
ili 
= 
28, 
=i 
los 


Ore NWR UKDNAIWNWO OFRNWARUNAANITIWUVUOKYNWEUANAIWUO DEHN WAUA~I DO 


Value 


.007 
.007 
.008 
.009 
.010 
-011 
.012 
.013 
.015 
.016 
.018 
"020 
#022 
.024 
027 
029 
092 
.036 
4059 
.043 
.047 
20527 
¿05% 
063 
. 069 
.076 
.083 
T091 
. 100 
=109 
ORS, 
Peo 
. 142 
. 154 
. 168 
152 
„198 
.214 
0251 
250 
69 


© OS O OIO OLO OF CO COC OIO Clo OC © COO Ceo OC Oe OF Oe Ce 2 CO OOOO 


Figure IV.7 continued 


Table of Values 


Total 1 or 
Total 2 
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20, 
=. 
=0, 
08 
=). 
=o 
20, 


NN NN NNNNDNNN hh 0000000000000 
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classified into the three levels of care thus generate 
demands for nursing services as approximated by the 
tri-level time estimates derived from McKnight's Col- 
orado Study [41]. Allocation and assignment of resources 
to meet these demands is the subject of the models pro- 
posed in Chapter III. 

The factors identified by the methodology used 
here seem to be in agreement with the findings of other 
researchers as reported in Chapter II with those few 
differences that exist being noted. In general, the com- 
mon items tend to be of the ADL (Activities of Daily 
Living) and psycho-social variety, with the various 
identified medically defined conditions seeming to depend 
on the sample used. In terms of comparison of the pre- 
dictive and recognitive powers of the system used here 
with previous work, only the results of Parker and Boyd 
[53] might provide a rough benchmark. With the total 
sample of 623 patients scored on the CPAI and all 60 
variables in use, a comparison of their discriminant 
analysis findings for what we have termed Int. B, Int. A, 
and SNF patients reveals an overall prediction rate of 
.73 and recognition rate of .79 with about 18% of misses 
on these groups coming at more than one step off the main 
diagonal of the classification matrix. Note, however, 


that with a total of 6 levels being fitted, it rightly 
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could be said that their model had more Opportunities to 
miss than ours, which is "locked in" for the middle 
category. Nevertheless, the 12- and 9-variable problems 
give a reasonably close result to that reported by Parker 
and Boyd so as to be encouraging. 

It should be emphasized that the patient classi- 
fication system proposed here can only be as accurate a 
reflection of patient health status as allowed by the 
data used in its derivation. Therefore, it is recom- 
mended that replication of the study be undertaken, uti- 
lizing a similar methodological approach with additional 
data as they become available. 

Our purpose has been to demonstrate a reasonably 
Simple methodology and approach to the derivation of 
nursing service demand data for use in allocation models. 
In addition, however, the classification scheme can be 
used to assist the appropriate agencies in patient level 
of care placement. That is, the results of application 
of the instrument can serve to "flag" patients who defi- 
nitely require certain levels of nursing care. This 
information, in turn, Can be combined with other inputs 
to arrive at an acceptable decision. 

The necessity for additional research in the 
area of subset selection for multivariate methods empha- 


sized by Fisher, et al. [23] is reiterated here because 
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of its obvious importance in both studies of this 


nature and a variety of other applications. 





Chapter V: Mathematical Programming Solution 
Methodology: Parametric Mixed-Integer 


Linear Programming 


Y I Introduction 
In Chapter III we noted that the examination of 
solutions of the Basic Staffing Model with alterations 
in the parameters provides valuable information for both 
the facility administrator and the director of nursing. 
Specifically, we saw that by extending the Basic Staff- 
ing Model to the forms of Model I and Model II we could 
interpret the results obtained by parametric changes in 
the two models to yield insights into the following four 
aspects of the administrator's problem: 
l. Alterations in the staff mix as a result of 
changing service levels; 
2. Alterations in the staff mix as a result of 
changing budget limits; 
3. Changes in legal staffing requirements; 
4. Alterations in the priority-related limits 
on patient demands. 


The Basic Staffing Model, 


MEE 
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Nas, 


was therefore modified in such a way that the first 


and second of the four aspects above were subsumed in 


Model I, while Model II was designed to relate to the 


remaining considerations. We restate these two models 


below for purposes of completeness as (V.2) and (V.3). 


MAX ZES C... X., +Z Èp. B. 
Er ijk ijk "x jkfjk 
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Si x.., -a.n. 0, 
4k 13k W 


DaB, 2° 
ik 
L; €N, 4 U; and ınteger, 
L., £€ re U: E 
jk = Pik Í “jk 
Span > O 


V4,k (V.2) 


Ni 
N j,k 


Wi,j,k 





193 


MAX UY Cc... X%..,+ FSp- . 
ae mee Ly PK Ex 
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We have previously indicated that solutions to 
each model with alterations in:the service level (S) 
the budget, and the upper and lower bounds via (4,0), 
respectively, could be accomplished by resolving each 
altered problem from scratch. In this chapter, however, 
we will propose two algorithms, one for each model form, 
which will allow for the efficient examination of a 
number of alternative values of S and the budget in 
Model I and the pair (4,8) in Model II. This will be 
accomplished by incorporation of the concepts of para- 


metric programming into an enumerative mixed-integer 
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linear programming algorithm. 

In general, parametric linear programming deals 
with the study of the optimal solutions to a family of 
mathematical programs related by a parametric representa- 
tion of one or more model data components. Included 
within the general heading of parametric programming are 
the class of techniques known as "postoptimality analy- 
sis." The main feature of these techniques is that the 
examination of the objective function as a function of 
the parameters takes place after some previous problem 
has been solved to optimality with a specific set of 
parameters. In this way, a minimum of effort relative 
to complete problem resolution need be expended in find- 
ing the new optima. The algorithms to be proposed in 
this chapter will gain their efficiency by incorporation 
of the ideas of postoptimality analysis. 

For review purposes, the basic ideas of para- 
metric linear programming will first be discussed, fol- 
lowed by an examination of mixed-integer linear program- 
ming (MILP) solution techniques. In conjunction with 
this latter topic, the literature survey of the para- 
metric MILP problems will be presented. Then a basic 
new parametric MILP algorithm will be given, and its 
applicability to Models I and II will be discussed. 


Finally, algorithms designed to solve the generic class 
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of problems represented by the structures of the two 
models will be stated, followed by the presentation of 
computational results for randomly generated example 
problems. We shall postpone to Chapter VI an example 
of the use of both algorithms in connection with the 


solution of the facility administrator's problem. 


V.2 Parametric Linear Programming 

The parametric MILP algorithm proposed is based 
on the MILP solution technique known as branch and bound 
(cf. Garfinkel and Neuhauser [25]). Branch and bound 
uses linear programming quite extensively, and since we 
will incorporate the parametric aspects of the new algo- 
rithm in the LP routine, we briefly review some of the 
major concepts of parametric linear programming. For the 
purposes of the exposition that follows, note that in 
general the upper and lower bounds on program variables 
could be explicitly included as model constraints. In 
this manner, the parametric representation of Model I and 
Model II can both be subsumed in the parametrization of 
the right-hand-sides (requirements vector). 

After addition of the appropriate slack and sur- 
plus variables we may write the general linear program 


(LP) as 
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MAX cx 
(LP) See (V.4) 
x20 


where the constraint matrix A is (m x n), the objective 
reacio nt vCctOnewcmaio (1 x n), the vector of vari- 


ap les x is (Ax), and the requirements vector b is 


(m x 1). For any set of basic variables Xp We have 
= =) = lite 
xX, = B b- B NX, > 0 VO) 


where p`} is the inverse of the basis matrix and is 
(m x m), N is the matrix of non-basic columns of A and 


is (m x n-m), and x the vector of non-basic variables, 


N’ 
is (n-m x 1). If B is an optimal basis, then ee 0. 
Suppose now that a change in the vector b is 
under consideration. Such a change may include altera- 
tions in the entire vector or a subset of its elements. 
In general, any alteration in the parameters of (LP) can 


affect either the optimality of a given solution or its 


feasibility. The optimality and feasibility criteria 


are: 
! , -] 
(optimalit Z.-C.=c_B a.-c.>0, A (V.6) 
Es y) sd O = : 
(feasibility) x.=B *b > 0 (v.7) 


=B 
Beyıousiyra change ın b has no effect in (V.6), but can 


possibly affect (V.7). As demonstrated by Hadley [28], 
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if we write the altered requirements vector, pr as 

b* = b + 0r (V.8) 
where for convenience we choose @6[0,1] and r to yield 
he desired range of consideration in -the -elements of b 
under study, then (V.7) becomes 


Btp* = B+ (b+er) (V.9) 


= 


=B b+ 08 tr > 


[0 


Thus, XA will be non-negative for certain values of @. 
ma fact; if B tr is non-negative, then the feasibility 
condition will be met for any 9€[0,1]. Writing the 
vector Bor as y, suppose that some elements of y are 


negative. Then based on (V.9) we can find that value of 


0, called Ou: such that some An £0 for 8> Our where 
i 


H- 





a (V.10) 


Then for 928 7 the optimality condition is satisfied, 
but the feasibility condition is violated. For any 
specific value of 676, the optimal solution can be 
found (if one exists) using the dual simplex algorithm, 
starting from the optimal basis obtained earlier. Thus, 


8 ,8, are of 


I: De > k 


interest, the feasibility condition (V.9) can be checked 


if only certain values of 6; say 9 


for these values. If the feasibility condition is vio- 


lated, we use the previous basis to reoptimize as 
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before. 

An interesting property of the type of post- 
optimality analysis under consideration here is given 
without proof in Theorem 1, first postulated by 
Dantzig [21]. 

Theorem 1: For the problem (LP), the objective 
function Z = MAX c x is a piecewise- 
linear, concave function of the 
requirements Vector b. 

We observe that such a strong result cannot be obtained 
in the case of the MILP (for further analysis of MILP 
continuity properties, see Radke [55]). 

A further consideration concerns the practicality 
of simultaneous parameterization of subsets of the 
requirements vector in terms of two or more parameters. 


For example, consider the problem posed below: 


oe 

Sec AX} = by + 8 Ly 
A,X, = b, +A TH, (V411) 
zn 
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The feasibility condition (V.9) requires that 


sg |2 -1 185 
bo Ar 


A O) 


Partitioning pē to be compatible with the row partition 


of the parametric alterations vector as 


-1 
B = ID, Do} , (V. IS 
we may write (V.12) as 
De + @D,r, +AD.r, > 0 rd”) 
= 1-1 22 ==" i 


The obvious difficulty is that each one of the above 
equations in the two unknowns @ andA makes the deriva- 
tion of a condition analogous to (V.10) difficult. This 
result has implications for our study of Model II, where 
Simultaneous parameterization of the upper and lower 
bounds is proposed. The combinatorial nature of such 
parameterizations may be resolved by fixing one parameter 
tor pt thus reducing (V.14) to the single parameter 
case as in (V.9). 

In this section we have briefly outlined some 
important results of parametric linear programming. 
These results, with modifications, will form the basis 
for the parametric branch and bound algorithms to follow. 
Note that due to the nature of Models I and II we have 


emly considered alterations in the requirements vector b. 
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There exist equally well-developed methodologies for 
dealing with the parameterization of the objective func- 
Sion COCMmuctenES vector c and thezgelements of the con- 
straint matrix A. For a detailed explanation of these 


results the reader may consult Hadley [28]. 


V.3 Parametric Integer and Mixed Integer Linear Program- 
ming 

When all or a subset of the variables of a linear 
program are constrained to be integers, there is a con- 
siderable increase in difficulty in solving the problem. 
It is well known that this is essentially due to the loss 
of the properties of the LP feasible region upon which 
the simplex algorithm is based. That is, there is no 
longer any guarantee that the optimal solution will lie 
at an extreme point of the feasible region. Obviously, 
therefore, the postoptimality analysis results applicable 
in the LP case cannot be extended to the case of the 
integer linear program (ILP) or MILP. Under some very 
special circumstances, however, the application of the 
Simplex algorithm and subsequent LP postoptimality analy- 
Sis will be appropriate. If the sufficient condition of 
total unimodularity (cf. Garfinkel and Neuhauser [25], of 
the ILP or MILP constraint matrix A is met, LP results 


may be applied to yield feasible solutions. Then, too, 
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it may be the case that A possesses some other special 
structure which allows the use of LP. 

Because of the general inappropriateness of the 
simplex method for use in the MILP, various other solu- 
tion techniques have been developed (cf. [25], [60]). 
These techniques essentially fall into three major 
classes: 

1. Cutting planes; 

2. Benders partitioning; 

3. Branch and bound. 

Based on the theory of each of these methodologies, 
several authors (26], [1], [6], [7], [59], [49]) have 
proposed extensions to deal with the questions of post- 
optimality analysis with varying degrees of success. We 
shall therefore briefly review each solution technique, 
indicating also the subsequent efforts in the realm of 
postoptimality analysis. Note that the efforts to be dis- 
cussed in connection with cutting planes and Benders par- 
titioning are concerned with closing the "duality gap" 

and the concept of "sensitivity analysis," both to be 
defined below. While not directly applicable to our work, 
the review of such efforts serves not only to give his- 
torical perspective, but also to indicate how postopti- 
mality analysis might be carried out for small changes 


in the requirements vector. We therefore consider the 
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techniques of sensitivity analysis 4S a subclass of the 


techniques of postoptimality analysis. 


V.3.1 Cutting Planes 

The cutting plane approach, primarily due to 
Gomory, is based on the philosophy of adding additional 
constraints, called "cuts," on the feasible region of the 
LP relaxation of the MILP. The cut constraint formation 
is based on information contained in any row of the sim- 
plex tableau (called the "source row") and such constraints 
are appended to the tableau. In general, the additional 
constraint will cause primal infeasibility, necessitating 
the application of the dual simplex algorithm. It is 
guaranteed that no mixed-integer feasible solution will 
be excluded by application of the cuts. For further de- 
tails on the specific methodology, the reader is referred 
to Garfinkel and Nemhauser [25] or Salkin [60]. The 
mathematics of cutting plane algorithms are quite elegant, 
leading one to believe the methods might be computation- 
ally attractive. This assumption has not been borne out 
in practice, however, where it is often the case that the 
finite number of steps needed for algorithmic termination 
are, in fact, unachievable in a practical sense. 

Because the cutting plane method is based on the 


solution of the LP relaxation of the MILP, with the 
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constraint set expanded by the addition of cuts, it is 
natural to conjecture that the resulting optimal LP 
might be useful in performing postoptimality analysis. 
One such attempt was made by Gomory and Baumol [26] for 
the ILP, but the analysis was somewhat different from 
what we have to this point called postoptimality analysis. 
The Gomory and Baumol approach is based on examination of 
a "Sensitivity analysis" of the final simplex tableau at 
termination of the cutting plane algorithm; that is, how 
sensitive is the optimal objective function value to 
marginal changes in the requirements vector. Such informa- 
tion may be obtained from an examination of the dual vari- 
ables, available as the non-basic slack variables of the 
primal LP simplex tableau or from solution of the dual 
linear program. A well-known result with respect to the 
primal and dual LPs is given without proof in Theorem 2. 
Theorem 2: If the primal and dual linear pro- 
grams possess feasible solutions x 
and y respectively, then x and y are 
also optimal to their respective problems, 
and further, 
cx = yb WELS) 


where yb is the dual objective function. 


In the case of the MILP, Theorem 2 has no direct inter- 


pretation since the dual for the MILP is not defined. 
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Nevertheless, in the solution of the LP relaxation of the 
MILP, (V.15) is true. In general, however, 
umge Sp on 
Thus the difference between the objective function values 
is called the "duality gap." If such a gap exists, the 
results of LP sensitivity analysis cannot be applied. 
The Gomory and Baumol approach is to use the dual varia- 
bles associated with the cut constraints to recompute 
the dual variables associated with the original constraints. 
By "imputing back" the effect ofthe dual vari- 
ables obtained using the cuts to the duals associated 
with the original constraints, the relation (V.15) was 
shown to be true for active (or tight) original con- 
straints only. Although in the usual LP problem inactive 
constraints have a zero dual variable associated with 
them, this property may or may not be true as a result of 
recomputation. Further, the recomputed duals are non- 
unique, depending upon the order of cut generation. Thus, 
the Gomory and Baumol approach does not completely solve 
the sensitivity analysis question. 

Some years later, Alcaly and Klevorick [1], 
recognizing the shortcomings of the Gomory and Baumol 
proposal, asserted a system of "generalized dual prices” 
which remedied the problems with inactive constraints. 


Nevertheless, the duality gap was shown to exist even 
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under this new methodology, leaving ILP sensitivity 


analysis an open question. 


V.3.2 Benders Partitioning 

The Benders partitioning method of MILP solution 
relies on the use of information obtained by iteratively 
solving two problems partitioned from the original MILP. 
Following the exposition given in Garfinkel and Nemhauser 
[25], we write the initial MILP as 


MAX C,X + CY 


s.t. A,X HA yezib (V.17) 


and integer 
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The Benders LP problem is 


MAX Cy X A 0 A,X) 
u 
x > 0 (v.18) 
Sos u A, > Cy» 
u>O, 


and is obtained by taking the dual of (V.17) with x 
fixed at some value, thus reducing the MILP to an LP. 


The second Benders problem, an ILP, is given as 
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MAX u 
t t 
S.L u 4c. x + y (b-A.x) ,Wu ET (v.19) 
a> len ee = 
De y 1 (b-A, x) ' Y y eo 


x > 0 and integer 


where T = (ue m is an extreme point of the constraint 


set of (V.18)) and Q = (y? u“ 


+ ey, om> 0 Usman extreme 
ray for some u € T). Benders' algorithm calls for the 
Seeicdtave sapplicataen OF (V.19) to determine an x to be 
used in (V.18), from which either an optimal solution is 
obtained or the sets T and/or Q are altered. Problem 
(V.19) is again solved, and so on. The algorithm can be 
shown to be finite. 

Based on the concepts of the partitioning method, 
Balas [7,8] proposed a constructive approach for solving 
the duality gap problem, thereby allowing resolution of 
the MILP sensitivity analysis question. As the initial 


step, a dual problem to (V.17) is defined as follows: 


MAX MIN c,x + u(b-A, x) 
x u 
Sry te. u A, “v=, (V.20) 
u>0 
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where Ny is the index set of x and N, the index set of 
Wess SOMmceroOn Of the MILPF(V.20) yıelds v, from 
which a vector of "subsides," s may be derived. These 
Subsidies, in turn, lead to the construction of the 
following special LP, formed from the relaxation of 
ely): 

MAX (c,+s)x + Coy 


Sree oe 


¡+A Y £b 2 


X, y >0, 


Balas demonstrates that 


(Cc, +s) x ES ub , (V.22) 


and thus the duality gap is closed. This is because an 
optimal solution to (V.17) will also be optimal to 
(V.21), and the dual variables of (V.21) are the proper 
values to use in a sensitivity analysis of the original 
MILP. The connection with Benders partitioning method 
is clear when (V.20) is examined. Assuming we fix x to 
be "admissible," (V.20) is the dual LP analogous to 
(V.18). The appropriate analog for (V.19) may also be 
derived. 

In spite of the fact that the MILP sensitivity 
analysis question appears to be answered, this method 
nevertheless has its drawbacks. Of major concern is that 
solution of (V.21) will not, in general, yield an MILP 


feasible solution. The reorientation of the objective 
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Function simply assures that the MILP solution may be an 
alternate optimum. Any change in b would be reflected 

in the solution for v, thereby changing Se na e, even 

if we could guarantee finding the proper solution, (V.21) 
may not be helpful in solving the postoptimality analysis 
problem. Attempts in connection with our study to find 
aomecriicrent Scheme for updating vy could not be found, 

any such alteration requiring a complete resolution of 
(V.20), which is tantamount to resolving the MILP from the 


outset with the parametric change incorporated. 


V.3.3 Branch and Bound 

The method of Branch and Bound is probably the 
most widely applied MILP solution due. Recognizing 
the essentially combinatorial nature of the MILP, this 
method seeks to identify those solutions which are domi- 
nated by others and to remove them from further consider- 
ation. Dominance is inferred from upper and lower 
bounds on the objective function obtained via the solu- 
tion of LP relaxations. 

If solution of the relaxation of the original 
MILP does not yield an MILP feasible solution (if it does, 
we are done), an integer-restricted variable in the 
optimal solution which did not terminate integer is 


chosen. Calling this variable xË, two new LPs are formed 
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with constraint sets 


A 
[E 


R = (( YA, x + Aay 4b; xy 20, x, >ix}] +1) 


R 


> (oy) ax + Any é b; x,y > 0, x; £ [x*]) 


This partitioning insures that no feasible solutions will 


1 or R, is chosen for 


solution, and the method proceeds in this fashion. An 


be excluded. The LP with either R 


upper bound Z. at problem i is obtained from the feasi- 
ble LP objective value. A lower bound, Z.: can be found 
from a MILP feasible solution. The maximal such lower 
bound is called Zo’ For any problem with constraint set 
Ry we may curtail further consideration of it and its 


successor problems if 


Zs = Lis or (V23) 
Za Z (V.24) 
ne 


The condition (V.23) implies that no better MILP feasible 
solution can be found, while (V.24) yields the conclusion 
that the best known MILP feasible solution is better 
(has higher objective function value) than any succeed- 
ing LP relaxation. Both conditions result from the fact 
that the successor to a given problem is more constrained 
than the problem itself. 

The choice of variables upon which to partition, 
the order of problem consideration, and the improvement 


of bounds are all subjects of algorithmic improvement and 
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will not be considered here. Again, the reader is 
referred to the texts mentioned earlier. 

The parametric methodologies we have discussed 
in connection with the cutting plane and Benders parti- 
tioning algorithms have been concerned with the question 
of obtaining sensitivity analysis results for the MILP 
by resorting to some modified LP form. Utilizing the 
branch and bound algorithm, however, two authors, Rood- 
man [59] and Nauss [49], have been able to provide in- 
Sights into the more general question of MILP postopti- 
mality analysis. 

The method of Roodman [59] is essentially a modi- 
fication of the branch and bound algorithm designed to 
extract sufficient information such that 

l. The critical values for parametric repre- 

sentation of the requirements vector and/or 
the objective function coefficients may be 
determined; and 

2. Efficient problem resolution after a specific 

parametric alteration in the requirements 
vector and/or objective function coefficients 
is possible. 
A crucial factor in the application of the method is the 
initial statement of the MILP in the form of the most 


restrictive problem to be considered. That is, for the 
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problem (V.17), c* and b are taken to be the minimum 
possible values. In addition, the MILP must be trans- 
formed so that the integer-restricted variables are 
binary, and the continuous variables are defined on the 
interval [0,1]. 

Initially, the restricted MILP is solved to 
optimality using branch and bound. The reason for 
fathoming (deleting from further consideration a problem 
and its successors) each node is noted as being due to 
one of the following three conditions: 

l. Bounds; 

2. Mmtegrality; 

3. Infeasibility. 

For conditions (1) and (2), critical values of the para- 
meters involved in the objective function and require- 
ments vector such that the fathoming condition is just 
violated are determined. For condition (3), a critical 
value for the requirements vector such that the problem 
would become feasible is found. All such information for 
every fathomed node must be stored. The minimum, taken 
over all nodes, of the first two of these values gives 
the critical values for the objective function and 
requirements vector parameters such that below these 
values the current optimal solution remains optimal. 


For a proposed parametric change, the appropriate 
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critical value is considered for each node (depending on 
the fathoming condition) to determine if it is violated. 
If so, these problems, with che change incorporated, 
must be resolved, and the branch and bound algorithm 
continued from these nodes. 

In retrospect, Roodman's method captures the 
essence of the MILP postoptimality analysis problem. 
Unfortunately, it has several limitations, among which 
are the following: 

1. The requirement for initial transformation 

of the MILP; 

2. Potentially large amounts of information to 
be stored (e.g., fathoming conditions, 
critical values, etc., for every node); 

3. LPs to be resolved must be started from 
scratch (if, for example, the optimal bases 
are retained, this greatly increases storage 

requirements); 

4. No use is made of LP parametric programming 
techniques. 

When taken together, it appears that these factors 
detract from the computational efficiency of the algo- 
mortnm. 

In a recent study, Nauss [49] considered the 


parametric integer linear programming problem (PILP) for 
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the following two cases: 


l. Parameterization over a finite number of 
points; 
2. Varying one parameter continuously over a 


specified range. 
The general (PILP) is proposed as 
MIN (et£,)x 


(PILP) s.t. (A+D,)x >b tr (V.25) 


k 
x > 0 and integer; ka, kt 

The solution of PILP is based on the branch and bound 
algorithm, modified by "problem dependent techniques" 
and "solution priorities" to be discussed below. The 
ultimate application of Nauss' analysis is to each of 
three specific PILP problems: the generalized assign- 
ment problem, the 0-1 knapsack problem, and the capaci- 
tated facility location problem. 

Initially, a rudimentary and general algorithm 


is given for the solution of (V.25), Defining 


Py A = PILP indexed on k with R; implying a given 
si 
restriction on the problem at node i (Ro is 


the original constraint set). 
The method proceeds as in the branch and bound algorithm, 
but all k* problems are, in general, considered at every 
node. If a fathoming condition is encountered for the 


m problem, consideration of this problem may be 
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deleted at successor nodes. Several open questions exist 
in connection with the PILP algorithm, among them 


l. The form of the relaxation of PL R to employ, 
E 


2. Any initial setup of PILP to take maximum 
advantage of problem structure, 


3. The branching rule to invoke, and how to 
apply it to the remaining problems at a 
node. 


No specific procedures are advanced to answer these ques- 
tions. 

Nauss proposes several problem dependent tech- 
niques to improve algorithmic efficiency for specific 
problem structures. The first of these is "problem reduc- 
tion," attained by pegging as many variables as possible 
to their optimal solution values and eliminating con- 
straints found to be dominated by others. This result is 
most powerful when the reduction is applicable to all 
problems of the PILP. "Feasibility recovery" involves 
the use of an existing optimal solution to the co prob- 
lem to generate a feasible solution (and upper bound) for 
the ria: problem. Third, "bounding problem reoptimiza- 
tion" is indicated in algorithms using LP as the relaxa- 
tion. The author proposes that LP parametric program- 
ming techniques such as those discussed in Section V.2 
might be a way to solve the relaxations of all problems 


of the PILP. Finally, the author suggests "wide range 


bounding" as a technique involving the use of the 
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Lagrangean dual problems to generate upper bounds for 
the problems of PILP. 

Solution priorities deal with the distribution 
of effort in solving the PILP; that is, how we should go 
about solving the k* problems. The "purely serial" 
approach calls for solution of one problem to optimali- 
ty, then gleaning as much information as possible from 
it to assist in solving the next problem, and so on. 
Information such as the separation of variables, upper 
and lower bounds, E Lonsy etc., can be 
helpful, Ang the course of enumeration for the next 
problem is similar to that of its predecessor (note that 
Roodman's algorithm is an example of the purely serial 
approach). The "lexicographically serial" approach is 
based on the continued solution of problem 1 to those 
nodes where a fathoming condition is encountered for it, 
then switching to problem 2, etc., until all k* prob- 
lems are considered. The final solution priority, the 
"parallel" approach, has all k* problems, each with the 
same Rir under consideration at every node. When a 
fathoming condition for problem k is met at a node, the 
problem is not considered at successors to that node. 
Ultimately, when all problems remaining at a node have 
been fathomed, the node itself is fathomed. Nauss states 


that, "the parallel approach...relies on the assumption 
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that solving a series of closely related problens at a 
given node can be accomplished relatively efficiently." 
[49, p. 38] He recommends the first two priorities over 
the parallel approach. 
An interesting result due to Noltemeier [51] and 
cited by Nauss is given as Theorem 3: 
Theorem 3: (Noltemerer) Given that c, A, and 
DC TC ESeE quantities, then the 
PILP parameterized in the require- 
ments vector only may be transformed 
into an equivalent PILP of the form 


MIN cx 


N 
Gh 
D 

|“ 


2 +t 


(Vv. 26) 
X a0 anden teger 


where 1S an integer vector con- 


k 
formable with b for k=1,2,...,K*. 

As a result of this theorem, continuous parameterization 

of the requirements vector may be reduced to the discrete 

case for the PILP. This result is not true, in general, 

in the parametric mixed-integer linear programming Case. 
The work of Nauss is a major effort in parametric 

programming. Little attention is given, however, to the 


mixed-integer case. In addition, several open questions 


pertaining to algorithmic implementation remain. These 
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questions were answered for the three problem forms 
mentioned earlier, but not for general PILP. As a final 
note, Nauss also recognized the difficult combinatorial 
problem implied by simultaneous parameterization of pro- 
gram data with different parameters, but did not recom- 


mend procedures to alleviate the difficulties. 


V.4 A Proposed MILP Parametric Programming Algorithm 
We now propose an enumerative algorithm for solv- 


ing the following parametric mixed-integer linear program 


(PMILP). 
MAX T ur 9 
sit A,x + A,y = b, (V27) 
(PMILP) x > 0 and integer 
yeno 


Kaliz aK 

We call the individual problems of (V.27)(P, In and the 

LP relaxations of these problems (RP). A parallel solu- 
tion approach will be taken. Further, we assume that as 
the index k increases, the problems of the PMILP become 
more restrictive. Note that the form of (V.27) implies 
that the continuous parameterization of b has been re- 
duced to the discrete case. In the absence of the result 


of Theorem 3, the discretization is certainly arbitrary, 
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and as such may be chosen in accordance with the problem 
setting. 
In our modified branch and bound algorithm, we 


wish to consider (Beal for all k at any node j. Let 


ex y$) be a solution to the R relaxed problem at 


J 
node j. Then if we write the constraint set of (RP, ) J 


(including additional restrictions due to branching) as 


E and that of (P) J as S$, itis immediate that 


k _k k k k k 
x or . E E . => .y . € S . r Nice 28 
(xj Yj) j xjr Yj? j : 


since, due to the nature of the LP relaxation 


px > gK 


V.29 
J = 


We can write 


2 (x, Y) = co, + Coy 30) 


from which we can infer the following relationships for 


our problem: 


k k ; k k . i 
(x*™, yk if (x*", y* is optimal to (P,)J 
2 (x? Y; ) ES v2 ) p (PL 
q XK = -oif gK = Ø 
J J 


co if (P,) J is unbounded 


Due to (V.28) and (V.29) we can obtain an upper bound, 


Z5, on 24" by solving (RP, ) J to yield 
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k k j 
22 = z(x9%, yo") it (x; Y; ) is optimal to Ree 
-k _ M , u 
23 = oo if Ts p 


co if (RP, ) 7 is unbounded 
With these concepts in mind, we now state a basic branch 
and bound algorithm for PMILP. 
Basic Algorithm 
Step l: (Initialization) Start at vertex j=l 
with 2 = 00 and zh=- 00 y k. Set e 
Y k, where sK is the initial constraint 
set for problem k, p=l, and go to step 4. 
Step 2: (Branching) If there do not exist any 


live vertices, go to step 7; otherwise 


select a live vertex j and go to step 4. 


Step 3: (Separation) Choose the lowest indexed 
problem k for which AS yo) $ s$ and 


partition sf into two new sets with 


corresponding live nodes by taking 


k k 
S341 = sen (xX lp e lo and 

1 1 

E Be): Ri E ok 
S542 EF Sj N (x “BL > en + DE 





oka. : , 
where Xp, İs some integer-constrained 
i 


variable such that en zZ eae Set 
1 i 
j=3+2 after partitioning. Choose the 


same partition Wk and go to step 2. 





Step 4: 


Step 3: 


Step 6: 
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(LP solution) Solve the set of 

LPs (RP)? WY k under consideration, say 
[ae EEV Er EEX]. 1£ H a feasible 
solution to (RP) set > =-0VYk > p, 
fathom the node, and go to step 2. In 
solving the problems (RP,) 7 for any sub- 
sequent k > p, Should an infeasible solu- 


san A 
tion be encountered set a =-aeyk > k. 
k ok 


Otherwise, set Z = = , update p as 

appropriate. 

Ciitegratrry test)” If (xe yo) is 

feasible to To Wk set a = ae 

On k k 

and z = MAX (z_, z.), fathom the node 
Oo =o" =j 

and go to step 2. If = y2“) are 


feasible to (P,) for a subset of 

l < k #4 k*, set bounds as above for 
these values of k, delete the corre- 
sponding problems from further consider- 
ation at successors of node j. 

(Bounds) If Zs £z° Vk, fathom the node 
J- and igo le step 2. If this condition 
is met for a subset of l1 < k Zk*, the 
corresponding problems may be deleted 
from further consideration at successors 


to node j. Go to step 3. 
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Step 7: (Completion) Terminate. For those 


values of k such that oe Bag 1S Cyn 


no feasible solution. Otherwise, for 


al Eso zÉ > -æ ' the solution 


which yielded the terminal value of = 


is optimal for the respective problem k. 


The following comments concerning the algorithm are use- 


PUT: 


It is convenient to choose the first avail- 

a oE enno demn backtracking In step 2. 
Identical separations are invoked Vk in 

step 3. 

The solution of the LPs at every node in 

step 4 should be accomplished in as efficient 
manner as possible. 

The increasingly restricted nature of the 
problems implied by increasing k does not 


allow for "feasibility recovery," but rapid 
detection of infeasible problems can enhance 
a “bounding problem reoptimization" approach. 


Exploration of this concept is the essence 


of Section VIS, 


Finiteness of the algorithm is assured if it is 


assumed that each element of x 1s bounded from above. 


This implies some maximum number of possible solutions 
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for each k, and hence the total number of solutions for 


all the problems put together is finite. 


V.5 The Upper Bounded Dual Simplex Algorithm 

Turning to a discussion of the efficient solution 
of the LP relaxations in the Basic Algorithm, several 
factors weigh on the decision concerning the use of 
specific techniques. These include: 

1. The desire to eliminate potentially large 
numbers of trivial upper and lower bound con- 
straints so as to reduce storage requirements 
as much as possible. 

2. The need for a method of regaining primal 
feasibility in a simplex tableau, having 
become infeasible due to a parametric change 
in the requirements vector and/or the addi- 
tion of constraints due to separation and 
branching. 

3. The expectation that the several problems 
under consideration at a node may be effi- 
ciently solved using LP parametric analysis. 

The concepts in item (1) may best be incorporated by use 
of the well known primal upper-bounded simplex algorithm 
(cf. Hadley [28], Zionts [86], Garfinkel and Nemhauser 


[25]). Conceptually, this method adds a bookkeeping 
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device to the simplex algorithm to deal with the upper 
bounds. Lower bounds may be handled by simple trans- 
formation of variables. The requirement of item (2) is 
best met by use of the dual simplex method, since the 
optimality conditions will remain satisfied. Both of 
these techniques are combined in the Upper Bounded Dual 
Simplex (UBDS) algorithm due to Wagner [76]. After exam- 
ining this method, we shall derive the rules required for 
parametric analysis. 

Because the UBDS algorithm is designed for upper- 
bounded variables only, lower bounds must be eliminated 
by a transformation of variables. Thus, after adding 


slacks and surplus variables as appropriate, we obtain 


the LP 
MAX CxX' 
sS. CMMA’ = b' (v31) 
Le xeu 
By defining 
x = xX tL 
pa~ b" ISSAL (Woe) 
47y -L 


and making the appropriate transformations, we obtain 


the LP given below as (V.33) 
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MAX cx 


A (V.33) 


Lex 2 u 
Assuming all x; Tavemanmuppers boundt-;, we define a 
"complementary variable x; such that 

x +x" =a. (V.34) 
Proceeding with the derivation as shown by Wagner, the 
system of constraints of (V.33) is written as 
n 
Ai By = b (V.35) 


a 
where P. are the columns of A. The essential substitu- 
tion to be made is as follows: If variable x, has its 
corresponding c,>0, then substitute the complementary 


variable into (V.35) and into the objective function, 
n 
z= SF C.X, (V.36) 
= 


The modified forms of (V.35) and (V.36) with, say, d 


variables having c> 0 are 


d n d 
z= >. (-c.)x" + 2 c.x. +2 c.a. (V.37) 
i=l 1 i=d+1 +? + j=1 7 1 
and 
d n d 
2 (-p.)x" E x.p=b- È piu: (V.38) 
i=l ©} ti=zgr1 A? To isi t? 


In establishing the initial tableau, the coefficients of 


the x; and x: Mm vests form the columns. An initial 
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basis is formed from the slack and/or artificial 


variables. Obtaining an initial dual feasible solu- 


tion 1S easy, Since all c, are non-positive. Thus if 
no artificials are present, Cp = 0 and 
° == . = = s = = — . 
Zj Cj CpB aj er cj > 0 (V 289 ) 


When artificials are used, the form of Cp is 


(0,...,0,-M,...,-M), where M is a large positive 
constant. It is not, therefore true, in general, that 


Zj > 0. In this case, for a non-basic j such that 


en 0, the corresponding column (including the entry 


FOr 23703) 1s transformed by multiplying by (-1), thus 


assuring 2,0; > 0. This is equivalent to the substi- 
tution of the complementary variable. The UBDS algorithm 
as stated by Wagner is given below: 


Step 1: (Exit rule) Choose the variable to 


leave the basis from 


x = MIN SE eee 
B oe 
r A eaa Oe = 
If all xX, are non-negative, stop. An 
i 


optimal solution has been obtained. 
Otherwise, go to step 2. 

Step 2: (Enter rule) Choose the variable to 
enter the basis from 


x, = MIN ees 1 Xo. 0 
Tale, 2. . apn J J 13 rJ 
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1f Ax, ~ 0, stop. The dual is un- 
bounded which implies the primal prob- 
lem has no feasible solution. Otherwise, 
go to step 3. 

Step 3: (Initial new tableau) Perform a trans- 
formation on the current tableau uSing 
the usual pivoting rules of the simplex 
method. Go to step 4. 

Step 4: (Final new tableau) For any basic vari- 


able, say Xp! which violates its upper 


t 
bound: 


l. Change the signs of all elements in 


the tableau row corresponding to Xp 
t 


except the coefficient 1 in the 


column denoting x 


B 
t 
2. Replace Xp by its complementary vari- 
E 

able. 

3. Label x as x! (or vice-versa as 
B B 
e t 

the case may be). 
4. Change the sign of the Ce correspond- 

ing to Xp Go to step l. 

E 


Upon algorithmic termination, the solution may 
be decoded by first noting that variables will belong to 
one of the following four mutually exclusive and collec- 


tively exhaustive classes: 
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l. Basic X; => x; = x; 
2. Non-basic x. => x. = 0 
1 1 
: we _ = _ u 
3. Basic Xx: > x; M: x; 
4. Non-basic x! => x. = e, 
i i 1 


Thus, aside from some rather easy transformations and 
decoding rules, the algorithm is much like the ordinary 
dual simplex algorithm. 

With the UBDS algorithm providing an efficient 
method for solving the general LP relaxation encountered 
at any node, we turn to the consideration of solving the 
closely related relaxed problems. This implies the need 
for a postoptimality analysis technique based on an opti- 
mal UBDS tableau. As demonstrated earlier in Section 
V.2, the basis for the total problem is necessary for sen- 
sitivity and/or parametric analysis. For this reason we 
initially examine the structure of upper bounded problens 
in a manner suggested by Hadley [28]. 

Beginning with the LP form (V.33), the so-called 
"expanded system" of constraints with upper bounds and 
complementary variables explicitly included is 


a o] fz b 
= X, x" > 0 r (V.40) 
w E 
where Ia is the identity matrix of order n and A is 


(m+n) x (mtn). 
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Hadley demonstrates the following three results. 
Result l1: A basis for the system (V.40) does 
not contain more than m variables 
which are not at their upper bounds. 
Proof: Assume there exist k>m such variables. 
Tnen there also exist k corresponding 
complementary variables in the basis. 
Also, there will be n-k additional vectors 
in the basis. Therefore, the basis will 
consist of n-k+k+k=n+k>n+m variables, 
which is impossible. QED 
Result 2: The basis matrix for (V.40) may be 


constructed in partitioned form as 


B O O R 
I I O O 
m m 
O O (v.41) 
O O O m 


where B, is (m+tn)x(m+tn), B is (mxm) 
and is the basis for (V.33) consisting 
of columns corresponding to the r 
variables not at their upper bounds. 
If r < m, the remaining m-r columns 
are filled with the columns corre- 


sponding to variables at their upper 
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bounds. Columns mtl to 2m contain 

the complementary vectors for the 
variables in columns 1 through m. The 
final k columns consist of the columns 
of remaining variables at their upper 
bounds, if any, with n-m-k being the 
remaining columns, if any. 

The inverse of the basis matrix By, 


may be written in partitioned form as 


Br? O O -BİR 
_pt I O BR 
m 
(A2) 
O O RO 
O O O ive 


Results 1-3 will be useful when attempting to 


determine specific rules for postoptimality analysis. 


Further, we may delineate two instances where AS will 


take on a special form. These are 


Case I: 


0 
Br eG 
p I 0 (V.43) 
m 
0 0 ı 
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Case II: n-m-k=0 
B+ 0 BER 
u er I BR (V.44) 
0 0 Tcl 
We propose Result 4 below. 
Result 4: Cases I and II are mutually exclusive. 


Proof: 


Without loss of generality, assume nen, 
Since n=m implies a unique solution and 
n<m implies Seneca or inconsistent 
constraints. Assume the converse of the 
proposition. Then it is impossible to 
construct a basis of dimension (m+n), 
Since the last row of (V.43) or (V.44) 
would consist of all zeroes, hence the 


result. QED. 


The remaining open question concerns the con- 


struction of B31 


from the optimal basis, D, provided by 


solution by the UBDS algorithm. Recall that at no time 


may a variable and its complement appear in the basis 


simultaneously. 


Result 5: B, 1 


as follows: 


Result 5 below provides the connection. 


may be constructed from an UBDS tableau 
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]. @Form B 1 from D by multiplying the elements 
of any row of De corresponding to a basic 
complementary variable by (-l). 

2». Form BR from UBDS tableau columns cor- 
responding to non-basic complementary varia- 
bles. 

3. The remaining n-k columns are determined 
according to results 2 and 3. 

Proof: (1) is simply the application of step 4 
of the UBDS algorithm to obtain the com- 
pact basis inverse in terms of the 
"original" variables x. (2) is immediate 
by construction of the partition in 
Result 2. (3) is obvious by construction 


of Ba 


V.6 Postoptimality Analysis of the Requirements Vector 
We may now proceed with the development of the 
rules for postoptimality analysis of the requirements 


vector. The general form of the problem to be solved is 





MAX cx 


Lasek <b 
a a Oland integer, JéN, 
x > oO, JEN, 


2.102 


(V.45) 


We assume, without loss of generality, that the upper 


and lower bounds on a jeN, are integer. After appli- 


cation of the transformation of variables in (V.32) to 


remove the lower bounds, where y = x - L, we obtain the 


following system of constraints: 

By Sar ea 

y Pruna integer", JEN, 
0, JEN 


he 
IV 


where 


[O 
D 
En 

Re, 


[cS 
l 

qn 
te) 


— 


amarre, the vector of alterations is such that 


The parameter 8 is chosen such that 86€[(0,1]. 


(V.46) 


(V.47) 


(V.48) 
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After suitably partitioning the rows of g, as 
(g5 2 95)” to conform with the most general case, 
Cy .4 Of Dar the feasibility condition (V.7) is 
BL (g+0r) >0 , (V. 49) 


which yields the following four relationships: 
1l il 


Lo ISO = BERA) + 0B r >0 (V.50.1) 

2. [gl - (Bog, - BtRgs)] - eB tr > 0 (V.50.2) 
2 

ar (V.50.3) 
3 

4. g, >0 . (V.50.4) 


Note that (V.50.3) and (V.50.4) simply require that 

ble upper «bounds.*U be greater than or equal to L. Assum- 
ing this to be the case, conditions (V.50.1) and (V.50.2) 
yield the following rules for determining the critical 


value of 80: 


ie eS 


(B g,). - (B “Rg). k 
Rule I: gt = MIN ld ee , (B Ly e 0 
Cc = ae 
- (B “r). 
zu 
(Wes) 
ih -] -]_ 3 
Rule II: g° = MIN BE. en. u nr... (Br). 7 0 
Cc -] — i 
(B r). 
— i 
l 2 


Rule III: 8 =-MIN (8,989), 
E C E 


where C); denotes the pe element of the vector result- 
ing from the matrix operation in parentheses. Computa- 


tionally, if for a given Or 6. > 8. the current basis 
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is not optimal, so reoptimization by the UBDS algorithm 
for this value of 0 is necessary. For discrete 8, the 
process is obviously finite, and for continuous 8 it is 
also finite since there exist only a finite number of 
bases. 

Rules I through III assist usin solving the LP 
relaxation by postoptimality analysis at any given node. 
Note that the additional constraints due to separation 
and branching to be appended to (V.46) may be written as 

eE (V.52) 
mere PF and Pare the bounds on y due to branching. At 
mie eproblem Outset ;soealeand E= U calefor all variables, 
and remain so throughout the algorithm for the continuous 
variables. Due to the manner in which separation and 
Peanchitiqutake plaee, the bounds on y given by (V.52) 
will always be tight as compared to the bounds implied by 
(V.47). In order to regain the “upper bounded only" form 
for (V.46) with (V.52) appended, we make another trans- 
fermation of variables by taking Z= y ~ P to yield 

Az = g + Or 


0£z¿éQ=F-P 


(SEN) 


E Bau zZandzinteger y JEN) 


Z, > 0 ; JEN 


J 2 


Note that in “decoding” from z to y, integer z <=> 


integer y. 





Zale 


The results of Theorems 4 and 5 given below are 
useful in performing the postoptimality analysis on the 
requirements vector. Theorem 4 is an extension of a 
result due to Nauss [49] for the PMILP (V.27). 
Definition: A PMILP of the form (V.27) iS monotone if 


F(S!) DF(S*) >... > F(sk*) 


, 


where F(s%) denotes the feasible region implied by the 
constraint set oo 
Theorem 4: Given the PMILP (V.27) where 
by =b + 6r, 9€[0,1], and if 
02 8, < 8, & l, then the PMILP is 
monotone and, further, the optimal 
solution value z(8) is a non- 
increasing function of 8. 
Proof: Monotonicity is assured by the choice of 
EV temetcld more restricted problems as 
8 increases, thus giving 
ey > ris 2). 
The non-increasing nature of z(8) is immediately 
obvious as a result of the monotone property. QED 
Note that no stronger result for z(8) is forthcoming as 
in the LP case, since discontinuities will exist, in 


general, in z(@) for certain values of 8 corresponding 


to discrete changes in the integer-restricted variables. 
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The result of Theorem 5 below provides a method 
for implicitly solving problems of the PMILP when only 
One Component OER iS parameterized. That is, in Model I 


(V.2) we consider the special case 


0 
where r, appears in the row corresponding to the element 
of the requirements vector to be parameterized. 
Theorem 5: Let Bix be optimal for K = K}. Also, 


let 
226%, = K} £K.. 
e P3x 2 1 


Mien the solution of which f $, are 
components is optimal for all K such 
that 


K, £K LK 


2 Ih 


Proof: Because K < K, 1 any solution feasible 
for K is also feasible for K, - Hence 
K 
optimal for F(S Ly > optimal for e. 
But since f ¥ Ts such chat 
jk 


2. = 
22 rn 





ele 


1S Optimal for Ki it is true that 
K K 
optimal for F(S Ls £ optimal for F(S zs < 


optimal for As a 


Hence, 
a K 
optimal for F(S `) = optimal for F(S ) 
fo all K, £K £ K}. QED 


As was emphasized in our discussion of branch 
and bound, two essential features are the effectiveness 
of the upper and lower bounds on z in reducing the scope 
of the enumeration and the efficiency in solving the LP 
relaxations. The major thrust of this work has been 
towards solving the LP relaxations of the several prob- 
lems under consideration by postoptimality analysis. The 
Basic Algorithm of Section V.4 is the cornerstone for 
Algorithm One presented below and Algorithm Two to be 
given in the next section. Thus, we state only the modi- 
fications to the Basic Algorithm, which occur in step 4. 
We assume that only a finite number of discrete values 
of 6 € [0,1] are of interest. 
Algorithm One (Postoptimality Analysis of the Require- 

ments Vector) 
Step 4': (LP Solution) At node j, obtain the solution 
to the initial problem under consideration 


(indexed by p) by the UBDS algorithm. If Le 
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feasible solution to (RP 2°, 

set Z5 =-00 ,Wk >p, fathom the node, and go 

to step 2. Otherwise, set p= p+ 1. 

a) Determine So from Rules I-III dE 
O go to (b). Otherwise, update Xp 
using (V.50) and check for possible appli- 
cation of the result of Theorem 5, index- 


ing p as needed. For those problems k 


for which an optimal solution has been 


found, set Z5 = En If all problems have 
been considered, go to step 5. Otherwise, 


i? > 8. ECO” to (b). 
b) Execute UBDS pivots to restore primal 
feasibility. If this cannot be done, set 
ze ==0 ,Wk > p, fathom the node, go to 
step 2. If feasibility is restored, go 
Eo (a). 
Algorithm One was coded by the author in FORTRAN 
IV, incorporating the features discussed in this section. 
Problems of the general nature of Model I (V.2) were 
Duedemly generated (e.g., A, b, C, L, and U); with the 
exception that A did not have the special structure or 
sparsity as in (V.2). Table V.1 indicates some sample 


results. Note that we have indicated the number of con- 


straints (m), total number of variables (n), the matrix 
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density, number of integer-restricted variables, per- 
centage of variables with large upper bounds, and the 
number of PMILP problems to be solved. All of these 
obviously affect the performance of the algorithm. Also 
reported are CPU time IBM(7094), number of nodes in the 
enumeration tree, the number of LPs solved (including 
solution from scratch and pivoting to restore feasibili- 
ty), and the number of times the postoptimality analysis 


routine was called into use. 


V.7 Postoptimality Analysis of the Upper and Lower 
Bounds 
Turning to the postoptimality analysis of upper 
and lower bounds, we state the following general problem 
Form: 
MAX CX 


s.t. Ax = b 


ne Cage UPA) (V.54) 
a > 0 and integer , JEN, 
Xy > 0 ; JEN, 


Eeer ore, we wish to transform (V.54) into an "upper 
bounded only" form in the variables y. After doing so, 


we have the following system of constraints: 
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Ay = g +9 dy = 9d, me 55) 
a O and integer , JEN- 
Y- a 0 , JEN, ’ 
where 
a dJi 
ee ae = AE (Vv. 3 60 
(U-L) lo 
di 
0 dai 
d, = x e a a (57) 
2 
Me 
U*-U dy 
and 
A a; 
d Swe ee ee ne = o (Yag) 
K 
Bi ds 


mie Upper and lower partitions of g, da and sd CONS use 


1° 2 
of m and n rows, respectively. Note that this is the 
same representation of the parameterization as given in 
Model II (V.3). 

Recall that in Section V.2 we discussed the prob- 
lems inherent in simultaneous parameterization for two 
parameters. These problems are certainly of concern 


here. In addition, however, we must consider the prob- 


lem of transforming a solution in y to one in x for the 
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integer-restricted variables. We cannot, in general 
assume that the lower bound on such variables will be 
integral, due to the parametric nature of the bound with 
respect to 9. Further, as will be seen below, the con- 
straints added to (V.55) due to branching will also be 
transformed by lower bounds, but must be integer, Thus, 
PeelSeeirst convenient to fix @ a priori at 6 = 9. Note 
that such an approach requires the solution of as many 
PMILP as there are discrete values ot 90. Next, for 
JEN}, Mimoug@ernetatainiceger X<s> integer y, for those 
values of 8 and j such that us + CU is not integer. 
the lower bound constraint, 

IL, + O(LE-L.)] + lex, , Jen, 
is substituted. Because of this, the right-hand side of 
(V.55) is no longer a linear function of 8. We instead 
may write (V.55) as 

Ay = g +tAd, - w (8) 


> 0 and integer, JEN (NED 9) 


J l 
> 0 , jeN, .** 


IZ 


IS 


J 
Proceeding as before, using the form (V.42) of 


oe, the feasibility condition for the expanded system 


of constraints, 


yn = By lg +A d, - w(0))>0 (V.60) 


yields the following four relationships. Note that 


MINE vector w(8) is formed from d, by appropriately 
modifying those elements of d> ae d by the 
Substitution of lower bounds. 
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9, ar and w, (8) have been partitioned to conform with 


the partition of Bos as 


i 21 1 
22 es u) 
2 
ga - g Ze 22 _ 2 
3 
g 23 3 
i Sa We (8)] 
=] 9421. 36 . -1,,23, ,-1 al 3 
l. (B g) B  Rg3) (AB “Rd, + B "w,(0)- B “R w, (0))>0 
N. 1273 Ze] al 3 21 
2. (g,-B “g,+B “Rg,)+(4B “Rd, +B w.(8) -B R w5 (8)+& d} 
- w> (8) >0 (V.61) 
Pec a? ee) > 0 
> 5 2% a 
4. 9) + as = w (0) 20 


Based on the relationships (V.61), we may define 


the following rules for determining critical values of. 


Rule I: 
-] -] 3 -1 -]. 3 
Bi: .- (B g91),7(B Rg5) ,+(B w,(8)) ,-(B Rw, (8)) = 
oo eo 
(B Rd ) i 
10% (Ei, > 0 
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Rule II: 
1 -1 -=1_ 3 -1 -1_ 3 1 
g3 (B g,),+(B "Rg,),+(B "w, (8)),-(B a a 
A “=MIN : DE 
E a s 
NS e O 
for (B Rd, +d] es 0 
Rule III: 
2 
o` il it 22 
-d 1 
T~ 
i 
Rule IV: 
3 
g - w3 (8) 
4 2» 2 
A” = MIN = = for d= 20 
23 ar 
-dq T 
1 
ai 
Rule V: 


l 2 3 4 
A = MIN er Day Ar oe) 

The critical values | obtained in Rule V is that 
value of A, given $ = 8, such that the current optimal 
basis will become infeasible, necessitating a change in 
the basic variables. Again, the critical values may be 
used for comparison with predetermined discrete X or to 
find all feasible bases of the problem. By an argument 
similar to that given in connection with parameterization 
of the requirements vector, the process is finite. 

when imbedding the postoptimality analysis tech- 


nique for upper and lower bounds into the Basic Algorithm, 
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the process is not as straightforward as in the case of 
the requirements vector. This is primarily because of 
the intimate connection between the process of separa- 
tion and branching and the parameterization involved. 
That is, the set of branching constraints may dominate 
or be dominated by the parameterized bound restrictions, 
the changeover often occurring within the scope of the 
parameterization at the node. In order to examine this 
phenomenon more closely, we begin by considering the 
system (V.59) which, in this case, has already had the 


lower bound transformation applied. We thus have 


Ay = g ta d} - w(@) (V.62) 
A e E 

A > 0 and integer , JEN, 

Y; > 0 ; JEN, 


There are two central issues to be examined with 
respect to the parametric bounding constraints and addi- 
tional branching constraints: 

l. Consistency between the parametric upper 
bounding constraints and lower bound con- 
straints due to branching. (If they are in- 
consistent, this implies infeasibility.) 

2. Dominance of either the parametric upper 


bound constraints or the upper bound 
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constraints due to branching, the latter 
implying non-parametricity of the upper 
bound over the range of for which the 
dominance holds. 

Without loss of generality, we chose to examine the 


applicable portion of constraint set (V.62) 


0 £ y 4g, thd, ~- w, (8) (V. 63) 
ne > 

Y- > 0 and integer, JEN, 

Yj > 0 r JEN, 


With O fixed at 6, we may derive the following limits 
one, 
l. Feasibility/Consistency: 


we 


2 
> Fd] - wo (8) 


must hold for consistency. In general, a‘ Z 0 since we 
consider more restricted problems as e increases, so 
that the value of A beyond which the problem would be 


infeasible is given by 


P. - g, + w, (86) 
j 2 Dias 2 
oa = MIN A lS ford ~~ OC (V.64) 
1 EN a? L: 
JEN, da J 


For the dominance question, when the parametric 
upper bound is applicable we say the bound is "released." 


When the fixed upper bound due to branching is tight, we Say 





227 


the bound is "pegged." 


II. Pegging/Releasing of Integer-Restricted Variables 
as a Function of%. 
The branching constraints will be tight if 
12) A 


1 
Thus, the upper bounds on all integer-restricted vari- 


ables will be pegged for those a 2X51 where 


FE em, 10 5 (v.65) 
ol, = MIN — ~] , ford; #0 
jEN d j 
27 ie 
J 
Further, fore 7 oe A,» each integer variable 


will have its upper bounds again made parametric for all 


A EI where T- is given by 


D. - q). + w, (0) 


: j j | 2 
y. = == AA , for jeN, and d #0 
J 2 1 ily. 
dy. J 
J 
We call mey "release points." The process defined by 


I and II above may be visualized as in Figure V.l, where 
(Dis the range of A for which all the upper bounds on 
integer-restricted variables are pegged, and (2) is the 
range of A for which the problem is infeasible. We 
assume the de have been rank ordered. Note that they are 
not necessarily distinct and not necessarily positive. 


An additional result is possible for the ILP, namely an 
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Figure V.1 


Tllustration of Consistency/Feasibility and 
Pegging/Releasing Rules 
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mohe Solution routine. That is, for all & £a, the 
problem is no longer parametric on& , so it is sufficient 
to solve the problem at «A =o, and by implication 
determine solutions for o a We demonstrate two 
additional results as Lemmas 1 and 2 for the rules pro- 
posed above. 

Lemma 1: For the feasibility and pegging scheme 


described above, 


Proof: Recall that, 1S such that 
_ 2 
= = 9). HALA w, (8), for some JEN, 
J J J 
anda, is such that 
Des + A ES - w, (8) for some JEN 
3%. De ER a 
3) J J 
Thus, PE E and a‘ 2 0 yield the result. 
J 
QED 
Lemma 2: If ds = 0, then E =oo and the upper 
J 
bound due to the branching restrictions is 
tight Wa . 
Proof: Immediate, since the assumption im- 


plies non-parametricity. 
Finally, we would like the expanded system (V.62) 
to be in "upper bounded only" form to facilitate solu- 
tion by the UBDS algorithm. Therefore, we apply the 


transformation of variables 
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2 E (V.66) 
to yield 
AZ =g tad, - wie) - AP (V.67) 
N Sn 
z > 0 and integer , JEN, 
2 > 0 A JEN, 


Note, however, that if we were to examine the con- 


ants implied by the lower partition of g, d and 


Je 


w(8) we would have 


-P £2 £95 tad, - wy le) (V.68) 


1 
But since P is necessarily non-negative, being a vector 
of lower bounds on necessarily non-negative variables, 
the left-hand inequality of (V.68) is dominated by 
a> 0, and hence, (V.67) is the correct general form. 
The objective functions corresponding to the constraint 
sets (V.62) and (V.67) are altered by constant terms 
which do not affect the optimization. 

As was the case in Section V.6, the developments 
of this section involve the modification of step 4 of 
the Basic Algorithm. Hence, we present the new step below, 
the remainder of Algorithm Two being the same as the 
Basic Algorithm. It is assumed that only a finite num- 
ber of discrete values of X€ [0,1] are of interest, these 


being indexed by 1,2,...,K*. 
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Algorithm Two: (Postoptimality Analysis of the Upper 


Step 4': 


and Lower Bounds) 

(LP Solution) At node j, obtain the solution 
to the initial problem under consideraticn 
(indexed by p, o and 8 are operable) 
from the UBDS algorithm. rf 7 a feasible solu- 
tion for (RP DD, set A =-0,WYk >p, fathom the 
node, and go to step 2. Otherwise, p=p+ 1 
and proceed as follows: 
ans ton ene problem 194: to One ii with 

SporopRlate transformations in P and F. Take 
care to insure integrality of lower bounds 
for integer-restricted variables. 


b) Compute & Ay and Yj’ JEN 


JS p 
c) If p = K* fathom the node and go to step 5. 
O set 2. =-@ ,Wk >p, fathom the 


node, and go to step 2. 


d) Otherwise update x 


B using the results of 


(Vaeli modified by the transform 2 = ye- P , 
choosing the proper form for the problem 
upper bound in accordance with Y}, JEN, (e.g., 


i 2 

either 9» +A dy -w, . (9) er Dj=P3" the latter case 
j mo 3 

requiring redefinition of g,, 8}, 


and w. (9). 
e) After updating, if primal feasibility exists, 


set 2. = ze, Be spe o to (C). Otherwise, 
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execute UBDS pivots to regain feasibility. 

If possible, set p > p+ 1l, de = a and go 

to (c). If not possible, set 2. =-0, Yk >p, 

fathom the node, and go to step 2. 
Occ desired”, a transformation frem z to X may be 
made in (e).) 

Algorithm Two was coded in FORTRAN IV and run for 

several randomly generated problems. In addition to the 
ether randomly generated model components, the vectors L* 


and U* were also generated. The*results of the computa- 


tional study are presented in Table V.2. 


Veo Conclusion 

In this chapter, we have developed and proposed 
methodologies for the efficient solution of the class of 
parametric mixed-integer linear programs typified by 
Model I and Model II. From a computational point of 
view, further improvements in algorithmic efficiency based 
on the special structures of (V.2) and (V.3) will be 
discussed in Chapter VI. 

One of the more interesting issues raised in this 
chapter concerns PMILP which are simultaneously parametric 
On two parameters. In Section V.7, we indicated how the 
bivariate case could be reduced to a univariate case by 


Bring Oma Priori and solving the PMILP for the values of 
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o& under consideration. When all values of ¢ have been 
considered, 8 may then be fixed at another value of 
interest, the algorithm reapplied, and so on. This 
essentially requires solving one PMILP problem for each 
value of 8 of interest. 

Recall, however, that in addition to fixing 0 
to reduce the problem to the univariate case, further 
modifications to the lower bounds on integer-restricted 
variables for a given ® were necessary so that integer 
<-> entegers=y. 9lt ts proposed that if there existed a 
dual LP algorithm that did not require a transformation 
of variables to remove non-zero lower bounds, but instead 
handled the lower bounds by a scheme analogous to that 
for the upper bounds, the following strategies might be 
employed. 

First, at any node we could obtain information 
relative to all LP relaxations implied by the combina- 
tions of & and 8 by fixing € at the first value of 
interest and solving the LP relaxations for all%. The 
parameter 8 could then be fixed at its next value, LP 
relaxations obtained for alle , and so on. This process 
may be visualized on the unit rectangle as proceeding on 
the line fromedA =0 to«xA =1 for each fixed 8, and is called 
the "line approach." Figure V.2 illustrates the pro- 


cess. Note that in this case, an additional set of 
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conditions for checking the consistency of the para- 
metric lower bounds with the upper bounds due to branch- 
ing and a pegging/releasing rule for 98 would be useful. 
Care would have to be taken, however, that the feasi- 
bility conditions derived for A and 9 were applied for 
the proper combinations of values of these parameters. 
Although the line approach would give a "one-pass" 
algorithm as opposed to the repeated application of 
Algorithm Two, it is not intuitively clear that the gain 
in computational efficiency would be substantial. This 
is due mainly to the relative appropriateness for the 
several combinations of “ and 8 of the branchings under- 
taken, 

A second approach that might be useful at a 
node fixes 8 at its initial value and solves the LP 
relaxations for all values of A of interest. The para- 
meter A is then fixed at its maximum value and LP relaxa- 
tions obtained for all 9 of interest. Then, fixing 8 
at its final value not previously encountered, consider- 
ation could be shifted toe, and so on. Note that we 
would never have to reconsider combinations of A and ® 
already examined. For this reason, this process is 
called the "Spiral approach," and is illustrated in 


Figure V.3. The principal merit of the spiral approach 
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is that when we begin parameterization with that para- 
meter previously held constant, an optimal basis is im- 
mediately available. This is not the case in the line 
approach, where the problem with a new value of 8 must 
be started from scratch. As was the case in the line 
approach, conditions on consistency and pegging/ 
releasing would have to be made available for both para- 
meters, again with the caveat on their proper applica- 
tion. In addition, a set of postoptimization rules 
analogous to those in Section 7 for & would be necessary 
for 8. The additional bookkeeping required in the spiral 
approach would appear to be a significant consideration. 
The development of the upper and lower bounded 
dual simplex algorithm and its incorporation into PMILP 
algorithms based on the line and spiral approaches 
appears to be an area for further research. Comparison 
of computational experience from such methods with that 
obtained from our proposed methodology would be useful. 
Finally, a direct attempt at using Algorithm Two in a 
line approach might prove to be an interesting applica- 


cion- 
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Line Approach to Reduction of Bivariate Parameterization 
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Spiral Approach to Reduction of Bivariate Parameterization 
(Illustrated for K* = 5) 
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Chapter VI. An Application of the Methodology 


VI.l Introduction 

Having developed a patient classification 
system for level of care required, the Basic Staffing 
Model, and extensions to the Basic Staffing Model, 
we turn to an example of the way in which the metho- 
dology can be applied. We shall demonstrate the efficacy 
of Models I and II for use in determining staffing 
patterns, the computational characteristics of algorithms 
for each model having been previously shown. Asa 
basis for the application, the nurse staffing problem 
and the underlying assumptions will be delineated, 
along with the necessary data elements. Some computa- 
tional considerations for problems of the form of 
Models I and II will also be discussed. Foliowing 
this, applications of Model I will be presented, with 
consideration devoted to variation in the total 
service level and changes in personnel budget restric- 
tions. Model II will be used to demonstrate the 
effects of alterations in the upper and lower bounds of 
care activity times on nursing service requirements. 
Finally, comments are offered on policy implications 


for institutional care. 
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VI.2 Problem Construction and Computational Considerations 

As the setting for our example, we have selected 
a representative proprietary long-term care facility 
in the Baltimore metropolitan area. This facility 
maintains 100 licensed beds for both Intermediate A 
and Skilled Nursing Care patients, and is generally 
recognized by professionals as providing good quality 
care. A facility with two levels of care, such as 
this one, is typical; only the much larger and more 
complex facilities have three or perhaps four levels, 
including chronic care. Note that the consideration 
of two levels of care instead of the three basic 
levels for which our models have been designed can be 
accomplished quite easily; such a problem is just a 
subset of the larger problem. Initially, we assume a 
patient mix of 50 Intermediate A care patients and 
50 Skilled Nursing Care patient as determined by the 
classification system of Chapter IV. 

Recall that in Chapter III certain data necessary 
for the specification of the staffing models were 
omitted. These include: 

1. Upper and lower bound restrictions on 
nursing staff by skill level; 

2. Salary costs by skill level; 

3. Personnel budget allowance. 


Since we have previously indicated that the methodology 


is best applied to the daily (24-hour) staffing problem, 
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items (1), (2), and (3) will be presented within this 
context. 

In setting the bounding restrictions on nursing 
staff, there are two major considerations. First, 
sufficient staff must be on duty each day to satisfy 
state regulations. Obviously, this constraint may be 
incorporated into the lower bounds. On the other hand, 
the exigencies of the local labor market for nursing 
personnel, especially licensed (R.N. and L.P.N.) nurses, 
set an upper bound on the number of nurses of each 
skill level that can be hired. Of course, budgetary 
restrictions may enforce a staff mix that is consi- 
derably below the upper bounds set in this manner. 

In order to set lower bounds, nursing home regulations 
for the State of Maryland were consulted.* These 
regulations specify, as a minimum, that for both 
Intermediate A and Skilled Care patients the ability 
to provide two hours of bedside care per patient day 
must be maintained. In addition, a staff-to-patient 
ratio of 1:25, again as a minimum, must be adhered to. 
Finally, for the assumed mix of Intermediate A and 
Skilled Care patients, two R.N.'s and three L.P.N.'s 
must be on duty in a supervisory/charge and care- 


providing capacity. Additional regulations concerning 


* Maryland State Department of Health Regulations 43602, 
43G06, and 43G07. 





241 


licensed nurse coverage by shift are provided, but 
are not essential to the models. Taking all of these 
factors into consideration, upper and lower bounds on 
nursing staff were set as shown in Table VI.l. It 
should be noted that the minimum mix represents a 


least-cost feasible solution. 


Table VI.l 
Upper and Lower Bounds on Nursing Personnel 


for the Example 








R.N. L.P.N. N.A. 
Upper Bound 8 20 40 
Lower Bound 2 3 23 


In computing the values in Table VI.l a 7 1/2-hour 
shift was used as the basis for determining compliance 
with minimum daily bedside care requirements. 

With respect to daily staffing costs, data 
were obtained from the facility under study as to costs 
for licensed and non-licensed nurses averaged over 
all three shifts, including fringe benefits. After 
application of an adjustment for the salary differential 
between R.N.'s and L.P.N.'s, the average daily staffing 
salary costs for each skill level were determined to 


be as shown in Table VI.2. 
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Table VI.2 
Average Daily Salary Cost of Nursing Personnel, 
by Skill Level, for the Example 


R.N. EmP „N. N.A. 





546.88 29215 2072 125 


The current average daily nursing personnel budget 
for the facility was determined to be approximately 
25.00. 

Because of the use of a 7 1/2-hour shift in 
our example facility, the results of Table III.8 
concerning nursing time availability for patient- 
centered activities per shift must be modified. The 
7 1/2-hour shift is derived by subtracting one-half 
hour for meals from the 8-hour shift. This being the 
case, the percentage allowance for meals was removed 
from the 8-hour shift data to give the percentages 
of time, per shift, devoted to non-patient-centered 
activities. By applying these percentages to the 
7 1/2-hour shift, the availabilities in minutes as 


shown in Table VI.3 were obtained. 
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Table VI.3 
Recomputed Availabilities, per 7 1/2-Hour Shift, 


in Minutes by Nursing Skill Level 


Availability 
Skill Level Percentage Per 7 1/2-Hours 
R.N. 82.4 370.8 
Loo N. 84.6 380.7 
N.A. 69.1 311.0 


With the remainder of the data for the problem 
as described in Chapter III, the models are thus 
specified. Before beginning the application, however, 
certain computational features of Models I and II are 
worthy of mention. Recall that in the description of 
the Upper Bounded Dual Simplex algorithm in Chapter V, 
no specific computational strategies were specified. 
The implication was that the algorithm would be 
executed on a complete simplex-type tableau. Examina- 
tion of the constraint sets of Models I and II reveals, 
however, that the corresponding constraint matrix 
contains quite a large number of zeroes; that is, 
it is sparse. Significant improvements in algorithmic 
efficiency are possible in such cases, due mainly to 
decreased computer storage requirements and numbers 
of calculations. By storing non-zero constraint 


entries in lists and using an upper-bounded dual form 
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of the revised simplex method (cf. Hadley [28]), the 
requirement for maintenance of a complete tableau is 
eliminated, thus improving efficiency. These aspects 
were incorporated into Algorithms One and Two when 


solving the example problems of this chapter. 


VI.3 Example of Model I for Service Level Variation 

We turn now to the application of Model I 
(III.2) to our example problem. In this section we 
will consider the effect on the objective function 
and staff mix of variations in the total number of 
minutes of patient-centered activities provided, S. 
Such a study assists the administrator in determining 
the sensitivity of the proposed model to varying 
amounts of service provided. By fixing the personnel 
budget at any given level, a determination can be made 
as to how much service such a budget would support, 
and where possible improvements could be made vis-a-vis 
staff mix and assignment. 

For our purposes, we have chosen to set the 
daily budget at a high level ($2000.00) to yield an 
indication of the ideal staff mix to provide any given 
level of service. Three proposed patient mixes, by 
classification, are examined: 50 Int. A and 50 Skilled, 
40 Int. A and 60 Skilled, and 60 Int. A and 40 Skilled 


patients. The range of S for consideration was chosen 
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in each case to approximate the difference between 
the maximum and minimum service levels obtained by 
summing the upper and lower bounds, respectively, 

for each patient mix. For the sake of brevity, we 
present results in tabular form, omitting detailed 
specification of the assignment patterns. Tables VI.4 

through VI.6 display the information, where Z indicates 


the objective function and S is constrained to be 


less than or equal to the value indicated. 


Table VI.4 


Model I Solution for Various Service Levels for a 
Patient Mix of 50 Int. A and 50 Skilled Patients 


Service Level(S) Staff Mix* S to Support 
(mins.) Z RN,LPN,NA Staff Mix 
16500 32419.4 8 19 20 1673.29 
15720 32419.4 8 19 20 LE 73:29 
14940 31395.4 8 17 19 1566.64 
14160 29835.4 8 15 18 1459.99 
13380 28275.4 8 14 18 1420.24 
12600 26384.4 8 13 16 1326719 


* Actual staffing; the lower bound of 23 for NA given 
in Table VI.1 may be relaxed due to sufficient staff 
to meet regulations. 
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ao heey ee 


Model I Solution for Various Service Levels for a 
Patient Mix of 40 Int. A and 60 Skilled Patients 


Service Level (S) Siendo ee Mees STO MS UPport 
(mins.) Z RN, LPN, NA Staff Mix 
17400 32174.4 8 20 20 171202 
16600 32174.4 8 20 20 1713.02 
15800 31574.6 8 20 18 1658.74 
15000 29974 .6 8 18 18 1579724 
14200 28374-6 8 16 17 1472.59 
13400 26483.6 Suelo iy 1432.84 


*Actual staffing; the lower bound of 23 for NA given 
in Table VI.1 may be relaxed due to sufficient staff 
to meet regulations. 


Table VI.6 


Model I Solution for Various Service Levels for a 
Patient Mix of 60 Int. A and 40 Skilled Patients 


Service Level (S) Staff- MTR" $ to Support 
(mins. ) Z RN, LPN, NA Staff Mix 
15500 31734.4 8 18 16 1524.94 
14560 31370.4 8 5 220 1514.29 
13620 237936. 2 8m 13719 1407.64 
12680 217916.2 8 LO 1328.14 


11740 and below--no feasible solution 
*Actual staffing; the lower bound of 23 for NA given 


in Table VI.1 may be relaxed due to sufficient staff 
to meet regulations. 


As might be anticipated, the cost ot providing 
care for a patient mix requiring predominantly Skilled 
care is uniformly greater than either of the costs 
there 


for other patient mixes examined. In addition, 
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is greater reliance on licensed nurses for care provi- 
sion in this case than for either of the other mixes. 
The normative nature of the model and data, as well 

as the relaxed budget constraint, are primarily 
responsible for the heavy reliance on R.N.'s in all 
three configurations. Two final points are of interest. 
In all three cases, the largest increase in the objec- 
tive function per dollar expended occurs at the lower 
end of the spectrum. At the upper end of the service 
level range for the 50-50 and 40-60 patient mixes, 
further increases in Z are not forthcoming beyond the 
S values of 15720 minutes and 16600 minutes, respec- 
tively. The indication here is that the true upper 
bound on S lies somewhere between the two largest 
values and/or the model is being forced into making 
assignments for which the appropriateness scores are 


zero. 


VI.4 Example of Model I for Budget Variation 

While the examination of changes in the optimal 
solution to the BSM as a function of the service level 
via Model I is, perhaps, more of academic interest, 
the use of Model I for testing the effects of altera- 
tions in the personnel budget can be of great value 
to an administrator. In this regard, Model I was 


applied for the case of a patient mix of 50 Int. A 
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and 50 Skilled care patients. The budget was chosen 
to range from $1200 to $1800 per day, in steps of 
$200. This range would include feasible staff mixes 
for both the maximum and minimum total service that 
can be provided under the upper and lower bounds in 
use. 

Results of this application are shown in detail 
in Figure VI.1. The output not only shows the objec- 
tive function, staff mix, and budget constraint, but 
also the optimal assignment pattern by care area, 
nursing skill level, and patient classification 
(A=Int. A and S=Skilled care). These results are 


Summarized. in Table VI.7 below. 


Table VI.7 


Model I Solution for Various Budget Levels for a 
Patient Mix of 50 Int. A and 50 Skilled Patients 


Total Service Staff Mix $ to Support 
Budget ($) Provided (min.) Z RN, LPN, NA Staff Mix 
1800 16194 32419.4 8 19 24 7781029 
1600 15829.9 31156.9 IS 1595.74 
1400 15926, 23 28027.4 8 10 23 1396.99 
1200 12463 21404.6 2 12 23 PU L 


Several comments on this result are in order. 
First, the increase in the budget from $1200 to $1400 
is obviously of high significance, as it precipatates 


erein the stati mix that allows a large increase 
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Figure VI.l 


Results of Application of Model I 
for Budget Variation 
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Frgure VI.l (cont'd.) 
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Figure VI.l (cont'd.) 
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in the objective function. The increase appears to 
yield not only more appropriate assignments, but 

also to increase the time allocated to a number of 
high priority care areas (viz.,4,7,10,22, and 23). 

It is also of note that as additional funding becomes 
available the majority of additional nursing time 

is allocated to Skilled patients. The policy impli- 
cations of this phenomenon will be discussed in the 
next section. Finally, examination of the allocation 
of nursing time between high and low priority care 
area/classification combinations* reveals that under 
the exigencies of decreased funding, allocation to 
high priority combinations decreases at a slower rate 
than for low priority combinations. As an illustra- 
tion, consider care area 1 for Skilled patients. The 
allocation remains at the upper bound (475 minutes) 
until funding is constrained to $1200, at which time 
it decreases to 450 minutes. The allocation to care 
area 2 for Int. A patients, on the other hand, begins 


to decrease at a budget level of $1600. 
VI.5 Example of Model II 
As the final stage in our example, we present 


an application of Model II (111.3). Having determined 


* See Table III.9 
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the effects of changes in the service level and 
budget constraints under the upper and lower bounds 
postulated in Chapter III, we would like to determine 
the effect of alterations in these bounds on the 
staff mix, nursing time allocations, assignment 
pattern, and necessary budget. In order to place this 
study in a realistic context, we will analyze the 
possibility of decreasing the normative bounds by 
various percentages in order to obtain a feasible 
solution at the current facility daily nursing 
personnel budget. 

The results of the application of Model II 
are shown in Figure VI.2. The maximum budget was 
set at $900. Theta at a value of zero indicates a 
30% decrease for each care area/classification 
lower bound, while theta at one indicates a 15% 
decrease. The parameter alpha at zero yields the 
original upper bounds, while alpha equals one implies 
a 15% decrease in all upper bounds. As can be observed 
in Figure VI.2., theta equals zero and alpha equals 
one yield a staff mix that may be maintained under 
the current facility budget. This is at the expense 
of a decrease in total service provided of 3055.5 
minutes and a decrease in the objective function of 


4128.7 units. In terms of increased budgetary 
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Figure VIv2 
the Upper and Lower Bounds 


Results of Application of Model II to 
in 


Variation 
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allocations, the proportional increase in the objec- 
tive function is not as great in moving from a budget 
of $900 to $1200 as that for an increase from $1200 

to $1400; on the other hand, the proportional increase 
in total service provided is greater. This would 

seem to indicate that relatively inappropriate 
assignments outweigh the increased allocation of time 
to high priority combinations, and it is not until 
further increases in the budget are made that more 
appropriate assignments are possible. 

An interesting feature of this analysis concerns 
reimbursement rate policies for Medicare/Medicaid 
patients. At current rates, the facility may provide 
a total of 9407.5 minutes of patient-centered care to 
its patient mix, allocated as 3615.8 minutes to Int. 

A and 5791.7 minutes to Skilled care patients. Under 

an increase in personnel budget of $310.87 per day 

a total service level of 12463 minutes can be maintained, 
allocated as 4472.5 minutes to Int. A and 7990.5 

minutes to Skilled care patients. The most significant 
increase in care provided is obviously for Skilled 
patients. For policy purposes, then, what would be 

the proportional increase in daily reimbursement 

rates for each level of care if it were decided to 
pursue the higher budget policy? Using the assumed 


patient mix, the increase should be apportioned as 
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$1.74 per day for each Int. A patient and $4.48 per 
day for each Skilled care patient. Carrying the 
analysis one step further, the policy of increasing 
personnel budget from $1200 to $1400 would require an 
additional assessment of $2.62 per day for each Int. 

A patient and $1.41 per day for each Skilled care 
patient. This type of resulto can be of great value to 
administrators and regulatory agencies alike, because 
it places in clear perspective the costs of providing 


additional care on a per capita basis. 


VI.6 Conclusion 

The sample application presented in this 
chapter gives an indication of the type of information 
made available to the facility administrator and 
director of nursing through use of the models and 
algorithms developed in this dissertation. Due to 
the exigencies of computer funding, the application 
may not have been as extensive as might have been 
desired; yet, with the exception of showing an example 
of Model II as applied to changes in the bounds on 
personnel, the range of possible model employment 
has been indicated. Use of the models to assess 
alternative courses of action with respect to increased 


service provided and budget has been stressed, but 
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it should be clear that the alternative application 
to service and budget cuts is immediate. In addition, 
by setting the budget constraint at various levels, 
the service level model can assist the administrator 
in recognizing suboptimal staffing patterns. Finally, 
further exploration into the policy aspects of the 
models could prove to be an interesting application 


of the methodology. 





Chapter VII. Summary, Recommendations, and Extensions 


Vill Se introduction 

This dissertation represents an effort to 
respond to the changing environment of and demand for 
long-term care. The approach taken has been essentially 
structure-oriented at the facility level, and is 
specifically aimed at providing managerial decision- 
makers with information that is basic for planning a 
nurse staffing program that responds more effectively 
to patient needs. Such information is obtained from 
the solution of two parametric representations of a 
normative staffing model. The model itself is based 
upon a patient classification system that has been 
derived from a comprehensive patient assessment instru- 
ment, and has been shown to bear a direct relationship 
to demand for nursing services. Moreover, as a direct 
result of the above efforts, certain theoretical advances 
have been made in the area of parametric mixed-integer 
linear programming. 

In this chapter we will briefly summarize the 
findings of our study, provide additional insights into 
areas of concern, suggest procedures for model implemen- 
tation, and recommend areas for further research. 
Finally, some possible extensions of the current effort 


will be presented. 
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VII.2 Previous Work 

Several studies concerned with patient classi- 
fication and nurse staffing in both acute and long-term 
care settings were discussed in Chapter II. It was 
observed that the early work of Connor [18] established 
the concept of patient classification based on the 
degree of patient independence in what are now called 
activities of daily living (ADL). Young's [84] 
regression-based pediatric patient classification 
system demonstrated the feasibility of using a statis- 
tical methodology to derive meaningful prediction of 
demand for nursing services. As was the case in both 
the Connor and Young studies, the classification 
systems developed by Poland, et al. [54] and the Hospital 
Association of New York State [30] for acute care 
patients attempted to link patient dependency status 
to nursing service demands. Typically, however, the 
predictions were given in terms of total nursing care 
time demanded for a patient mix on a nursing unit, with 
no specification as to the partitioning of the demand 
by nursing skill level. 

With respect to long-term care patient classifi- 
cation, it was also found in most cases that the 
relationship between a patient's health status and a 
detailed specification of demand for nursing services 


was not provided. The classification systems of Katz, 
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et al. [36] and Burack [13,50] serve as typical examples 
of this point. Nevertheless, these two systems proved 
to be significant; the former for its prescient 
recognition of the concept of ADL, and the latter for 
its emphasis on goal-oriented care. The RAPIDS system 
of Salmon, et al. [61,62] addressed the issues of 
patient needs and placement vis-a-vis level of care; 
but it should be noted that acceptable results were 
reported only with regard to the placement question. 
The two studies by Parker [52,53] represent an initial 
effort towards the use of various statistical techniques 
in the derivation of long-term care patient classifi- 
cation. Through use of Bayes' Theorem on patient 
assessment data from an abbreviated list of binary 
health status indicators, and through the use of 
cluster and discriminant analysis on data derived from 
assessments using the CPAI, insights were obtained into 
the feasibility of developing effective classification 
systems, and key indicators, as a basis for further 
analysis. Finally, the Colorado study by McKnight 
[41,42] quantified the relationship between a tri-level 
patient classification system and existing patterns 

of nursing service care activities. Data from this 
study have been a key ingredient in our presentation. 


Several proposed nurse staffing models were 
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surveyed, among which the mathematical programming 
approaches of Wolfe [80] and Liebman [38,39] were found 
to be most relevant to this work. Wolfe's integer 
programming model, although developed for acute care 
units, was nevertheless of interest because of its 

use of groupings of nursing tasks into task complexes. 
Although such groupings were suggested in order to 
reduce the model size for computational ease, they may 
have been too broadly drawn to allow for the incorpor- 
ation of specific allocations and assignments of staff 
to individual care activities. Extensions to Wolfe's 
model proposed by Shuman [66] were shown to provide 
additional features that brought the original model 
closer to a representation of reality. The work of 
Liebman was presented as the only known application, 
previous to this dissertation, of mathematical program- 
ming techniques to long-term care personnel allocations. 
In this study, the insights provided as to the nature 
of objective function forms and coefficients proved to 
be invaluable; the use of psychometric measurement 
techniques and determinations of preferred staffing 
patterns (actual vs. predicted) formed a large part of 
the Liebman study and appeared to be most appropriate 
for specifying an objective function. The proposed 


models, however, were structured in such a way as to 
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make the presentation of alternative courses of action 
difficult, with a consequent increase in computational 
effort required in order to develop such information. 
As a final note, the extensive application of 
computer simulation to long-term care nurse staffing, 
assignment, and allocation was surveyed. The studies 
OF Turner, et al. [74], McKnight and Steorts [43], and 
Hundert [31] were shown to provide estimates of staffing 
and assignment patterns based on examination of 
resources predicted to meet a series of randomly gener- 
ated patient demands. Although simulation models have 
become more sophisticated over the years in their 
quantification of systems, their results are subject 
to statistical error. It was also observed that 
inadequate provisions were incorporated in these models 
for a parametric examination of model data and model 


results. 


VII.3 The Proposed Methodology 

In retrospect, it can be seen that few of the 
studies surveyed in the long-term care area are capable 
of relating patient needs, as reflected in health 
status, to the provision and assignment of adequate 
staff to meet those needs. Recall that in Chapter III 
it was suggested that the statements of the administra- 


tor's staffing problem, and the aggregate budgeting 
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problem, given by Baloff, et al. [8] capture the essence 
of the type of information necessary to plan an 
effective staffing program and the reasons why the 
derivation of such information is non-trivial. The 
unified approach presented in this dissertation is 
believed to embody a rational technique for developing 
such information, in a manner useful to facility 
administrators and directors of nursing. 

As the initial step in the derivation, it was 
demonstrated that the care area performance time data 
from McKnight's Colorado study [41,42] could be 
interpreted as performance times, per 24-hour day, for 
"typical" or "average" patients belonging to each of 
three homogeneous classification groups. Further, it 
was remarked that the groupings, designated as the 
minimum, moderate, and maximum care categories of the 
Colorado study, are analogous to the Intermediate B, 
Intermediate A, and Skilled levels of care. From these 
data, the set of normative average time estimates were 
derived. 

The next step was the development of an easily 
applied patient classification system that could group 
patients into categories representing the three levels 
of care indicated above. Major criteria for group 


identification were the functioning ability and 





268 


psychosocial status of the patient, as well as medical 
status as reflected by certain medically-defined condi- 
tions. After extracting a representative sample of 
patient assessments based on the Collaborative Patient 
Assessment Instrument (CPAI) and 37 CPAI variables 
indicating overall patient status as defined above, a 
regression-based method due to Walker and Duncan [78] 
was applied. This method was held to be appropriate 
because of its ability to use an ordered, polychotomous 
response variable; that is, use of the appropriate 

level of care for a patient as adjudged by a professional 
panel. Further information concerning the assumptions 
and derivation of the Walker and Duncan method were 

also discussed. After analysis of the 37-variable 
problem, five subsets of the original variables contain- 
ing 12 variables and one subset containing 9 variables 
were chosen for further study by recourse to results 

of previous Studies and by backward elimination of 
several of the 37 variables. One such 12-variable 
subset was shown to possess good prediction and recog- 
nition powers, and based on comparison of several 
measures of performance with other subsets, was chosen 
as a reasonable answer to the subset selection problem. 
This subset, in turn, became the basis for an implement- 


able patient classification form proposed in Chapter IV. 
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It was noted that this patient classification system 
could serve not only for nursing service demand predic- 
tions, but also as a tool for use by officials 
responsible for legal placement of patients. 

Having thus developed a readily-applied patient 
level of care classification system and having shown 
how nursing service demands for "typical" patients of 
each category could be estimated, it was proposed that 
average unit or facility-wide demand could be predicted 
on a daily basis. This was accomplished by first 
determining the numbers of patients in each category, 
then multiplying each care area/classification average 
demand by the appropriate group membership. 

Determination of an optimal nurse staffing mix 
and allocation and assignment patterns could then be 
addressed. The Basic Staffing Model (BSM) proposed in 
Chapter III incorporates those features recognized as 
being most germane to the administrator's problem. A 
synthesis of both appropriate assignment of nursing 
skill levels to care area/classification level demands 
and assured allocation of nursing resources to high 
priority demands was proposed as the basis for an objec- 
tive function to be maximized. Constraints that provided 
for adherence to legal staffing guidelines, satisfaction 
of bounded demands by care area and patient classification, 


and recognition of the availability of nursing resources 
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were incorporated. The BSM was formulated as a mixed- 
integer linear program, and, based on the assumptions 
and data proposed in Chapter III, provided a 24-hour 
staff mix, allocation, and assignment pattern for a 
ünit or facility. 

Recognizing that presentation of a single 
optimal solution can be of dubious benefit to the 
administrator, an extension to the BSM was derived 
that allows for the analysis of alternative courses 
of action with respect to the total service provided 
and the personnel budget allocated. These aspects were 
incorporated in Model I. As an additional feature, 

a method for inspection of model sensitivity to the 
upper and lower bounds on demands as well as alterations 
to staffing limitations was proposed and synthesized 
into Model II. Models I and II were both cast as 
parametric mixed-integer linear programs, and algorithms 
for their efficient solution were derived in Chapter V. 
Computational results for general problems in tke 
generic classes of both models were also reported. 

As a final step in the presentation, a repre- 
sentative application of the methodology was given in 
Chapter VI. Using data obtained from an existing 
long-term care facility, the response of the objective 
function, staff mix, and allocation and assignment 


patterns to alterations in the total service provided, 
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personnel budget, and upper and lower bounds on demands 
was analyzed. In general, it was evident that the 
model displayed a fairly high degree of sensitivity 

to its parameters, thereby reiterating the necessity 
for analysis of this sort. The information derived in 
the analysis should prove to be valuable for the 
facility administrator and the director of nursing, 
both for planning a staffing program and for indicating 
how increased amounts of care (and their cost) can be 


evaluated. 


VII.4 Recommendations and Extensions 

Before proceeding further, a major point must 
be emphasized. Although the models and methodologies 
set forth in this dissertation are deeply rooted in 
concepts and data obtained from field studies of 
long-term care facilities, some of the proposals 
remain, by nature, Speculative at this point. Therefore, 
before undertaking further theoretical work on patient 
classification systems and nurse staffing models, it 
is recommended that some form of implementation, 
testing, and evaluation of the work proposed herein 
be attempted. In this way, further experience and 
feedback can provide information that would be essential 
for extending and perhaps modifying the procedures 


suggested in this study. 
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Turning to specifics, some factors with respect 
to the CPAI data used in the derivation of the patient 
classification system are worthy of comment. The 
CPAI is currently undergoing revision and field- 
testing, a major consideration being an attempt to 
better quantify the psychological and social functioning 
status of a patient. That such factors are of impor- 
tance in planning for and providing care for long-term 
care patients is well known. A second point is that 
due to a lack of data in the proper form, the special 
care factors listed on the last page of the CPAI* were 
not included in the current study. Intuitively, it 
would appear that both the expanded psychosocial and 
special care factors might affect the results presented 
here; the expectation is that inclusion of such factors 
would serve to increase the predictive capability of 
the model. Therefore, it is recommended that upon 
derivation of a representative sample of patient 
assessments using the revised CPAI, the analysis of 
Chapter IV be reapplied. Comparison of results could 
be highly informative. 

Several recommendations concerning data used in 
the staffing models are in order. No new and extensive 


nursing activity analysis data for long-term care has 


* See Appendix A. 
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been produced since McKnight's Colorado study, although 
such an undertaking is currently underway at The 

Johns Hopkins University. Data obtained with respect 
to amounts of care provided by level of care and by 
existing allocation of nursing time would be invaluable 
as a supplement to those data made available for this 
dissertation. Subject to interpretation of the new 
results in a manner analogous to that proposed here, 
the models, in their current forms, would readily accept 
such information. As a further recommendation, data 
obtained by subjective estimation (e.g., appropriate- 
ness scores, upper and lower bounds on demands, and 
priorities) should be replicated, perhaps by more 
extensive use of psychometric measurement techniques 
such as the Q-sort. 

Recall that in considering the allocation of 
nursing time to those tasks that are non-patient- 
centered, a per-shift allowance (by skill level) was 
subtracted from the time on duty during a normal shift. 
In this way, time available to devote to patient- 
centered activities was obtained. As a possible 
extension to the models, the grouping of non-patient- 
centered tasks into clusters similar to care areas 
might be attempted, with the added feature that 
appropriateness scores and performance times be 


attached to each such cluster. Assuming the availability 


AMT 





274 


of such additional data, the proposed models could 
become more useful to managerial decision-makers. 

A final comment concerns the inclusion of such 
activities as physical therapy, occupational therapy, 
and recreational therapy in the models. What little 
time was observed as being devoted to these activities 
in the nursing homes of the Colorado study was included 
in the performance times used here and allocated to 
nursing staff. Should additional data become available, 
these activities and the specialized personnel charged 
with their direction could easily be included by adding 
appropriate variables and constraints in extended 
versions of our models. 

In retrospect, the unified approach to patient 
classification and nurse staffing presented in this 
dissertation represents a first step towards relating 
patient health status to the managerial decisions 
faced by facility administrators and directors of 
nursing. Innovations in meaningful and readily applied 
classification and staffing models have been stressed, 
as well as demonstration of the methodology. It is 
believed that such an approach will prove itself in 
meeting the changes in the long-term care system almost 


certain to occur in Succeeding years. 
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M-t52 


PATIENT CLASSIFICATION 


FORM 


SOCIO-DEMOGRAPHIC DATA 





(Check Those Which 


MA ce cal A O E A 
g T T es ee A [| 
Month Day GOA 
3. BIRTHPLACE __ USA. SPECIFY STATE 
(Check One) 
— OTHER., SPECIFY COUNTRY a 
u, SEX MATE 
(Check One) |_| 
__ FEMALE 
5, — WHITE __ AMER. INDIAN — CHINESE — OTHER 
(Check Ons) SPECIFY) 
__NEGRO __ JAPANESE —_FILLIPINO 
6. MARITAL STATUS __ SINGLE ___ WIDOWED _ SEPARATED 
(Check One) 
—_MARRIED —~ DIVORCED HOW MANY YEARS? i 
7. FELIGIOUS PREFERENCE __ NONE —— JEWISH — OTHER 
(Check One) (SPECIFY) 
— CATHOLIC __, PROTESTANT 
8. RESIDENCE ADDRESS Weber. "7 ie o Y “tata Poy Coda a ae 
9, PAT {EXT LOCATION RESIDENTIAL HeaLtH Cage Facitity 
N __ PRIVATE RESIDENCE __ DOMICILIARY/PERSONAL CARE __ CHRONIC DIS./REHAB. HOSP, 
INTERMEDIATE CARE _. MENTAL HOSPITAL 
— REYTED ROOM(S — KURSING HOME __ OTHER SFECIALTY HOSP, 
COMMERCIAL __ EXTENDED CARE — SHORT IERM ACUTE HOSP, 
10. LENGTH OF TIME 1E RESIDENTIAL LE HEALIH CARE FACILITY. 
AT LOCATION SPECIFY YRS, MOS. SPECIFY TAIS AUMISSION 
DATE CF: 
li. LIVING ARRANGEMENTS ___ OWN HOM ___ ALONE 
__ ANOTHER S HOME ___WITH SPOUSE 
(Check Those Which — PAYING  __ NOT PAYING WITH OTHER, WHO? 
Apply) —— RENTED ROOM(S) __HE3ZLTH rRZLATED Pr 
SPECIFY TYPE 
u 12. LIVING CHILDREN SPECIFY BY NUMBER: SONS: DAUGHTERS; a 
13. EDUCATION —— NO SCHOOL ING __ TRADE, TECHNICAL. VOCATICNAL SCHOOLING 
. SPECIAL EDUCATION. SPECIFY YEARS _ SOME UNDERGRADUATE COLLEGE. SPECIFY YRS. ___ 
(Cheek Those Which GRADES COMPLETED. SPECIFY — COLLEGE COMPLETED 
Appty) — HIGH SCHOOL DIPLOMA —_ GRADUATE COLLEGE 
14, USUAL OCCUPATICN SPECIFY: OR CHECK ONE: __ MILITARY 
: ___. HOUSEWIFE __ NEVER EMPLOYED 
15. EMPLOYMENT STATUS — EMPLOYED __WORKING, HRS. /WEEK 
Apply) ~ RECEIVING PAY 
— RETIRED 
— UNEMPLOYED m 
16, FAMILY INCOME — LESS THAN $3,900 
ec be 
$7.06 - $9.999 i 
17. NEALTH CARE COVERAGE — NONE 
(Chach Those Which __ SELF-PaY 
Apply) — NO CHARGE 
— MEDICARE 
— MEDICAID | 7 
DIAGNCSIS(ES) ee _ 2 
(List) 


JNU OPERATIONS RESEARCH, 6 HANYS, JUNC, 1972 
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DATE OF CLASSIFICATION 

INTERVIEWER 
RESEARCH GROUP 

PATIENT NUMBER 













INSIDE HOME/FACILITY 
OUTSIDE HGME/FACILITY 


— RECEIVING PENSION 
— NEVER IN LABOR MARKET 


—..$5.000 - $6.939 
—..$]5.000 Prus 


— 53.000 - 54.999 
— $10.090 - $14.599 


LIST ALL OTHER TYPES OF HEALTH INSURABCE PATIENT MAS 


— mann (nn 


lor 4 
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n-130 
INTERVIEWER DATE OF CLASSIFICATION 
RESEARCH GROUP PATIENT RUMBER 


FUNCTIONING STATUS 


, 





ITEMS 
18. MOBILITY LEVEL (cuech eed Complete es Applicable) 


GOES OUTSIOE HOUSE/FACILITY MOVES ADOUT HOUSE/FACILITY CONF ERED 
W/O KFL? __H/0 HELP CHAIR 
—NITM HELP. DESCRIBE __MITW HELP. DESCRIBE — BED 








nee 


DESCRIBE ANY DEVICES USED 























19. TRANSFERRING (coach aed Complete ee applicable) DESCRIBE ANY DEVICES USED 
TRANSFERS DOES NOT TRANSFER 
W/O HELP __DCME BY OTHERS 
MITH HELP, DESCRIBE __BEDFAST 
20. WALKING (Check aed Cosplate an Appliceble) DESCRIBE ANY DEVICES USED 
WALKS DOES NOT WALK 
W/O WELP ___BEDFAST 
WITH HELP. DESCRIBE — _CRAIRFEAST 
21. WHEELING (Chach ced Complate ea Applicehla) DESCRIBE ANY DEVICES USED 
WHEELS DOES NOT WHEEL 
W/O HELP WALKS 
___WITW HELP. DESCRIBE ~_CHAIREAST 
BEDFAST k 
IS WHEELED BY OTHERS 
22. STAIR CLIMBING (Check ees Complata ea Applicedle) DESCRIBE ANY VEYICES USED 
GOES UP ANN DOWN STAIRS DOES NOT sn STAERS 
W/O HELP — GOES UP/DOWN Sun 
WITH HELP. DESCRIBE ___ GOES UP/DOWN ONE/TWO STEPS 
__USES RAMP FOR ONE/TWO STEPS 
___USES ELEVATOR/CHAIR LIFT 
23. BATHING (Check and Coeplate se Applicable) DESCRIBE ANY DEVICES USED 
BATHES WHERE 
— W/O HELP — BED 
___WITH HELP. DESCRIBE ___ SINK 
—_TUB 
SHOWER 


__-!/$ BATHED bY OTHERS 


DPESSIhH (Check sed Complete es Applicable) 


DRESSES KIND OF ORESS 
___W/O HELP STREET CLOTHES 
— WITH HELP. DESCRIBE .— ROBE L PJ/GOWN 
— SLIPPERS 
___ SHOES 


~. IS DRESSED BY OTWERS 
EATIAG/FEEDING (chach end Complete ae Appliceble) 


EATS WHERE IS FED 
— W/O HELP BEO —_ORALLY BY OTHERS 
— WITH HELP. DESCRIBE __CHAIR IN ROOM JUBE FED 


DINING ROOM PARENTERALLY FED 























26. TOILETING (Check and Completa aa Applicable) 
USES TOILET ROOA WHEN SUBSTITUTES 
—W/0 HELP ALL THE TIME BEOFAN/URINAL 
WITH HELP. DESCRIBE DAY ONLY — COMMODE 
NEVER 
27. BOWEL FUNCTION (Chack and Complete ala Appliceble) 
KO PROBLEM OSTOMY 
~_IMPACTION __SELF CARE 
—_ INVOLUNTARY LOSS MOT SELF CARE 
28. BLADDER FUNCTION (Chack sed Complete en Applicable) 
__KO PROBLEM OSTOMY 
— RETENTION -~ INDWELLIHG CATHETER —_ SELF CARE 
—_ INVOLUNTARY LOSS EXTERNAL DEVICE. HOT SELF CARE 
OE SCRIBE 
29, ORIENTATION (Check Oee enad Complete ss Appiicable) 
—_ ORIENTED IF eT USE THE TABLE TO INOICATE 
__DISORIENTED by (v ) THE AREA(S) AND FREQUENCY, 
COMATOSE 
X 
30, COMMUNICATION OF MEEDS (Chece One 3nd Cumplere es Applicabla) 
—_WERBALLY 
___NONVERBALLY. SPECIFY HOW: __ AS DA > A 
DOES NOT COMMUNICATE 
31. BEHAVIOR PATTERN (chech eed Complete as Applicadla) 
APPROPRIATE BEHAVIOR _ THAPPROPRIATE BEMAVIOR 
— WAHDERING. PASSIVE 
—_ ABUSIVE. AGGRESSIVE 
-— OTHER, DESCRIBE 
JNU OPERATIONS RESEARCH, 6 MAWYS, Juwe, 1972 


— — m o l —™ 





DESCRIBE AMY DEVICES USEO 





OESCRIBE ANY DEVICES USED 


DESCRIBE ANY DEVICES USED 


DESCRIBE ANY CONTROL PROGRAM BEING USEO 


DESCRIBE ANY CONTROL PROGRAM BEING USED 




























ONCE A WEEK OR LESS 
PORE OFTEN THAN OKCE A WEEE 


2 oF & 














He t48 





RESEAKCH GROUP 
IMPAIRMENTS 
32. 


SIGHT (Checb sud Caupleco ss Appticebdle) 
NO IMPAIRMENT 
LEGALLY BLIND 


___IMPAIRMENT, DESCRIBE 2, 


HEARING (Chact ead Complete es Applicedis) 
— KMO IMPAIRMENT 
—— DOES NOT HEAR 
— IMPAIRMENT. DESCRIBE 


SPEECH (Checd aad Completa sa Applicebie) 
—_NO IMPAIRMENT 
DOES MOT SPEAK 
IMPAIRMENT, DESCRIBE 





HIP FRACTURE 
WITH PROSTHESIS 
——MITH REPAIR 


JOINT MOTION (coecth Ove sad Cospiste es Applicable) 
— MO IMPAIRMENT 


JOINT SITE (R. L. B) 


MISSING LIMBS (Checa Ono end Complete ae Applicedie 
— NONE MISSING —MISSING, USE TABLE TO SPEC 


RT... 


MISSING PART(S) 





PARALYSIS (Chech One snd Compinte as Applicable) 
__ NONE 
—__PARACYSIS. DESCRIBE TYPE AND LOCATION 


VENTITION (Chech sod Complete as Applicedie) 
— NG TEETH MISSING 
___ SOME TEETH MISSING 
~~ EDENTULOUS 


MEDICAL STATUS ITEMS _ 
40. RISK FACTOR MEASUREMENTS Specify chu Resdinge ena D 


HEIGHT 
WEIGHT 


BLOOD PRESSURE 
















CIGARETTE SMOKING (Check One end Compirte as Applic 
NEVER SMOKED ~ EX SMOKER 


BE 






ndicated the laluına 


YES 
ALC DHOL ISa 


ANG INAJ MI 
ARTHRITIS 
CARDIAC ARRMYTHA A 
CONGESTIVE WT, FAILUME 


DECUBITUS ULCERS 


HYPERTENSION 


MALIGNANCY 


MENTAL ILLNESS 
WEURDLOGICAL DISORDERS 
RESPIRATORY DISEASE (CHR, ) 





JNU OPERATIONS RESEARCH, & NAWYS, JuvE, 1091 


» FRACTURES AND DISLOCATIONS (coech ead Complate ae Applicable) 
—__ ONE 


__DPISLOCATION(S). DESCRIBE 


— IMPAIRMENT, USE TABLE TD SPECIFY SITE,JOINT AND SIDECS) AND CHECR (9) COLLIN FOR TYPE 











READING 4 DATE 


2. MEDICALLY DEFINED CONDITIONS ei either Bo or Yes Coluen for esth conditina Itated. 


TYPE/LOCATION (AS APPLICABLE) 


DATE OF CLASSIFICATION 


PATIEHT NUMBER 





TYPE OF COMPENSATION 
GLASSES 
CONTACT LENS 
LARGE PRINT 
DIMER. DESCRIBE 





TYPE OF COMPENSATIDN 
JOUD VOICE 

___ SHOUT ING 

HEARING AID 

__LIF READING 

___DTHER, DESCRIBE 


TYPE DF COMPENSATION 
—_ WRITES 

— GESTURES 

—_SICN LANGUAGE 
—JTHER. DESCRIBE 








OTHER FRACTURE(S). DESCRIBE 














) 
IFY MISSING PARTÍSIAND CHECK [WY) COLUMN FOR PROSTHESIS 


LEFT, BDTH 


—JUSCULAR WEAKNESS, WASTING OR ATROPHY 


TYPE DF COMPENSATIDN 
—_PARTIAL PLATE 

— COMPLETE UPPER PLATE 
—_COMPLETE LOWER PLATE 

— OTHER APPLIANCE., DESCRIBE 


USES FOR CHEWING 
YES 


— 


HO 





— SPECIAL DIET. DESCRIBE 





etes in che Spacea Provides Below: 


BLOOD CHOLES. 







READING 





ALBUM | NURE A 


shie? 
PRESENT SMOKER 





SPECIFY NUMBER PER DAY; 
(PAST/PRESENT) 


Also apecify e. 
tios for Zach Cosdlitioco Ip the Last Tus Columss Beieu:) 


DURATION (vas./mos.) 





-. 


304 


278 
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INTERVIEWER 


DATE OF CLASSIFICATION Lo) 
RESEARCH GROUP PATIENT NUMBER 
SPECIAL CARE FACTORS (Check «ed Completa es Applicebla) 


ITEM SPECIFY TYPE RE QUENCY 
(aS APPLICABLE) (LoaY. ImEEK, ETC.) 





— DECUBITI CARE 


— DRESSINGS (STERILE, LARGE, BULKY) 


2 INHALATION THERAPY 
1PPB-O2-NEBULIZER) 
= IAR TGA TION {EXCLUDING BLADDER 
CATHETER) SPECIFY LESION 


OSTOMY CARE 
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